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Chapitre 1 



Théorie des poutres 



L'objectif de ce premier chapitre est de mettre en place et définir toutes les 
notions de base en mécanique des milieux continus permettant d'aborder 
les chapitres suivants traitant de la mécanique des structures, plus commu- 
nément appelée Résistance des Matériaux. 

1.1 PRINCIPES DE BASE EN RÉSISTANCE DES 
MATÉRIAUX 

1.1.1 La notion de contrainte 

Si un solide est en équilibre sous l'action d'un ensemble de forces, de 
couples et de liaisons, ce dernier se déformera. La contrainte est l'objet 
mathématique permettant de quantifier les tensions internes à la matière. 
Pour définir la notion de contrainte, il suffit de procéder par la méthode des 
coupures virtuelles du solide étudié. En un point M, isolons une partie du 
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solide, défini par un plan de coupure orienté par le vecteur normal sortant 
au solide Ht. 




Figure 1.1 

En chaque point M de la surface de coupure, il faut remplacer la par- 
tie du solide manquant par une densité surfacique d'effort sur la coupure 
représentant l'action de ce dernier sur le solide isolé. Cette densité d'effort, 
définie localement en un point M et orientée par une normale sortante ~n 
est appelée le vecteur contrainte T(M, ~rt). Le vecteur contrainte dépend 
linéairement du vecteur unitaire ~rt . Il existe donc localement un opérateur 
linéaire reliant le vecteur contrainte sur un plan à sa normale, c'est le tenseur 
des contraintes cr, symétrique du second ordre. Il vient ainsi, 

~T(M,~n) = a(M).~n 

La matrice du tenseur des contraintes, relative à la base Çe x , ~e y , ~e z ) 
prend la forme suivante : 





&XX 


a xy 


Vxz 


a = 




(Tyy 


Cyz 




_ &ZX 


°zy 


Czz 



Il est usuel de représenter graphiquement le tenseur des contraintes dans le 
plan de Mohr, permettant de séparer les contraintes normales des contraintes 
de cisaillement. Si on désigne par r la contrainte de cisaillement (portée par 
un vecteur t ) et <j n la contrainte normale, on peut décomposer le vecteur 
contrainte en deux contributions : 

T (M, Ht) = cr n ~n + r~t 
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Il existe un repère particulier dans lequel le tenseur des contraintes est dia- 
gonal, c'est le repère principal des contraintes. En notant a\, cr 2 et cr 3 les 3 
contraintes principales, avec <j\ < cr 2 < cr 3 , le domaine d'admissibilité du 
vecteur contrainte est défini par les 3 inéquations suivantes, définissant un 
ensemble de cercles. 



r 2 + (o- n 
r 2 + (a n 



Cr2)((T n 
CT\)((T n 



0-3) > 0 
o- 2 ) > 0 



Dans le plan de Mohr, l'admissibilité de la contrainte est visualisée par 
la zone grisée ci-après. 




Figure 1 .2 



Pour un état plan de contrainte, si on désigne par 6, l'angle entre le repère 
principal des contraintes et le repère dans lequel la matrice du tenseur des 
contraintes est exprimée, la relation entre les contraintes principales et les 
différents termes du tenseur des contraintes est 



0~xx = 



ai + cri a 1 cr\ 



(7 



2 2 
ai + a 2 a 2 — ai 



yy 



2 2 
a 2 - ai 



. cos(2<9) 
. cos(20) 



. sin(2<9) 
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1.1.2 La déformation 

La déformation est la variation relative de longueur d'un solide lorsque ce 
dernier est soumis à une action extérieure. Le tenseur des déformations, sous 
l'hypothèse des Petites Perturbations est la partie symétrique du gradient du 
champ de déplacement U (u xi u y , u z ). 



1 -> 7 
6 = -(grad£/ + grad 



U) 



Dans la base (e x , e y , e z ), le tenseur symétrique du second ordre se cal- 
cule de la manière suivante : 





^xy 




€yx 


e yy 


e yz 


*zx 


e zy 





du x 
dx 

1 (dll x dUy V 

2V dy + ax ' 



a y 

>^ az 



a* 





awj, x 


1 / du x , 
2^ dz + 


a« z 
a* 






2^ az + 


a« z 










a« z \ 
a? ' 


a« z 




az 
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1.1.3 La loi de comportement 

La relation liant le tenseur des contraintes au tenseur des déformations est 
la loi de comportement. Sous l'Hypothèse de Petites Perturbations, pour un 
matériau homogène, linéaire et isotrope, la loi de comportement est un opé- 
rateur linéaire du 4ième ordre. La loi d'élasticité, ou loi de Hooke introduit 
2 paramètres matériaux : le module d'Young E et le coefficient de Poisson 
v. La relation tensorielle s'exprime comme suit : 

l + v v 
e = cr tr cr . Id 

E E 

Id est le tenseur identité, tr cr désigne la trace de cr. À l'inverse, la relation 
exprimant la contrainte à la déformation, en fonction des coefficients de 
Lamé (A et jul) est : 

a = 2/ul€ + Xtre Ad 

Les paramètres d'élasticité et les coefficients de Lamé sont liés par les rela- 
tions : 

A _ V A 

(1 + v){\ - 2v) 
E 

M " 207^0 

En écrivant les tenseurs des contraintes et déformations en vecteur colonne 
de taille 6, on peut exprimer la loi de Hooke de manière matricielle, plus 
simple d'interprétation. 





' 2/x + A 


A 


A 


0 


0 


0 " 






(Tyy 


A 


2jul + A 


A 


0 


0 


0 




e yy 


a zz > = 


A 


A 


2jul + A 


0 


0 


0 


< 


^zz > 


Cyz 


0 


0 


0 




0 


0 


2tyz 


Czx 


0 


0 


0 


0 


fi 


0 




2e Z x 




0 


0 


0 


0 


0 


A* - 




< ^ e xy À 



Dans un cas de chargement bidimensionnel, la loi de comportement se 
retrouve sous forme simplifiée et condensée. Dans un cas d'hypothèse de 
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contraintes planes : 



V.xx 

(Tyy 
O-rv 



(1 - 1/2) 



" 1 


V 


0 " 






1 


0 




V 








1 - V 


{ e yy 


0 


0 


2 J 


{ 2e xy 



Sous l'hypothèse de déformation plane : 



(1 + p)(l - 2v) 



1 - v 
v 

0 



v 

1 - V 



0 
0 

\-2v 



2e Y 



1.1.4 Définitions et hypothèses en mécanique des structures 

La théorie des poutres consiste à associer à la mécanique des milieux conti- 
nus des hypothèses statiques, géométriques et cinématiques permettant de 
réduire la taille du problème à étudier. 



a) La géométrie et le matériau 

Une poutre est un élément particulier de structure décrit par une surface 
plane 2 appelée section droite de centre de gravité G. La ligne moyenne de 
la poutre T est formée par les différentes positions du centre de gravité G de 
la poutre lorsque lorsque l'on parcourt cette dernière selon toute sa longueur 
L. Si 2 est petit devant L alors l'état de contrainte a ainsi que le champ de 
déplacement U du solide pourront être approximés en fonction de quantités 
exprimées uniquement le long de T. 

Si T est une droite, la poutre est dite droite. Si y est dans un plan, la 
poutre est dite plane. Dans les autres cas, la poutre est dite courbe. 

Le matériau constituant la poutre est homogène, isotrope élastique 
linéaire. 
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b) Les hypothèses cinématiques 

Les déplacements, rotations et déformations sont supposés petits. La poutre, 
de part ses équations d'équilibre, sera donc étudiée dans sa configuration de 
référence. C'est l' hypothèses des Petites Perturbations. 

L'hypothèse d'Euler-Bernoulli suppose que, pour toutes transformations 
géométriques, les sections droites d'une poutre restent planes et perpendicu- 
laires à la fibre moyenne. Cette hypothèse permet de définir le déplacement 
U p de tout point P, d'une section droite uniquement en fonction du dépla- 
cement du centre de gravite U G et de la rotation de la section droite © (s) à 
l'abscisse curviligne s. 

U P (s) = U G (s) + ©0) A GP 

Dans le cas d'une poutre droite, l'abscisse curviligne est confondue avec 
l'axe x, le champ de déplacement U (u x ,u y , u z ) de tout point P(x, y, z) peut 
ainsi être obtenu par : 

u x (x, y, z) = u x + 6 y (x).z - O z (x).y 
u y (x,y,z) = u y - 6 x (x).z 
u z (x,y,z) = u z + 6 x (x).y 

Dans la suite, nous définirons la courbure \ P ar X — d®/dx. 
Les déformations non nulles sont données par : 

du x (x) 

e xx (x,y,z) = ^ +Xy(x)-z - XzW-y 



du y (x) 

2e xy (x, y, z) = 0 z (x) - XxW-z 



du z (x) 

2e xz (x,y,z) = + 6 y (x) + x x (x).y 

Sous l'hypothèse que les déformations de cisaillement sont nulles (hypo- 
thèse d'Euler-Bernoulli), on peut définir la rotation 6 d'une section droite 
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y 




Figure 1.4 



par : 



O z (x) 



du y (x) 
dx 

du z (x) 
dx 



si de plus, la poutre ne subit pas de torsion (6 x (x) = 0), alors le tenseur de 
déformations est uniaxial et il reste : 



du x (x) 

e xx {x,y,z) = ^ +x y (x).z - XzM-y 



dx 

2e xy (x,y,z) 0 
2e xz (x,y,z) = 0 



c) Efforts de cohésion et hypothèses statiques 

L'hypothèse de Barré de Saint Venant suppose que pour une section droite 
quelconque, suffisamment éloignée du point d'application des efforts exté- 
rieurs sur une poutre, les effets de ce même chargement peuvent être rem- 
placés par un torseur équivalent s 'appliquant à cette section droite. Pour une 
section droite quelconque, le tenseur des contraintes en un point P prendra 
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la forme suivante : 



a 



xz 



(7 = 



<J 



zx 



0 
0 



0 
0 



Le torseur des contraintes généralisées, par application du principe fon- 
damental de la statique est défini par une composante en efforts R et une 
composante en moments M. 



N est l'effort normal, V y et V z sont respectivement les efforts tranchants 
selon les directions y et z. M x est le moment de torsion, M y et M z sont 
les moments de flexion autour des directions y et z. L'expression de ces 
actions de liaison s'obtient par intégration des contraintes dans la section 
droite. 



1.1.5 Équations d'équilibre d'un élément de poutre 

En isolant un tronçon de poutre, on peut établir les équations différentielles 
d'équilibre d'une poutre en terme de contraintes généralisées. Pour une 



R =Nx+V y y+V z z 
M = M x ~x + M y ~y + M z ~z 



M x = J (<rx Z y ~ (TxyZ)dydz 
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flexion dans le plan (x, y), 3 équations d'équilibre sont obtenues. 

dN(x) 



dx 
dV{x) 

dx 
dM(x) 

dx 



= 0 
-p(x) 
-V(x) 



y 



-V(x) 



p(x) V(x + dx) 



À< -Af(*) ^ ^ ^ ^ ^ À< ^M(x + dx) 



-N(x) 



N(x + dx) 



x x + dx 

Figure 1 .5 



1.2 ÉTUDES DES POUTRES SOUS DIVERSES 
SOLLICITATIONS 

1.2.1 Lois de comportement généralisées pour les poutres 

En mécanique des structures, on distingue certaines actions de base ainsi 
que leurs combinaisons. Par définition, 

- traction ou compression : N ^ 0, V y = V z = 0 et M x = M y = M z = 0 

- flexion pure : M y ^ 0 ou M z ^ 0, N = 0, V y = V z = 0 et M x = 0 

- flexion simple : M y ^ 0 et V z + 0 ou M z ^ 0 et V y + 0, N = 0 et 
M x = 0 



1.2 É tudes des poutres sous diverses sollicita tions 
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- flexion composée : N ^ 0, M y ^ 0 et V z ^ 0 ou M z ^ 0 et V y ^ 0, 
M x = 0 

- flexion déviée : ^ 0, V z ^ 0 et M z ^ 0, ^ 0, Af = 0 et M x = 0 

- flexion composée déviée : N ^ 0, M y / 0, V z ^ 0 et M z ^ 0, Vy 7^ 0, 
M, = 0 

- torsion : M v ^ 0, N = 0, V y = V z = 0 et M y = M z = 0 

Les effets des 3 actions de base sont étudiés dans les paragraphes sui- 
vants. 



a) L'effort normal dans une poutre 

On considère une poutre droite (barre) isostatique de module d'Young E, 
de longueur L et de section S(x) soumise à un effort normal N QXt . L'effort 
normal est constant le long de la barre et vaut N(x) = N ext . La contrainte 
est constante le long de la barre et dans la section droite et vaut : a xx = 
N x (x)/S. Par la loi de comportement, e xx = a xx /E = N(x)/ ES(x). La 
déformation latérale, dans les directions z et y vaut e yy = e zz = —ve xx . 
Considérant la relation entre le déplacement axial et la déformation axiale, 
nous pouvons exprimer la loi de comportement généralisée d'une poutre 
sous un effort axial centré : 

du Y 

N = ES(x)-^ 
ox 

Le déplacement en tout point d'une barrre s'obtient par intégration : 

f L N 
u(x) = / ^7 z dx 



b) Le moment fléchissant dans une poutre 

Considérant une flexion plane selon la direction ~y, l'hypothèse d'Euler- 
Bernoulli permet d'obtenir le profil des déformations axiales le long de la 
section droite : 

d0 7 
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Utilisant la loi de comportement uniaxiale reliant la contrainte et la défor- 
mation ainsi que la définition d'une contrainte généralisée de flexion, nous 
pouvons exprimer le moment fléchissant dans une section 2 d'abscisse x 
par : 



Ex z • y ds 



Soit, en introduisant la définition du moment d'inertie quadratique I z = 
y 2 ds, la loi de comportement généralisée d'une poutre en flexion : 

M z = EI zXz 

Pour deux sections infiniment voisines, sous l'hypothèse des petits déplace- 
ments et petites déformations, on peut assimiler la courbure de la poutre à 
la dérivée de la rotation, soit à la dérivée seconde du déplacement : 



M 7 = EL 



d 2 u y (x) 
dx 2 



La contrainte normale, en fonction du moment fléchissant, est obtenue en 
tout point de la poutre par : 

M z (x)y 



<r xx (x,y) 



L 




€ XX (V) = XzV 

Figure 1.6 



(v) 



M z v 



Les distributions de la déformation et de la contrainte normale le long de 
la section sont linéaires. 



1 .2 Études des poutres sous diverses sollicitations 
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c) Le moment de torsion dans une poutre 

Considérons une poutre de section circulaire (pleine ou é vidée), de rayon R 
soumise à un moment de Torsion constant M x . La contrainte de cisaillement 
en tout point P, repéré par sa distance r au centre de gravité G de la section 
est obtenue par : 

M x r 

T = 

avec / le moment quadratique de torsion (/ = ttR 4 /2, pour une section 
pleine et / = tt(R\ — R\)/2 pour une section évidée de rayon intérieur R2). 



A y 




< >\ 



dx 

Figure 1.7 

Le déplacement relatif de deux sections voisines distantes de dx est une 
rotation dO autour de l'axe x . En introduisant le module de cisaillement du 
matériau G, la loi de comportement généralisée d'une poutre en torsion est : 

d6 M x 
dx GJ 

Pour toute section droite autre que circulaire, l'expression de la contrainte 
et de la loi de comportement généralisée doit prendre en compte le gauchis- 
sement des sections. 

• Section rectangulaire 

En considérant une section rectangulaire de hauteur h et d'épaisseur e (avec 
e < h), la contrainte de cisaillement maximum vaut T max = M x /k\he 2 et la 
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loi de comportement généralisée s'exprime par : 

d6 M Y 



dx Gk2he 3 

Les valeurs de k\ et k 2 peuvent être estimées à l'aide du tableau suivant : 



h/e 


1 


1.5 


2.0 


2.5 


3 


4 


5 


10 


oo 




0.208 


0.231 


0.246 


0.258 


0.267 


0.282 


0.291 


0.312 


1/3 


k 2 


0.141 


0.196 


0.229 


0.249 


0.263 


0.281 


0.291 


0.312 


1/3 



• Profilés minces ouverts 

Pour des poutres dont les sections droites s'apparentent à un I ou un U, 
les contraintes maximales peuvent être approximées de la manière suivante, 
pour les contraintes de cisaillement maximum dans l'âme : 



TU 



Dans l'aile du profilé : 



h\e\ + 2h 2 e\ 

= 3M x e 2 
h\e\ + 2/^2^2 



La loi de comportement généralisée s'exprime par : 



d0 



3M X 



dx G(hie\ + 2h 2 eD 



3 e 2 



<?2 



e 2 



h 2 

Figure 1.8 



1.2 É tudes des poutres sous diverses sollicita tions 
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• Tube mince ouvert 

De manière analogue aux sections rectangulaires étroites, lorsque le rapport 
h/e tend vers l'infini : 

_ 3M X dO 3M X 
max he 2 dx Ghe 3 

• Tube mince fermé 

Pour un tube mince fermé, le flux de cisaillement le long du tube est 
constant, de même donc pour la contrainte de cisaillement. En dénommant 
T la ligne moyenne du tube, et S l'aire délimitée par cette ligne moyenne, 
on peut calculer la contrainte de cisaillement par : 

= Mx_ 
7 2Se 

et la loi de comportement généralisée : 

dO _ M x 
dx GJ 



avec / 



= AS 2 / J ds/e 



1.2.2 Poutre en flexion simple 

En flexion simple, viennent se superposer aux contraintes normales engen- 
drées par le moment fléchissant, les contraintes de cisaillement liées à la 
présence d'un effort tranchant non-nul à l'abscisse x de la poutre. Pour 
une section droite quelconque, de moment quadratique I z , la contrainte de 
cisaillement, en tout point y de cette section s'obtient en figure 1.9. 
b(y) est la largeur de la poutre à la côte y et H(y) est le moment statique de 
la section au dessus de la côte y. Pour une section rectangulaire (largeur b et 
hauteur h), le profil des contraintes de cisaillement est parabolique, nul aux 
bords et maximum au centre : 

*>-^0-£) 
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r(y) = 



V y (x)H{y) 
hHy) 



Figure 1 .9 



Le diagramme des contraintes reste linéaire dans la section : (r xx (x, y) = 
-M z (x)y/I z . 

1.2.3 Poutre en flexion déviée 

Lorsque l'axe du moment fléchissant extérieur ne correspond pas avec l'un 
des axes principaux d'inertie de la section, il est nécessaire de décomposer 
ce dernier en deux composantes, selon les axes principaux d'inertie de la 
section afin de procéder au calcul des contraintes par superposition. Ainsi, 
la contrainte normale pourra être obtenue par : 



Si une section droite est repérée par un sytème de coordonnées x'y'z, Les 
axes principaux d'inertie d'une section droite portée par la normale ~x' sont 
obtenues par une rotation d'un angle a autour de x' tel que : 



1.2.4 Poutre en flexion composée 

Lorsqu'aux actions de flexion se rajoute une composante d'effort normal, 
nous parlerons de flexion composée. La contrainte normale est obtenue par 



o-xxM 



M y (x)z _ M z (x)y 
Iy h 



tan 2a = 




1.2 É tudes des poutres sous diverses sollicita tions 
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superposition des différents termes : 

N M y (x)z M z (x)y 

O-xxW = — + — 

O ly 1 Z 

Dans certains cas, il peut être intéressant de remarquer que le sytème d'effort 
est équivalent à un effort unique appliqué dans le plan de la section appelé 
le « centre de pression ». La position de ce centre de pression est y\ = 
—M z /Netzi = My/N. Sile centre de pression est à l'intérieur « du noyau 
central » de la section, alors toutes les contraintes normales sont de même 
signe, quelque soit le point de la section droite considéré. À l'inverse, si 
le centre de pression est en dehors du noyau central, alors les contraintes 
normales de part et d'autres de la section droite seront de signes contraires. 
La limite du noyau central est obtenue en faisant coïncider l'axe neutre de 
la poutre avec le contour de la section droite. 

Prenons l'exemple d'une poutre rectangulaire ne subissant qu'un moment 
de flexion M z et un effort normal N. La position du centre de pression repé- 
rée par ses coordonnées z\ — 0 et y\ = —M z /N permet de calculer la 
contrainte en tout point, sachant que I z = bh 3 / 12 : 

N_ 6M Z 
Œxx ~bh + ~bh^ 

Pour les bords les plus sollicités : 

= f ( 1+6 t) 

<r xx (- yi )=j (l-6|) 

La frontière du noyau central, dans ce cas, est obtenue en imposant <j xx = 0 
en y = y\ et y = — yi, soit h/6 et — h/ 6. La même étude peut être effectuée 
pour une flexion selon z due à un moment M y permettant de définir pour 
limite du noyau central b/6 et — b/6. La forme du noyau central pour une 
section rectangulaire est un losange centré sur G, de dimension h/3 et b/3. 



Chapitre 2 

Caractéristiques des sections 



2.1 PRÉAMBULE 

Les calculs des différentes intégrales se font dans un repère cartésien ou 
dans un repère cylindrique. 





Coordonnées du point 


Surface infinitésimale 
ds 


Repère cartésien 


<y,z) 


dydz 


Repère cylindrique 


(r, 0) 


rdrdd 



z = r cos 6 
y = r sin 6 

2 2 2 

r 1 = y l + z 




2.2 Définitions 



2.2 DEFINITIONS 

2.2.1 Surface 

Surface de la section (2) : 




Figure 2.2 

2.2.2 Centre de gravité 

Position du centre de gravité de la section (2) par rapport à l'axe Oy 



ff%y ds = noz 

JLds S 



[m] 



Position du centre de gravité de la section (2) par rapport à l'axe Oz 

[m] 



zg 



Sk zds n 



Oy 



Sh dS 



2.2.3 Moment statique 

Moment statique de la section (2) par rapport à l'axe Oy : 



Oy 



-IL 



zds 



[m 3 ] 
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Moment statique de la section (2) par rapport à l'axe Oz : 



H 



Oz 



IL yds 



[m 3 ] 



Propriétés 

- Le moment statique d'une surface est nul par rapport à un axe passant 
par son centre de gravité. 

- Le moment statique d'une surface est égal au produit de sa surface par 
la distance de son centre de gravité : 

Ho y = S zg 7~Loz = S yc 

2.2.4 Moment d'inertie 

Moment d'inertie (ou moment quadratique) de la section (2) par rapport à 
l'axe Oy : 



[m 4 ] 



Moment d'inertie de la section (2) par rapport à l'axe Oz : 



Oz 



IL 



y 2 ds 



[m 4 ] 



Propriétés 

- Les moments d'inertie sont toujours positifs. 

- Les moments d'inertie sont minimaux au centre de gravité de la section 
(voir le théorème d'Huygens en 2.3.1 pour l'explication). 



2.2.5 Produit d'inertie 

Produit d'inertie de la section (2) par rapport aux axes Oy et Oz : 



l Oyx 



IL yzds 



[m 4 ] 



2.2 Définitions 
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2.2.6 Moment polaire 

Moment polaire de la section (2) par rapport à l'axe Ox : 



•°=IL 



r z ds 



[m 4 ] 



Propriétés 

- Le moment polaire est toujours positif. 

- Relation entre le moment polaire et les moments quadratiques : 



lo 



JJ^ 2ds = JJ^y 2 + z 2 )ds = i 0y + ioz 



2.2.7 Axes principaux d'inertie 

Le repère formé par les axes principaux d'inertie GY et GZ est tel que : 

- Le produit d'inertie est nul dans ce repère. 

- Les moments d'inertie sont respectivement maximal par rapport à un 
axe et minimal par rapport à l'autre. Ces moments sont alors appelés 
moments principaux d'inertie. 

Propriétés 

- Si la section possède un axe de symétrie, cet axe est un axe principal 
d y inertie. 

- Si la section possède deux axes de symétrie, ces axes sont des axes prin- 
cipaux d'inertie. 

2.2.8 Rayon de giration 

Rayon de giration d'une surface (2) par rapport à l'axe Oy : 
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Rayon de giration d'une surface (2) par rapport à l'axe Oz : 




2.3 THÉORÈMES ET PROPRIÉTÉS 

2.3.1 Théorème de Huygens 

Connaissant les inerties d'une surface dans un repère (G,y,z), le théorème 
de Huygens permet d'exprimer les inerties dans un repère d'axes parallèles 
(0,y,z). La réciproque est naturellement vraie. 

loy = iGy + S Zq 

ioz = iGz + Syl 

hyz = iGyz + $ y G ZG 

2.3.2 Changement de repère 

a) Repère orienté par un angle a 

On considère le repère (0£i, Oy\) orienté par un angle a par rapport au 
repère (0£, Oy). Les coordonnées d'un point dans ce nouveau repère sont : 

yi = — sin(a)z + cos(a)y 

Z\ = cos(a)z + sin(a)y 

Les moments d'inertie d'une section dans le nouveau repère s'écrivent : 

iGyx = \Vgz + hy) + \^Gz ~ lGy)C0S(2a) + I Gyz sin(2ûf) 

kzi = \( j g z + h y ) ~ \^ l Gz ~ lG y )cos(2a) - I Gyz sin(2a) 
Le produit d'inertie est : 

I GyiZl = \{I Gz ~ J Gy) sin(2û0 + I Gyz cos(2a) 



2.3 Th éorèmes et propriétés 
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b) Repère principal d'inertie 



L'angle a entre le repère principal d'inertie OZ et 07 et le repère (0£, Oy) 
est : 



tan(2a) 



27, 



Gyz 



iGy — iGz 

Les moments d'inertie principaux, respectivement maximal et minimal, 
sont : 

1 



ÏGy + Y/ {Igz ~ iGy) + 47, 



Gyz 



7min = - ( Igz V (7gz — iGy) + 47 G> , Z 




Figure 2.3 



2.3.3 Décomposition d'une surface 

Une section (2) peut être décomposée en n sections élémentaires (Si), i G 
{1, ..,n}. Les caractéristiques géométriques de la section (2) peuvent être 
calculées à partir de la somme algébrique des caractéristiques géométriques 
des sections (Si), en utilisant si nécessaire le théorème de Huygens pour le 
transport des inerties. 

Exemple. Calcul des caractéristiques d'une cornière en L 



(2) = (5 1 ) + (<S 2 )-(<S 3 ) 
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e 

.G 


h 


0 


> 

e 








►1 

oi 




o 






0 




(2) 




(<$i) 


(5 2 ) 


(S3) 


Section 


e(h + b — e) 




be 


e 2 


Moment 
statique 
lOz 


|(/7 2 + be-e 2 ) 


eh 2 
~2~ 


be 2 
2 


e 3 
2 


Moment 
d'inertie 
lOz 


|(/7 3 +be 2 -e 3 ) 


e/7 3 


be 3 


e 4 
3 



Tableau 2.1 



Coordonnées du centre de gravité de la surface (2) : 



y? 



K { oz + W ( £ 2) - H^f _h 2 + be-e 2 



S<S) 



2(h + b- e) 



7 (X) 
Z G 



n Oy + ^Oy 
S® 



<iy(«S 3 ) 



he + b 2 - e 2 
2(h + b- e) 
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2.4 CARACTÉRISTIQUES DES PRINCIPALES 
SECTIONS 



Géométrie 


Section 


Centre de 
gravité 


Moment statique 


Moment quadratique 


0 


•G 

a 




S = a 2 


a 

ZS= 2 
a 

ys= 2 


ri-Oz = ~Y 

n 0y = — 


, a 4 
/Gz= Î2 

, * 4 


0 


*G 


11' 


S = A 2 -a 2 


A 

Zo= 2 
A 

yo = j 


A 3 -a 2 A 

rtOz — ^ 

A* -a 2 A 

n °y= 12 


A 4 - a 4 

lGz= 12 

A" -a" 
'Gy= 12 


0 


•G 

b 


I" 

± 


S = bh 


b 

Z *=2 
h 

y <=2 


ni bh 2 

ri-Oz = ~Y~ 

n 0y = — 


bh 3 
/Gz = l2" 
hb 3 


0 


G 

< B > 


h H 


S = BH-bh 


B 

ZG= 2 
H 

ys= 2 


BH 2 - bhH 

rioz = ^ 

B 2 H - bhB 

n 0y = 2 


BH 3 - bh 3 
l Gz = 12 

HB 3 - hb 3 
<Gy- 12 




G 


h H 


S = BH-h(B-b) 


B 

Zs= 2 
H 


BH 2 - hH{B - b) 

rLOz — ^ 

B 2 H - hB{B - b) 


fiH 3 - h 3 (B - b) 
Igz- 12 

B 3 H - h(B 3 - b 3 ) 










n 0y = 2 


'Gy= 12 






S = BH-h(B-b) 


B 

ZC= 2 

y e = ^f 


BH 2 - h 2 (B - b) 
rtoz = 2 

B 2 H - hB{B - b) 


/6z=fiW 3-,3 (fi _ b) yis 

fi 3 H - MB 3 - b 3 ) 


0 




G 

? , 


h H 








n 0y = 2 


>Gy= 12 
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Géométrie 


Section 


Centre de 
gravité 


Moment statique 


Moment quadratique 


0 


t b 'h 
.G fi 


/ H 


S = BH-h(B-b) 


z s = ^f- 


(B-b)h' 2 +bH 2 
rtoz = 2 

B 2 h' + hb 2 

H 0 y = 2 












S = BH-h(B-b) 


H 




BH 3 - h 3 (fi - b) 
/Gz= 12 

2/ib 3 + (H — h)(2B - b) 2 
= 24 


0 




• G / 


1 H 






0 


< b 1 


} 


S= ^bh 


b 

Za = ï 
h 


-H ^ 2 


bb 3 
/Gz = ^6" 
hb 3 


° 

1 


3 

d 




s - 770(2 

4 


d 

Ze= 2 
d 

y e = j 


/w Td 3 

-H 770/3 


ttD 4 
/Gz= "64" 

ttD 4 
/( * =_ 64" 


1 


D 


"à 


s _ 77(D 2 -d 2 ) 

4 


D 

Ze= 2 
D 

VG =2 


tt(D 3 - d 2 D) 
H 0 z= g 

77(D 3 - d 2 D) 


77(D 4 - d 4 ) 

' Gz = 64 

77(D 4 - d 4 ) 

/G ' = 64 








^ 0 


d 


\ 


s 770/2 

~ 8 


d 

ZG= 2 
2d 

Vg= =- 

37T 


d 3 

Hoz = Î2 
^ = 16" 


d 4 /7T 8 \ 
GZ "16U Stt) 

77d 4 

Gy " 128 



2.5 Exemple : caractéristiques d'une section en T 
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Centre de 
gravité 



Moment statique 



Moment quadratique 



^ o 



7TA- 2 



3 77 

4r 



l G z = 





(tt 8 




V8 9tt 




^77 8 




V8 977 



z l *GJ 

O 



1~Loz - 



7-Loy 



ab 2 



ira 2 b 



l G z- 



irab 3 




S = 2eh- e l 



T~(-Oz '- 



2eh 2 - e 2 h 



2eh 2 - e 2 h 



Igz-- 



eh 3 + e 3 h — e 4 
12 

eh 3 + e 3 h — e 4 




y g 



4 

h 3 



l Gz = 

iGy = 



O 
Oh 

J 

I 



2.5 EXEMPLE : CARACTÉRISTIQUES D'UNE SECTION 
ENT 

Si- 



o 



G 



zg 



>S 2 



e = 10 mm 
if = 120 mm 
L = 120 mm 
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1. Surface : 



S = S (Sl) + S (S2) = e L + e (H - e) 
= 10 x 120 + 10 x (120 - 10) = 2300 mm 2 

2. Moments statiques : 



rZ =L />y=// rz=L/2+e/2 py=H-e 

z dy dz+ I I z dy dz 

Jz=L/2-e/2 Jy=0 



rz=L ry=H 
J z =0 Jy=H- 



120 



[y] 



120 
110 ' 



J 0 



65 



[y] l o l ° = 138 000 mm 3 



55 



H 0z = H ( £ + H ( £ 

rz=L ry—H nz—L/2+e/2 r-y—H—e 

= / y dy dz+ / y dy dz 

J z =0 Jy=H-e Jz=L/2-e/2 J y=0 



Mo 



120 



2 



-, 120 



110 



2 



-, 110 



198 500 mm 3 



3. Centre de gravité : 

7~Loz 
y g = — ^r- 



86, 3 mm 



Zg = = 60 mm 



4. Moments d'inertie : 
- Au point O : 

lOy — In,, + I 



Oy ^ 1 Oy 

çz—L ry=H 

J Z =0 Jy = H-i 



z dy dz + 



pZ=L/2+e/2 ry=H-e 

/ / z 1 dy dz 

Jz=L/2-e/2 Jy=0 



120 



[y] 



120 
no • 



j 0 



65 



M™ = 9, 73 x 10 6 mm 4 



J 55 



Exemple : caractéristiques d'une section en T 
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L 0z ~ *0z + L 0z 

r z=L r y=H 



r>z,=L py—n pz—L/2+e/2 r-y—H—e 

/ / y 2 dy dz+ \ \ y 2 dy dz 

Jz=0 Jy=H-e Jz=L/2-e/2 J y=0 



120 



= lz] 0 

Au point G : 



120 



no 



no 



20, 3 x 10 6 mm 4 



J 0 



I Gy = Io y ~ z 2 G S = 1, 45 x 10 6 mm 4 
Igz = loz ~ y 2 G$ = 3, 19 x 10 6 mm 4 



Produit d'inertie : 



t _ t(Si) , 

7/Ov7 — + L 



Oyz — 1 Oyz + 1 Oyz 

rz=L ry—H nz—L/2+e/2 py—H—e 

= / yz dy dz+ / yz dy dz 

J z =0 J y =H-e Jz=L/2-e/2 J y=0 



V 


120 


\y 2 } 


120 


V 


65 




_2_ 


0 


_ 2 _ 


+ 

no 


_2_ 


55 


_ 2 _ 



z=L/2-e/2 Jy=0 
110 



11,9 x 10 6 mm 4 



Rayon de giration : 



io y = y -^7- = 65, 0 mm 



l Oz 




94, 0 mm 



Chapitre 3 



Théorèmes Généraux-Méthodes 
énergétiques 



3.1 PRINCIPE DES TRAVAUX VIRTUELS - PTV 

Soit un système de poutres en équilibre sollicité par des actions extérieures 
et tel que les liaisons soient parfaites. 
Actions extérieures : 

- forces de surface (ponctuelles ou réparties), 

- forces à distance, 

- forces de liaison (à priori inconnues). 
On s'intéresse aux : 

- efforts de liaison, 

- déplacements, 

- déformations, 

- contraintes. 



3.1 Principe des travaux virtuels - PTV 
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Restrictions : 

- pas d'effet de la température, 

- état initial non sollicité. 

Pour déterminer l'état du système, on peut utiliser le premier principe de 
la thermodynamique (structure élastique, pas de dissipation). 

Dans un référentiel galiléen et pour toute partition du système considéré, 
la somme du travail des forces extérieures T ext et du travail des forces inté- 
rieures Tim calculée dans un champ de déplacement virtuel quelconque est 
égale à zéro. 

Wp ext + Wr int = 0, Vw*virtuel (3.1) 

Pour utiliser ce principe, il faut : 

- définir ce que signifie u* 

- expliciter le W£ 

3.1.1 Champ de déplacement virtuel 

A priori, u* est quelconque mais dans la pratique on utilisera uniquement : 
Champ cinématiquement admissible : champ de déplacement virtuel qui 

respecte les conditions de liaisons, u* C.A. est continu, dérivable et tel que 

la déformée existe. 

Champ rigidifiant : champ de déplacement virtuel tel que la déformée est 

nulle en tout point de la structure. 

3.1.2 Définition du travail des forces dans le champ de 
déplacement virtuel 

Travail des forces extérieures : si on note F t , M t les forces et moments 
appliqués au point G t et ïï{, co* les déplacements et rotations de la section 
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en G/, nous obtenons alors : 

n 

w>„ = £{3*r+^(3r} 0.2) 

i=\ 

Travail des forces intérieures : lié au travail de a dans le champ de défor- 
mation virtuel : e* déformation axiale, \* courbure, y* déformation de 
cisaillement, 0* angle de torsion. 

W* rint = - I (Ne* + M\* + Vy* + M x 0*) (3.3) 

J structure 



3.2 ÉGALITÉ DE CLAPEYRON 

Le principe des travaux virtuels (PTV) est valable quelque soit w*, en parti- 
culier on peut choisir pour champ virtuel le champ réel de déplacement. On 
obtient alors : 

W^+W^ =0 (3.4) 

yV^ ext : les efforts et les déplacements ne sont pas indépendants, il faut 
donc prendre en compte l'histoire du chargement. 

n 1 f 1 

y^T ext = Y, 2^ iiïi + + j i^ Uds (3 - 5) 

—Wjr int : potentiel interne W ou énergie de déformation élastique. 

W W f /N 2 M 2 V 2 M 2 \ ^ 

Jstructure^ES EI GSr GjJ 

avec S R : section réduite en cisaillement. 

L'égalité de Clapeyron ne dépend que de l'état initial et de l'état final de 
la structure. 



3.3 Théorème de réciprocité de Maxwell-Betti 
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3.3 THÉORÈME DE RÉCIPROCITÉ 
DE MAXWELL-BETTI 

Ce théorème se déduit de l'égalité de Clapeyron en considérant deux sys- 
tèmes de chargement appliqués à la même structure. 

Le travail produit par un système de chargement S\ sur une structure dans 
le champ de déplacement dû à un système de chargement S2 est égal au 
travail du système de chargement S 2 dans le champ de déplacement dû à S\ . 

Cette égalité est donnée en produit scalaire. On ne pourra donc connaître 
que la composante du déplacement suivant la direction d'application de la 
force. 

Le déplacement (ou la rotation) produit en / par une force (ou couple) 
unitaire agissant en j est égal au déplacement (ou la rotation) produit en j 
par la force (ou couple) unitaire agissant en i. 



3.4 THEOREME DE CASTIGLIANO 

Le théorème de Castigliano établit une relation entre les déplacements et 
le potentiel interne. Pour l'établir, on part de l'égalité de Clapeyron et on 
calcule plus explicitement le travail des forces extérieures. 

La dérivée partielle du potentiel interne par rapport à une action quel- 
conque est égale au déplacement du point d'application de cette action 
mesurée algébriquement sur la ligne d'action de celui-ci. 

Pour une force ponctuelle F*, le déplacement 8 k est ainsi : 

wr 8k (3 - 7) 

Pour un moment ponctuel M k , la rotation (o k est ainsi : 

dW 

wtt- = <*>k (3.8) 
dM k 



En conséquence, si l'on souhaite calculer le déplacement 8 (ou rotation 
co) d'une section 2 d'une poutre dans une direction donnée, on applique 
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une force fictive F* (ou moment fictif M*) dans la section 2 suivant cette 
direction. On aura alors : 

dW 



ô = 



dF* 

dW 
<9M* 



F*=0 



M*=0 



(3.9) 



3.5 THEOREME DE MENABREA 

Dans un système hyperstatique sur appui rigide, les réactions hyperstatiques 
dues aux liaisons surabondantes ne travaillent pas pendant la déformation 
du système. Les dérivées partielles du potentiel par rapport aux réactions 
hyperstatiques R t sont donc nulles : 

dW 

wr 0 (3 - 10) 



3.6 THEOREME DE MULLER-BRESLAU : 
FORMULE DE MOHR 

Pour calculer le déplacement (ou la rotation) en un point P d'une structure 
(S) isostatique sollicitée en flexion par un système de chargement extérieur, 
on procède de la manière suivante : 

- on calcule les efforts intérieurs associés au chargement extérieur (Le. 
M, N, V). 

- on applique une charge unitaire l P au point P dans le sens du 
déplacement recherché et on détermine les efforts intérieurs associés 
(M, N, V). 

Le déplacement (ou la rotation) recherché ô P s'obtient par l'intégrale 
suivante : 

f (MM NN VV \ , 



3.6 Théorème de Mûller-Breslau : Formule de Mohr 
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Si le système est hyperstatique, on démontre que les efforts intérieurs 
M, TV et V peuvent être calculés sur n'importe quelle structure isostatique 
déduite de S en supprimant les liaisons hyperstatiques. 

Dans la plupart des cas, les deux derniers termes de l'intégrale peuvent 
être négligés et l'on obtient : 



> Calcul des intégrales de Mohr par une méthode simplifiée 

Il s'agit d'une méthode simple lorsqu'un des diagrammes est linéaire (avec 
El constant), voir la figure 3.1. 




(3.12) 



Mi(x) 




B : aire sous la courbe 



x l lx G; 



x 



M,(x) 



A 




x l x Gj 



X 



Figure 3.1 



On obtient alors : 




(3.13) 



avec A la valeur dans le diagramme Mj(x) (linéaire) au niveau du centre de 
gravité x^ du diagramme M t (x) et B l'aire sous la courbe M t (x). 
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3.7 LIGNES D'INFLUENCE 



Considérons une poutre à plan moyen chargée dans son plan par des forces 
verticales. On appelle fonction d'influence f(a) l'expression en une section 
droite 2 donnée d'un effet (noté E) donné, engendré par une cause de valeur 
unité se déplaçant sur la poutre. 
Causes : 

- force ponctuelle, ensemble de forces ponctuelles, 

- couple ponctuel, 

- déplacements imposés. 



- déformation (déplacement), 

- actions de liaison, 

- efforts internes (N,V, M, M x ). 

La représentation graphique de la fonction d'influence est appelée ligne 
d'influence de l'effet E dans la section 2. 

Soit une poutre sur deux appuis soumise à une charge ponctuelle 



Effets : 



(figure 3.2) : 



Y 



X 




a 




Figure 3.2 



3.7 Lignes d'influence 
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Les fonctions d'influence de V et de M en une section 2 en x sont : 

a-l 



x < a 



x > a 



V(x,a)-. 

M(x. a) = x 
a 



l 



l-a 



(3.14) 



V(x,a)-. 



l 



M(x. a) = a 



l 



Notons que, dans ces expressions, la variable est a tandis que x est fixée. 
Nous obtenons les lignes d'influences suivantes (figure 3.3) : 

V(a) 
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3.7.1 Effet d'un ensemble de charges 

L'effet E d'un système de charges concentrées P\, P 2 . . . P n d'abscisses a\, 
ct2 • • • ol u est égal, en vertu du principe de superposition, à : 



L'effet d'une charge répartie p(a) entre deux points A et B est égal à : 



3.7.2 Lignes d'influence des déformations 

Le théorème de Maxwell-Betti montre que : 

- la ligne d'influence du déplacement vertical d'une section 2 d'abs- 
cisse x sous l'action d'une charge verticale unité d'abscisse a est la 
ligne représentative d'un déplacement vertical de la section 2' d'abs- 
cisse a sous l'action d'une force verticale unité placée dans la section 
d'abscisse x ; 

- la ligne d'influence de la rotation d'une section 2 d'abscisse x sous 
l'action d'une charge verticale unité d'abscisse a est la ligne repré- 
sentative d'un déplacement vertical de la section 2' d'abscisse a sous 
l'action d'un couple unité placé dans la section d'abscisse x. 



E = P x f{a x ) + P 2 f(a 2 ) + . . . + Pnfia n ) 



(3.15) 




(3.16) 



Chapitre 4 

Systèmes isostatiques 



L'objectif de ce chapitre est de donner les solutions de base en terme de 
réactions d'appuis, diagramme de contraintes généralisées et déplacement, 
pour des systèmes de poutres isostatiques. 

4.1 DÉFINITIONS 

4.1.1 Systèmes isostatiq ues 

Un système est dit isostatique si toutes les réactions d'appuis ainsi que les 
contraintes généralisées en tout point de la structure peuvent être détermi- 
nées par la seule utilisation des équations d'équilibre. Les caractéristiques 
géométriques de la structure (inertie, section, matériau) n'entrent à aucun 
moment dans le calcul de la distribution des efforts. Ces dernières ne 
seront introduites que dans les chapitres suivants dans l'étude des systèmes 
hyperstatiques. 
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4.1.2 Efforts et conditions de liaisons 

Les structures étudiées ci-après peuvent être soumises à des efforts et 
couples ponctuels (F et C) ainsi qu'à des densités d'efforts Çp), appelées 
aussi charges linéiques. 

Nous ne considérerons dans ce chapitre, dans un souci de simplicité 
d'écriture, que des structures planes dans le plan (x ,~y). Les rotations et 
moments fléchissants sont donc définis autour de l'axe ~z. La structure est 
en équilibre sous l'action de ces efforts et de conditions de liaisons. Les prin- 
cipales liaisons entre la structure et son environnement sont résumées dans 
le tableau ci-dessous. Nous donnons les réactions d'appuis produites ainsi 
que les degrés de liberté en déplacement bloqués ou autorisés. Les efforts 
verticaux sont notés V, les horizontaux H et les moments fléchissants M. 
Les déplacements axiaux sont notés u, transversaux v et les rotations 6. 



Type de liaisons 


Schémas 


Réactions d'appuis 


Degrés de liberté 


Rotule 




H / 0, V ^ 0, 
M = 0 


u = 0, v = 0, 0 ^ 0 


Appui simple 




V / 0, H = 0, 
M = 0 


v = 0,u^0, 0/0 


Encastrement 


y\ 
y 
y 
y 
y 
y 


V / 0, H ^ 0, 
M/ 0 


v = 0,u = 0,6 = 0 



4.1.3 Exemple 

Un exemple simple est traité ici, permettant dans un cadre isostatique de 
comprendre les principales étapes de détermination des efforts et déplace- 
ments dans la structure. Soit une poutre isostique de longueur /, de moment 
quadratique / et de module d' Young E supportant une charge concentrée au 
tiers de sa longueur comme indiqué sur la figure ci dessous. 



4.1 Définitions 
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1/3 



Figure 4.1 



> Détermination des réactions d'appuis 

Nous noterons H les réactions d'appuis verticales en un point donné, V les 
verticales et M les moments fléchissants. Les équations d'équilibre en effort 

^ ^ F/ x et ^ F/y^j nous indiquent que : 

H A = H B =0 
V A + V B -P=0 

L'équation d'équilibre en moment, écrite au point A nous permet de 
déterminer les réactions verticales en A et B. 

v - 2P 

V B - — 
P 

> Détermination des contraintes généralisées 

Dans la détermination le long de la poutre des contraintes généralisées 
(effort normal, effort tranchant et moment fléchissant), nous appliquerons 
la méthode des coupures. Par application du principe fondamental de la 
statique, cette méthode permet de calculer le torseur des efforts intérieurs 
en tout point d'une section droite. Pour une section comprise entre le point 
A et le point d'application de l'effort extérieur P : 

l'équilibre en effort de la section d'abscisse x, selon ~x nous indique que : 



N(x) = 0 
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V(x) 



Af(jc) 



Figure 4.2 

l'équilibre en effort selon ~y nous permet de déterminer l'effort tran- 
chant : 

V(x) = -V A = -j 

l'équilibre en moment autour de ~z nous permet de déterminer le moment 
fléchissant : 

Px 

M(x) = V A x = — 

De même pour les sections de poutre situées entre le point d'application de 
P et B. 

N(x) = 0 



V(x) 



V A + P 



2P 



M(x) = V A x-p[x- j y — - — 

Les diagrammes des efforts tranchants et moments fléchissants ainsi calcu- 
lés sont tracés sur la figure ci-après : 



V(x) 



M(x) 



2P/3 




-P/3 



Figure 4.3 



4.2 Poutre sur deux appuis 
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> Calcul des déplacements 

La flèche de la poutre est calculée en utilsant la loi de comportement reliant 
la courbure (dérivée seconde de la flèche : v"(x)) au moment fléchissant : 
M(x) = EIv"{x). Ainsi pour la partie de la poutre entre le point A et le 
point d'application de P, 

v(x) 



px 
ÏSË1 



Pour la partie restante de la poutre : 

p(l - x) 



v(x) 



1SEI 



x(2l — x) 



4^ 
~9~ 



4.2 POUTRE SUR DEUX APPUIS 

Dans ce paragraphe, nous nous proposons de passer en revue les différents 
cas de charges possibles rencontrés pour une poutre simplement appuyée 
sur ses appuis. 

4.2.1 Cas d'une charge concentrée 

i a | b | Pb Pa 

V A = — 0a = -—77^ - a)W ~ à) 



1 



p i 6Eir 



A ^~i^__^^ B Pa n Pa , 2 2N 



Figure 4.4 

> Calcul des contraintes généralisées 

Pour 0 < x < a, 

Pb 
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Pour a < x < /, 




M(x) = Pa 




Le moment maximum est en a et vaut M, 



Pab 



Max 



/ 



> Calcul de la déformée 



Pour 0 < x < a, 



v(x) 



Pbx 
6EII 



[l 2 -b 2 -x 2 ] 



Pour a < x < /, 



v(x) 



Pa(l - x) 
6EII 



[x(2l 



x) — a 2 



Le déplacement au centre de la poutre en x =1/2 vaut : 

Pa(3l 2 -4a) 
v(l 2) = ^ J - 

4.2.2 Cas d'un convoi de charges ponctuelles : 
théorème de Barré 

Un convoi est un système de charges ponctuelles susceptibles de se dépla- 
cer dans son ensemble le long de la poutre. Si la détermination des efforts 
tranchants maximaux aux appuis de la poutre ne pose pas de problème, la 
localisation du moment fléchissant maximum, fonction de la position du 
convoi sur la poutre doit être menée avec précision. Le théorème de Barré 
permet de déterminer la position du convoi fournissant le moment fléchis- 
sant maximal : le moment fléchissant est maximum au droit d'une charge 
lorsque cette charge et la résultante R des charges du convoi se trouvent 
dans des sections symétriques par rapport au milieu de la poutre. 



4.2 Poutre sur deux appuis 
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pi, p: 



P3 



I / 1 

Figure 4.5 

Prenons l'exemple de deux charges sur une poutre, distante d'une lon- 
gueur a entre elles, le moment maximal est égal à : 



4.2.3 Cas d'une charge uniformément répartie 

L'effort p exprimé en N /m est uniformément réparti sur toute la longueur 
de la poutre. 



v B = p - e B = ^- 

2 24EI 

Figure 4.6 



Calcul des contraintes généralisées 

V(x) = -p ' 1 
M(x) 



& \2 

px(l — x) 



2 



© Le moment maximum est en / /2 et vaut M Max 
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> Calcul de la déformée 

vM = -^j[l 3 -Vx 2 +x 3 ] 
Le déplacement est maximal au centre de la poutre en x =1/2 et vaut : 

w(Z/ 2) = -J£^ 
v ' ' 384£7 

4.2.4 Cas d'une charge répartie partielle 

f§ — g* Vb = p^^> ^ + .) 2 _ 2 _^ 2 

u = ► 2/ 48£/Z 



/ 



Figure 4.7 



> Calcul des contraintes généralisées 

Pour 0 < x < a, 

V(x) = + 2c) M(x) = ^-(b + 2c)x 

Pour a < x < b + a, 

V(x) = -^-(b + 2c) + p(x-a) 

pb (x — a) 2 

M{x)=^{b + 2c)x-p K —^ 

Pour b + a < x < l, 

pb 

V(X)= 2Ï (2a + b) 
M(x) = ^(2a + b)(l-x) 



4.2 Poutre sur deux appuis 
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b(b + 2e) 

Le moment maximum est obtenu pour x = a + — et vaut : 

M Max = ^ + 2c ) [ b ( b + 2c ) + 4al ] 
> Calcul de la déformée 

La flèche maximale est obtenue entre a et a + b et s'exprime par : 

u(jc) = ^— \b(b + 2c)x [4(/ 2 - x 2 ) - (b + 2c) 2 - Z? 2 1 + 2/(jt - a) 4 ] 

AxEIl J J 

4.2.5 Cas d'une charge répartie partielle proche d'un appui 

n pa / a\ pa 2 (2l — a) 2 



24E/Z 



™2 ^2/ 9 /2 _ „2 



pa 1 pa z (2l z - a 1 ) 



21 2AEII 
Figure 4.8 

2 > Calcul des contraintes généralisées 

a Pour 0 < x < a, 

S V(x) = -^[a(2/-a)-2/x] 



% Pour a < x < l, 



M(x) = ^- [a(2/ - a) - /x] 



MCI = , "" !< ' " *> 



o 

-a - - ' 2/ 

i Le moment maximum est obtenu pour x = — (21 — a) et vaut : 
-g 2/ 

| ^ 2 (2/-a) 2 

© m Max — 
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> Calcul de la déformée 

Pour 0 < x < a, 



v(x) = - ^ X [a 2 (2l - a) 2 - 2a(2l - a)x 2 + lx 3 ] 



2AEII 

Pour a < x < /, 



= _ pa (/ x) , _ 2 _ 2 . 
v 7 24£7Z L J 



4.2.6 Cas d'une charge triangulaire 

Nous supposons que a < b et que P est la valeur de l'effort le plus important 
au sommet du triangle. 

i a i b , p p 

y A = -Q+b) 6 A = (l + b)(ll 2 -3b 2 ) 

6 360£/ V ; 



V B = -(l+a) B B = P (l+a)(ll 2 -3a 2 ) 
6 360£/ V ; 



Figure 4.9 



> Calcul des contraintes généralisées 

Pour 0 < x < a, 

V(x) = --^-\a(l + b)-3x 2 ] 
6a L J 

Px 

M(x) = — \a(l + b)- 3x 2 } 
6a L J 

Pour a < x < l, 

V(x) = -~ (3x 2 - 6lx + 2l 2 + a 2 ) 
6b v 7 

M{x) = ^-(l - x)(2lx - x 2 - a 2 ) 
6b v 7 



4.2 Poutre sur deux appuis 
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• , b(l + o) 
Le moment maximum est obtenu pour x =1 — \ — et vaut : 



M Max = -7zd+a)\ 



> Calcul de la déformée 

Pour 0 < x < a, 

v(x) = -^H^ [3^ 4 + ail + Z7)(7/ 2 - 3Z> 2 - 10x 2 )] 

Pour a < x <U 

v(x) = ~ ll~ E f b [3(/ - *) 4 + *(/ + ^)(7/ 2 - 3a 2 - 10(Z - x) 2 ] 

Si la charge triangulaire est symétrique, le moment maximum est obtenu en 
PI 

x = 1/2 et vaut M Ma .x = ~r • 
6 

PI 3 

La flèche maximum est égale à vil 12) = — — — , et les rotations aux 

60EI 

5 PI 2 

appuis, o A = -0 B = - — . 



4.2.7 Cas d'une charge triangulaire monotone 

P est la valeur de l'effort le plus important au sommet du triangle. 

_^rnP PI _ 1P1 3 



PI 8P/ 3 

v B = — e B = 

3 360EI 

Figure 4.10 
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> Calcul des contraintes généralisées 

P 



V(x) = -^(l 2 -3x 2 ) 
M(x)=^(l 2 -x 2 ) 

l 2P1 2 

Le moment maximum est obtenu pour x = — ■= et vaut M Max = — . Au 

v3 18y3 

milieu de la poutre, en x = 1/2, le moment est égal à M (1/2) = PI 2 / 16. 
> Calcul de la déformée 



»«=-3^(' 2 -' 2 )(" 2 - 3 ' 2 ) 
Le déplacement maximal est obtenu en x = 0.519/ et vaut : 

1PP 



VMa 



16.6EI 



4.2.8 Cas d'une charge triangulaire anti symétrique 

P est la valeur de l'effort le plus important au sommet du triangle. 

P PI PI 3 



B 6 360£/ 

V fi = 0 B = 

6 360£/ 

Figure 4.1 1 



4.2 Poutre sur deux appuis 



> Calcul des contraintes généralisées 

V(x) = -^z(6x 2 -6lx + l 2 ) 
61 v y 

M(x) = ^-(l - x)(l - 2x) 
61 

Le moment maximum est obtenu pour x = — - — / et vaut : 

6 

M Max -±-^. 

> Calcul de la déformée 



Px 

v(x) = ^^jj (6x 4 - 15lx 3 + I0l 2 x 2 - l 4 ) 



4.2.9 Cas d'une charge trapézoïdale symétrique 

P est la valeur de l'effort le plus important du trapèze. 

/iiiiiiiiiiv 2 24£/ 

AJ^111111^5 p{l-a) P . a o 

Figure 4.12 

> Calcul des contraintes généralisées 

Pour 0 < x < a, 

= - — (al -a 2 -x 2 ) 
2a v J 

M(x) = ^(3al - 3a 2 -x 2 ) 
6a v 7 
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Pour a < x < l — a, 



V(x) = ~Q-2x) 
M(x) = ^(3lx - 3x 2 -a 2 ) 



V(x) = -—[(l-x) 2 -a(l-a)] 



Pour l — a < x < l, 

P 

M(x) = ^-(l - x) \3a(l -a)-(l- x) 2 ] 
6a 1 J 

Le moment maximum est obtenu en milieu de poutre pour x = 1/2 et 
p 

vmtM Max = —(3l 2 -4a 2 ). 
> Calcul de la déformée 

La flèche maximale est obtenue en milieu de poutre et vaut : 
v(l 12) = — (5/ 2 -4a 2 ) 2 

4.2.10 Cas d'une charge parabolique 

P est la valeur de l'effort au sommet de la parabole. 

P pi p/3 



e A 



pi pp 

V B 



3 30EI 
Figure 4.13 



4.2 Poutre sur deux appuis 
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> Calcul des contraintes généralisées 

P 

P 



M(x) = ^2^(^ 3 -2/x 2 + / 3 ) 

Le moment maximum est obtenu en milieu de poutre pour x = 1/2 et 
5 PI 2 

vaut M Max = 

> Calcul de la déformée 

Px 

v(x) = (x 5 - 3lx 4 + 5l 3 x 2 - 3/ 5 ) 

La flèche maximale est obtenue en milieu de poutre et vaut : 

6AP1 4 



VQ/2) : 



516EI 

4.2.11 Cas d'un couple en un point quelconque 



/~t Va = j dA = -érS 12 -^ 



Figure 4.14 



j2 > Calcul des contraintes généralisées 

Oh 

3 Pour 0 < x < a, 

i C 
-g V(x) = -- 

§ 1 
Q Cx 
© M(x) = — 
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Pour a < x < l, 

C 

V( X ) = -j 



> Calcul de la déformée 

Pour 0 < x < a, 



Pour a < x < /, 



v(x) = -^jj [x 3 - 3lx 2 + (2/ 2 + 3a 2 )x - 3a 2 l] 

Le déplacement au point d'application du couple en x = a vaut : 

1 Cab{b — a) 



v(a) 



El 31 



4.2.12 Cas d'un couple à une extrémité 

r c ci 

c ( a b v A = - eA = — 

^ Z^L l 3« 

A Vb\ C Cl 

V B = -- #5 



Figure 4.1 5 



> Calcul des contraintes généralisées 



C 

V(x) = - 7 



4.2 Poutre sur deux appuis 



> Calcul de la déformée 



v{x)=-^-{l-x){2l-x) 
6EII 



Le déplacement au milieu de la poutre en x = 1/2 vaut : v(l/2) 

4.2.13 Cas d'un couple uniformément réparti 

C est le couple appliqué sur la poutre par unité de longueur. 



A , , , , F, , , , , , B 



jblCCCC C ((Uk v A = c 
eT~ e B \ 



e A = o 

v B = -c e B = o 



Figure 4.16 

> Calcul des contraintes généralisées 

V(x) = -C 
M(x) = 0 

> Calcul de la déformée 



v(x) = 0 
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4.3 POUTRE CONSOLE 

4.3.1 Cas d'une charge concentrée 




Va = P 



M A = Pa 



Figure 4.17 

> Calcul des contraintes généralisées 

Pour 0 < x < a, 

V(x) = -P 
M(x) P(x - a) 

Pour a < x < l, 

V(x) = 0 



M(x) = 0 



> Calcul de la déformée 

Pour 0 < x < a, 

v(x) - 

Pour a < x < l, 

v(x) - 

Le déplacement en B vaut v(l) 



Px z 
6ËÏ 

Pa 2 
6ËÏ 

2 

pa 
6ËI 



(3a — x) 



(3x — a) 



(31 — a) et sous le point d'applica- 



tion de la charge v(P) 



pa 
3Ë1' 



4.3 Poutre console 

4.3.2 Cas d'une charge uniformément répartie 

p est une charge linéique. 



P 




Figure 4.18 

> Calcul des contraintes généralisées 

V(x) = p(x - l) 
M(x) = -^(x-l) 2 

> Calcul de la déformée 

v(x) = — \(l - x) 4 + 4l 3 x 

Le déplacement en bout de poutre en B vaut vil) - 

4.3.3 Cas d'une charge triangulaire croissante 

P est la valeur maximum de l'effort. 




Figure 4.19 



8£7' 
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> Calcul des contraintes généralisées 



M(x) : 

> Calcul de la déformée 



P(l - x) 2 (2l + x) 
61 



v(x) 



Px 2 



(20/ 3 - 10/ 2 jc+x 3 ^ 



120EII 



Le déplacement en bout de poutre en B vaut v(l) 



IIP/ 4 
Ï2ÔË7' 



4.3.4 Cas d'une charge triangulaire décroissante 

P est la valeur maximum de l'effort. 



Va 



M a 



Figure 4.20 



> Calcul des contraintes généralisées 



pi 

6 



V(x) = --(l-x) 2 

M{x) = -?-{l-xf 
61 



4.3 Poutre console 



61 



> Calcul de la déformée 



v(x) 



120EII 



[Al 5 - 5l\l -x) + (l- x) 5 



Le déplacement en B vaut vil) = 



PI 4 
3ÔËJ' 



4.3.5 Cas d'un couple 




Figure 4.21 



o 

Oh 



> Calcul des contraintes généralisées 

pour x < a 

V(x) = 0 
M(x) = C 

pour a < x < l 

V(x) = 0 
M(x) = 0 

> Calcul de la déformée 

pour x < a 



v(x) 



Cx 2 
2ËÏ 
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pour a < x < l 



Ca 1 

Le déplacement au point d'application du couple vaut v(a) = — et 

2EI 

Ca ( a\ 

en bout de poutre en B : v(l) = — — yl — — j . 



4.4 ARC PARABOLIQUE ISOSTATIQUE 

4f 

Nous définissons un arc par sa fibre moyenne, à savoir : y = ~j~ x Q — x). 

4f 

L'angle d'inclinaison de la fibre moyenne est défini par : 9 = ~j~Q — 2x). 



l 

Figure 4.22 

4.4.1 Cas d'une charge uniformément répartie 

p est une charge linéique. 

N(x) = p — sin# 
V(x) = p ^x — Ç^j cos 6 
M(x) = PX(l ~ X) 



4.4 Arc parabolique isostatique 



P 

Va = V B = 
B 1 

Ha — Un — 0 



/ 

Figure 4.23 



Le déplacement au centre de la poutre est égal à v (1/2) = 
4.4.2 Cas d'une charge ponctuelle horizontale 

^ ^ V * = ~ v ' = -T Hb 

l J H A = 0 

Figure 4.24 

> Calcul des contraintes généralisées 

pour x < a 

Pb 

N(x) = — sin 0 
Pb 

V(x) = — cos 0 



pour a < x <l 



Pb 

N(x) = — smO 
Pb 

V(x) = — cos 0 
Pbx 

M(x) = P(y - b) 
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Le déplacement au centre de la poutre est égal à : 
v(l/2) = PQ-d) 



4f 2 / 3a 2a 2 12a 3 8a 4 \ 1 

' 2 - — + — + — 77-) + 



15EI V / l 1 P l 4 ES 



4.4.3 Cas d'une charge ponctuelle verticale 



P PQ - a) 

V A = — H A =H B =0 

B 1 



Pa 

J Vb = — 

l 1 

Figure 4.25 



> Calcul des contraintes généralisées 

pour x < a 

PQ - a) 
N(x) = — sin 0 

V(x) = - P(l - a) œs6 



pour a < x < l 



Pa 

N(x) = — sinO 

V(x) = - ^cos0 
Pa 

M(x) = —Q-x) 
Le déplacement au centre de la poutre (si a < 1/2) est égal à : 

vil/2) = -^JÏÂ ( /4 ~ 13 x ~ 8/2x2 + 201x3 ~ 12.8x 4 ). 



Chapitre 5 

Systèmes hyperstatiques 



5.1 GÉNÉRALITÉS 

5.1.1 Degré d'hyperstaticité H 

a) Treillis plans isostatiques et hyperstatiques 

Un treillis est isostatique s'il est possible de calculer toutes les actions de 
liaison (et les efforts dans toutes les barres) étant donnée une sollicitation 
extérieure quelconque. Pour cela, il faut faire le bilan des équations indé- 
pendantes et des inconnues indépendantes. 

Le treillis est décomposé en barres et en nœuds auxquels chaque barre est 
liée (figure 5.1). 



3 



3 




Figure 5.1 Décomposition du treillis en nœuds et barres 
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Puis on fait le bilan des équations : 

- on isole chaque barre : 3 équations x nb de barre (b). 

- chaque nœud impose une restriction selon le type de liaison : 



nœud 

p = 2 pour une articulation liant 2 barres 
p = 3 pour un encastrement liant 2 barres 
p = 1 pour un appui simple. 

- nombre d'inconnues : actions de liaison à l'extrémité de chaque barre, 

nœud 



action de liaison avec l'extérieur. 2 / 



J 



barre 



Dans notre cas le nombre d'équations est : 



b = 3 x 3 

Sa- = 2 + 2 + 2 = 6 




le nombre d'inconnues est : 




est liée 



> 15 inconnues 



nœud 1 , ext : 1 (appui simple) 
nœud 2, ext : 2 (articulation) 



Le treillis est donc isostatique. 

Cette formule est valable pour les treillis et aussi dans d'autres situations 
(cas des ossatures composées de poutres). Dans le cas des treillis elle se 
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simplifie car toutes les liaisons entre barres sont des articulations : 
nb de nœuds 

ÏNt = 2 x 

S// = i^Jj + 

chaque barre inconnues de liaisons 
à 2 articulations avec l'extérieur 



Pour un treillis isostatique : 

3b + 2n = 4b + m 



b = 2n — m 



b) Cas général 

Un système comprenant k solides possède dans l'espace 6k degrés de liberté 
(dans le plan, 3k). Ces solides sont cependant liés entre eux à différents 
points de jonction et fixés en un certain nombre de points d'appuis. 

Soit r le nombre de réactions d'appuis, c'est-à-dire de composantes de 
réactions inconnues aux différents points d'appui. Soit p le nombre de liai- 
sons internes indépendantes entre solides. On appelle degré d' hypers taticité 
la quantité : 

H = p + r — 6k dans l'espace 
H = p + r — 3k dans le plan 

Pour déterminer le nombre de liaisons internes (p), on peut utiliser les for- 
mules suivantes : 

- articulations ou rotule : si k solides concourent en une articulation 
p = 3(k — 1) dans l'espace et p = 2(k — 1) dans le plan. 

- encastrements : si k solides concourent en un encastrement 
p = 6(k — 1) dans l'espace et p = 3(k — 1) dans le plan. 

- liaisons mixtes : dénombrement à la main du nombre de conditions de 
liaisons. 
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5.1 .2 Méthode des forces 1 

a) Choix d'une structure isostatique associée 

Une structure est isostatique si les seules équations d'équilibre permettent 
de déterminer de façon unique les efforts intérieurs. Lorsque la structure 
présente des liaisons internes ou des appuis supplémentaires, la solution des 
équations d'équilibre est un espace affine, l'espace des efforts statiquement 
admissible SA. La dimension de cet espace est le degré d'hyperstaticité de 
la structure, c'est-à-dire le nombre de liaisons "en trop". Un élément de cet 
espace s'écrira : 



où M 0 est une solution particulière des équations d'équilibre et {M;, i = 
l,...w} est une base de l'espace vectoriel des efforts intérieurs auto- 
équilibrés (on ne considère que les efforts de flexion pour simplifier la 
présentation). 

Comme en mécanique des milieux continus, la résolution du problème 
nécessite dans ce cas l'introduction d'équations supplémentaires : le com- 
portement. En élasticité on montre, comme dans le cas continu, que la solu- 
tion du système minimise l'énergie complémentaire sur l'espace (principe 
des travaux virtuels). La méthode des forces consiste donc à : 

- paramétrer l'espace SA : on utilise pour cela la méthode des coupures, 

- écrire la minimisation de l'énergie complémentaire : on aboutit à un 
système linéaire faisant apparaître la matrice de souplesse. 

- résoudre le système et reconstruire la solution à partir du paramétrage 



Exemple. Soit une poutre (AB) de longueur 21 encastrée à gauche, simple- 
ment appuyée à droite, soumise en son centre à une charge verticale F diri- 
gée vers le bas. On utilisera la réaction verticale V B = X comme paramètre. 



n 




(5.1) 



i = \ 



de SA. 



1. avec l'aimable contribution de notre collègue Bruno Sudret (EDF) 
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Y 



À 






' A 




X 


< 


C 


-T 



21 

Figure 5.2 Poutre droite hyperstatique de degré 1 
P 



ha 



M x 



m A| 



Figure 5.3 Paramétrage de la poutre droite hyperstatique de degré 1 



L'équilibre s'écrit, avec les notations de la figure précédente : 

H A =0 
V A + V B -F = 0 
Ma + Fl — 2V B l = 0 



(5.2) 



On peut donc paramétrer le système de réactions statiquement admissibles 
par : 

H A =0 V A = F - X V B = X M A = (2X - F)l (5.3) 
Les efforts intérieurs statiquement admissibles sont alors : 

Af = 0 

V = X - F sur [0, /] V = X sur [/, 21] (5.4) 
M = F{x — l) + X(2l - x) sur [0, /] M = X(2l - x) sur [/, 21] 
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b) Méthode des coupures 

> Choix d'une structure isostatique associée 

Soit (S) une structure n fois hyperstatique soumise à un chargement exté- 
rieur (forces ponctuelles, densités linéiques, etc.). On lui associe une struc- 
ture isostatique (SX) en effectuant n coupures. Une coupure est obtenue : 

- soit en libérant une liaison d'appui. Sur la poutre précédente (h=l), on 
peut par exemple supprimer l'appui à l'extrémité de droite ou encore 
introduire une rotule sur l'extrémité de gauche. 

- soit en annulant une composante d'effort intérieur dans une section 
donnée : toujours sur le même exemple, on peut introduire une rotule 
pour annuler M sur n'importe quelle section x G]0, 2/[. 

On voit qu'il n'y a pas unicité de la structure isostatique associée. 

> Caractérisation de SA 

Pour obtenir une solution particulière M 0 des équations d'équilibre sur (S), 
il suffit de déterminer les efforts intérieurs sur (SX) dus au chargement exté- 
rieur . 

Pour obtenir une base de l'espace des efforts auto-équilibrés, il suffit de 
déterminer les efforts intérieurs M t sur (SX) dus à l'application de forces 
unitaires 1/ auto-équilibrées aux lèvres de la coupure. Les différents cas de 
charge à considérer en fonction de la coupure choisie sur la poutre h = 1 
sont donnés figure 5.4. 

c) Méthode des forces 

> Détermination des inconnues hyperstatiques - Théorème 
de Menabrea 

L'énergie complémentaire associée à un élément de SA s'écrit : 




(5.5) 
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sx 2 



ST 3 



4^ 



Figure 5.4 Cas de charge à considérer en fonction du choix de (SX) 



Soient (jci, . . .x n ) les valeurs des paramètres qui caractérisent la solution 
élastique dans l'espace SA. Par définition, ce sont les valeurs des efforts 
intérieurs ou des réactions d'appui relâchés par l'opération des coupures. 
On les appelle inconnues hyperstatiques. 

Le principe de minimisation de l'énergie complémentaire exprime que 
les valeurs (Xi, . . . X n ) minimisent (5.5). Les conditions de minimisation 
s'écrivent : 

dS* 

V/ = 1, . . . n, 



dxi 



0 



Xi=Xi 



(5.6) 



ce qui donne : 

V/ = 1, . . . n, 



J m -Êr dx+ ]: x >J m 



*£dx=0 (5.7) 



Ces équations constituent le théorème de Menabrea. Ce système peut se 
mettre sous la forme synthétique suivante : 



[5].X + 5o = 0 



(5.8) 
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où les termes de la matrice de souplesse [8] et du vecteur 8 0 sont donnés 
par : 



La matrice de souplesse est symétrique, définie positive et donc inversible. 
Ce qui signifie que le système (5.8) admet toujours une solution unique. 

Les termes de l'équation (5.9) peuvent être interprétés comme les ouver- 
tures des coupures pour les différents systèmes de chargement. Le théorème 
de Menabrea montre que la somme de ces ouvertures est nulle : les incon- 
nues hyperstatiques sont donc telles qu'appliquées à la structure (SX), elles 
referment les coupures (et redonnent donc la structure (S) initiale). 

> Calcul d'un déplacement sur une structure hyperstatique 

Le déplacement d'un point P de la structure peut se mettre sous la forme : 



où M et M désignent respectivement les moments fléchissants calculés sur 
(S), associés au chargement extérieur T et au chargement unitaire \ P . La 
méthode des coupures permet de mettre M sous la forme M = Mq + 



Xi M t . On démontre ainsi le théorème de réduction (ou de Pasternak) : 



pour calculer un déplacement Q P au point P, associé à un chargement exté- 
rieur T sur une structure hyperstatique , il suffit d'appliquer le théorème de 
Miiller-Breslau en appliquant l'effort unitaire lp à une structure isostatique 
associée quelconque. 




(5.9) 




(5.10) 



n 





avec Mo calculé pour une structure (SX) quelconque. 



(5.11) 



Remarque : La démarche a été illustrée en ne tenant compte que 
de l'énergie de flexion. Un calcul complet nécessite de prendre en 
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compte également l'énergie de traction-compression (et éventuelle- 
ment celle associée au cisaillement). On obtient de façon générale 
des termes de souplesse du type : 

r (^ + Mi + YiXL) dx (5 . 12) 

lJ J (S) \ El ES GS R J 

d) Guide pratique d'application 

- Calcul du degré d' hypers taticité de la structure n - Choix de n cou- 
pures. 

- Calcul de la matrice de souplesse [8] et du vecteur 8 0 . 

- Résolution du système linéaire. 

- Reconstruction des diagrammes d'efforts intérieurs. On notera qu'on 
peut très bien n'avoir pris en compte que l'énergie de flexion dans 
le calcul des inconnues hyperstatiques, alors qu'on déterminera par 
reconstruction tous les diagrammes, y compris ceux d'effort normal et 
d'effort tranchant. 

Exemple d'application. 

On considère la poutre de l'exemple précédent 5.2, soumise à une force 
verticale F dirigée vers le bas appliquée en son milieu. On utilise la structure 
(SX) de la figure 5.5. 




Figure 5.5 Application de la méthode des forces 
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À l'aide des intégrales de Mohr on obtient : 

(2/) 3 



Su = 
<5io = 



3EI 
(5F/) 3 
6EI 



(5.13) 



Ce qui donne pour l'inconnue hyperstatique X\ = 

diagramme des moments reconstruit, on utilise M 
donne le diagramme suivant : 



16 



:F. Pour obtenir le 



M 0 + X\M\ ce qui 




-3F//8 



Figure 5.6 Diagramme des moments reconstruit 

Pour calculer le déplacement vertical du point C, on applique une force 
unitaire dirigée vers le haut sur la structure isostatique associée, et on fait le 
produit des diagrammes (voir figure 5.7). 



SX 



■t 



M(x) 



Figure 5.7 Structure isostatique pour le calcul du déplacement 
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1F1 3 



On obtient Sr = — — • 

96EI 

5.1.3 Méthode des déplacements 

La méthode des déplacements est une méthode de résolution systématique 
des systèmes (iso)hyperstatiques. Elle débouche sur les techniques de réso- 
lution matricielle et la résolution numérique des problèmes (premier pas 
vers la méthode des éléments finis). Contrairement à la méthode des forces 
où les inconnues du problème sont les actions de liaisons hypers tatiques, 
les inconnues du problème seront les déplacements en certains points de la 
structure appelés nœuds. 



Exemple. 




>■ 



7 



4 



P 



3 



5 



8 




9 



1 

/777T7 



^7 



Figure 5.8 Exemple de structure hyperstatique 
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On souhaite résoudre le problème en calculant les déplacements aux nœuds 
1 à 9 sous l'effet du chargement { F t }. 

Pour cela, les inconnues du problème seront les déplacements aux nœuds 
1 à 9, certains de ces déplacements seront connus du fait des liaisons. Ceci 
sera à prendre en compte lors de la résolution. 

Les équations du problème sont les équations de la statique en chaque nœud 
de la struture. 

On utilisera les hypothèses suivantes : 

- structure plane chargée dans son plan, 

- petites déformations, 

- comportement élastique linéaire du matériau, 

- liaisons parfaites, poutres droites, section constante. 



a) Étude d'une poutre 

Prenons par exemple la poutre [5-8] de la figure 5.8 isolée du reste de la 
structure. 

Y 

V 58^ 
M 58^~ j 

Figure 5.9 Poutre [5-8] 



'85 
A- 



M, 



85 



8 N, 



85 



Bilan des actions : 

- le chargement extérieur p 

- les actions du reste de la structure sur cette poutre appliquées aux 
nœuds 5 et 8 : {Af 58 , V 58 , M 58 , Af 85 , V 85 , M 85 } 
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Attention : On travaille ici dans le repère local de la poutre étudiée qui 
peut être différent du repère global (poutre [1-2]). 

Objectifs : Exprimer les relations entre les actions appliquées aux nœuds 
de la poutre et les déplacements et rotations aux mêmes points. 

Déformée : Voir la figure 5.10 
Y 




Figure 5.10 Déformée de la poutre [5-8] 

Pour la barre [5-8], on peut écrire : 

( EA 

N 5 s —j-(u 5S - u S5 ) 

EA 

ddl supplémentaires 

Dans le cas d'une poutre, nous avons des ddl supplémentaires qui sont 
les déplacements perpendiculaires à la fibre moyenne et les rotations des 
nœuds. Pour le calcul des relations moment/rotation on peut procéder de la 
façon suivante : 

Posons : vs$ = = 0 (Figure 5.11) 

On va calculer les moments résultants des rotations imposées : 



r v(x) = v 5S 

{ M(x) = -M 58 + y 58 .x 
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Y 




Figure 5.1 1 Déformée de la poutre [5-8] avec u 58 = u 85 = 0 



avec le PFS : 



V 5S + V 85 = 0 

M 58 + M 85 -lV 5 s = 0 



Utilisons le théorème de Castigliano avec pour énergie élastique : 



M 2 , 

dx — - 

El 2 



i r l (-m 58 + y 58 x) 2 



On obtient alors les rotations aux appuis n 58 et H 8 5 en dérivant cette expres- 
sion par rapport à et M 8 5 respectivement : 



^58 
^85 



de 
dM 58 

de 
dM 85 



""/ïïO-T^^jCiZTO-ï)* 



ce qui donne après inversion (et en considérant £7 constant) : 

M 58 
M 85 



^58 + — ; — ^85 



/ 

4EI 



l 



l 



n 2EI n 

+ — ; — **58 



/ 



5.1 Généralités 
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En utilisant les équations d'équilibre on peut aussi calculer Vsg et Vg5 • 
6EI 



V 58 = 



/ 2 



-(n 58 + ii 8 5) = - Vt 
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Posons maintenant : vs% et ugs 7^ 0 (figure 5.12). 
Y 




Figure 5.1 2 Déformée de la poutre [5-8] avec u 58 / 0 et u 8 s / 0 



Cela revient à traiter le cas précédent en faisant une rotation de repère 
d'angle (u 85 - v 5S )/l. 

Bilan : on a obtenu les relations entre les efforts aux extrémités d'une 
poutre et les déplacements aux nœuds : 



EA 

Nu = —^(uij - u n ) 



Vu 



Ni 



6EI 



6EI 



Mu = 



4EI 



2EI , 

y 



12£7 

-71— <•> " "y) 



EA 

T 

6EI 



—i u ji ~ u n) 



Mu = 



_ 6EI \2EI 

12 % - — % + —p-( v ji - v n) 

4EI 2EI 6EI 

— n ji + - -^"("J'" - V U) 
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> Effet d'un chargement extérieur 

Les actions que l'on vient de calculer (Ny, Vy, M t j) sont les actions per- 
mettant d'obtenir les déplacements recherchés aux nœuds sans chargement 
extérieur. L'effet d'un chargement extérieur est donc inexistant aux nœuds 
(pas de déplacements supplémentaires). On considérera donc la poutre bi- 
encastrée soumise à un chargement extérieur (figure 5.13) : 



Figure 5.13 Effet d'un chargement extérieur sur la poutre [5-8] 



"58/ 



À ^85 



N 58 



"85 



Figure 5.14 Bilan des actions de liaisons pour un chargement extérieur sur la 

poutre [5-8] 

Il faut ensuite calculer les actions de liaisons dues à ce chargement exté- 
rieur (Ntj, Vfj, M® sur la figure 5.14) et les ajouter aux actions (Ny, Vy, 
Mij). 

Donnons le cas d'un effort ponctuel centré sur la poutre et le cas d'un 
effort réparti (figure 5.15) : 



P/2 

pl/ 5-A 



^-Pl/8 
Oj x 



, Pl/2 
Pl 2 /12^_ A 



Pl/2 
A ^-Pl 2 /12 



Figure 5.1 5 Deux cas de chargement sur une poutre bi-encastrée 



5.1 Généralités 
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La superposition de l'effet d'un chargement extérieur et des chargements 
aux extrémités de la poutre correspondant aux déplacements inconnus per- 
met d'écrire : 







EA/l 


0 


0 


-EA/l 


0 


0 


Uij 




V 


Vij 




0 


12EI/1 3 


6EI/1 2 


0 


-12EI/1 3 


6EI/Ï 2 


v ij 




V°- 


Mij 




0 


6EI/1 2 


4EI/1 


0 


-6EI/1 2 


2EI/1 


< 


> + < 


K 


Nji 




-EA/l 


0 


0 


EA/l 


0 


0 


Uji 






Vji 




0 


-12EI/1 3 


-6EI/1 2 


0 


12EI/1 3 


-6EI/Ï 2 


V ji 




v°. 






0 


6EI/1 2 


2EI/1 


0 


-6EI/1 2 


4EI/1 






M 0 - 

^ J l 



Ce qui donne sous forme condensée : 



^ =K Jt^ + ^ 

Action ext. aux nœuds Déplac^ aux nœuds Effet d'un chargé ext. 

sur la poutre 



La matrice de rigidité K est symétrique. Elle est indépendante du char- 
gement. 

Remarque : Les diagrammes des contraintes généralisées dans la 
poutre s'obtient par superposition. En l'absence de chargement 
extérieur : 

M(x) = -M^ + Vijx 
V(x) = -Vij 

Il suffit d'ajouter (par superposition) les diagrammes dus au charge- 
ment extérieur appliqué sur la poutre bi-encastrée. 



82 



5 • Systèmes hyperstatiques 



b) Equilibre des nœuds 

Considérons maintenant chaque nœud de la structure. Ceux-ci sont sollicités 
par les barres (action de la barre sur le nœud) + le chargement extérieur 
appliqué aux nœuds. 
Nœud 4 : figure 5.16 



-M 46 




Figure 5.16 Équilibre du nœud 4 

On obtient trois équations dans le repère global : 
Tt Tr Tt Tt "Tt : 2 équations dans le repère 

^ - ^43 - v 43 - n 46 - y 46 + f, = o glog J 

k — M 43 — M 4 6 =0 : 1 équation dans le repère glogal 

Les liaisons sont parfaites (nœuds rigides). Dans le repère global les 
déplacements (u, v, H) sont identiques pour toutes les barres attachées au 
nœud considéré (dans le cas d'un encastrement). On peut alors utiliser les 
relations entre les déplacements globaux et les actions extérieures obtenues 
pour la barre. On substitue les efforts dans les relations d'équilibre aux 
nœuds, pour obtenir un système d'équations. 

Pour n nœuds on obtient le système : 

(F} 

vecteur 3n du chargement ext. 
(sur les barres et les nœuds) 

Dans le cas le plus général, on obtient donc un système de 3n équations 
à 3n inconnues. 



^ {"w } 

matrice 3nx3n vecteur 3 n 



5.1 Généralités 



Remarques : 

- Le calcul de la matrice K s'effectue à partir de chaque matrice 
élémentaire. Celle-ci ne dépend donc pas du chargement appliqué 
et peut être utilisée pour d'autre cas de chargement ; seul { F } va 
changer. 

- {~u} et {F} contiennent des quantités connues et inconnues. Par 
exemple, au nœud 1 dans notre exemple les déplacements sont 
connus (ils sont nuls ici) tandis que les efforts extérieurs sont 
inconnus. A l'inverse, au nœud 2 les déplacements sont inconnus 
tandis que les efforts extérieurs sont connus. Ces conditions sont 
appelées « conditions aux limites » : statiques (efforts imposés) 
ou cinématiques (déplacements imposés). 

- On peut rajouter à volonté des nœuds dans la structure car pour 
chaque nœud on écrira alors 3 équations supplémentaires. La 
décomposition en nœuds et barres est donc « arbitraire ». Ceci 
s'appelle la « discrétisation » de la structure. 

- Si la discrétisation est arbitraire, peut-on choisir un nombre de 
nœuds très faible ? (voir figure 5.17) 




o 



î 

nrm 



Est-ce possible ? 



3 

rrrm 



Figure 5.17 Discrétisation minimale? 
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Oui mais il faut connaître la matrice de rigidité élémentaire pour la 
poutre de la figure 5.18 : 



Figure 5.18 Poutre dont il faut connaître la matrice de rigidité 

Dans la pratique, on ne calcule pas K pour toutes les géométries pos- 
sibles. On utilise les matrices de rigidité d'une poutre droite ce qui 
impose une discrétisation minimale avec 5 nœuds (voir figure 5.19). 



O 

2 



1 

////// 



3 

////// 



Figure 5.19 Discrétisation minimale avec 5 nœuds 



c) Méthodologie- Application 

1. Choix d'une discrétisation - Bilan du nombre de déplacements intro- 
duits (les ddl). 

2. Pour chaque barre, on écrit la matrice de rigidité et les équations : 
F = K~u + F 0 dans le repère local de la barre. 

3. Pour chaque nœud, on écrit dans le repère global à la structure les 
équations d'équilibres. 



5.2 Poutre droite à une travée 



85 



4. Passage des relations pour chaque barre dans le repère global et on 
forme le système d'équations algébriques : K~u = F . Identification 
des quantités connues et inconnues puis résolution. 

5.2 POUTRE DROITE À UNE TRAVÉE 

Dans cette partie, nous considérerons que : 

- la poutre est chargée dans son plan ; 

- la déformation de la poutre est alors également dans le plan de la 
poutre ; 

- la section et l'inertie de la poutre sont constantes {El constant). 

5.2.1 Encastrement élastique aux extrémités 

Soit la poutre de la figure 5.20, les appuis A et B sont supposés fixes tandis 
que la travée A B est encastrée élastiquement sur ses appuis. Les raideurs 
des appuis sont représentées par une constante positive k ce qui signifie que 
les rotations sont reliées aux moments par 6 A = k A M A et 6 B = k B M B . Si 
l'appui est un encastrement parfait alors k = oo tandis que pour un appui 
simple k = 0. 




Figure 5.20 Poutre à une travée avec des encastrements élastiques à ces 

extrémités 

Soient le moment fléchissant M 0 et l'effort tranchant Vb q ue produiraient 
dans la section d'abscisse x les charges appliquées si la poutre était isosta- 
tique. Le moment fléchissant M(x) et l'effort tranchant V(x) peuvent alors 
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être obtenus à l'aide des moments sur appuis par : 

M(x) = M 0 (x) + M A + M B j 

(5.14) 

On peut calculer les moments Ma et M B à partir des rotations 6 a et 6 B , 
rotations que subiraient les appuis A et B si la poutre était isostatique, à 
partir des relations suivantes : 

(a + k A )M A + f3M B = 6 A , fi M a + (y + k B )M B = -6 B (5.15) 



avec : 



1 f l (l-xf Â 
a = — — — - — dx 
l 2 Jo El 



1 ? x\ 



Ce qui donne : 



pe B + (y + k B )6 A 

M A - 



Mn = 



(a + k A )(y + k B )-l3 2 
/30 A + (a + k A )6 B 



(5.17) 



/ R — 

(a + k A )(y + k B )-f3 2 
Si la poutre est de section uniforme, nous avons les relations suivantes : 

a = * = ài> ? = èi (5 - 18) 



2 Poutre droite à une travée 
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.2.2 Formulaire d'une poutre simplement appuyée 
d'un côté et encastrée de l'autre 
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1 Poutre droite à une travée 

2.3 Formulaire d'une poutre bi-encastrée 
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1 



92 



5 • Systèmes hyperstatiques 




5.2 Poutre droite à une travée 
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5.2 Poutre droite à une travée 
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5.3 POUTRE CONTINUE 

5.3.1 Notations et définitions 

Une poutre continue (voir Figure 5.21) est une poutre reposant sur plus de 
deux appuis simples. Pour une poutre de n travées, on numérote les appuis 
de 0 à n. La travée i de portée U est la travée comprise entre les appuis 
A/_i et Ai, de moment quadratique I t suivant l'axe de flexion concerné, de 
module d'Young E. On appellera T\ le chargement extérieur sur la travée i. 

A 0 A l A i-1 A i A i+1 A n-1 A n 

\< —>\ \< >\< L " L ^| \< ~ >\ 

Figure 5.21 Poutre continue 

Le nombre des inconnues hyperstatiques est égal au nombre des appuis 
intermédiaires (n — 1). On peut choisir comme inconnues hyperstatiques les 
moments fléchissants sur les appuis intermédiaires. 



5.3.2 Poutre isostatique associée 

Une poutre continue comportant n travées peut être décomposée en n 
poutres isostatiques sur lesquelles s'appliquent les mêmes charges que sur 
la poutre continue avec en plus les moments aux appuis. Nous obtenons 
alors pour la travée i + 1 la figure 5.22 où : 

- Mi désigne le moment sur l'appui A t 

- Mi+\ désigne le moment sur l'appui A i+ \ 

- M i+ i(x) désigne le moment fléchissant dans la travée i + 1 de la poutre 
continue 

- M(y +i (x) désigne le moment fléchissant dans la travée i + 1 isostatique 
associée et chargée uniquement par Ti+\ sans les moments sur appuis 



5.3 Poutre continue 
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- Voj+i (x) désigne l' effort tranchant dans la travée i + 1 isostatique asso- 
ciée et chargée uniquement par J^+i sans les moments sur appuis 

- £1® d la rotation à droite de l'appui i dans la travée i + 1 isostatique 
associée et chargée uniquement par J^+i sans les moments sur appuis 

- fî? +1 la rotation à gauche de l' appui / + 1 dans la travée i + 1 isostatique 
associée et chargée uniquement par T\+\ sans les moments sur appuis 



YA 




Figure 5.22 Poutre continue 



Pour que nos poutres isostatiques associées se comportent comme la 
poutre continue d'origine, il faut écrire la continuité des rotations sur les 
différents appuis. Nous obtenons ainsi la formule des trois moments. 

5.3.3 Formule des trois moments 
a) Travées quelconques 

^{àr)Màr,^àt)Mèt) - (««) 

(5.19) 



b) Travées de mêmes inerties 

M^û + 2Mi(h + Z i+ i) + M !+1 / i+1 = 6EI fa° - n° 



(5.20) 
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c) Travées identiques 

M f _i + 4M; + M M = ^ (n% - n9^ (5.21) 

d) Travées quelconques avec dénivellations d'appui 

Les dénivellations d'appui (vt pour l'appui A t par exemple) sont comptées 
positivement vers le haut. 

M i J^) + 2M i (^ + ^-] + M i J k+l 



6EIiJ '\6EI, 6EI i+1 J "'\6EI M 

= (^-^) + ^-^i (5.22) 

5.3.4 Expression des sollicitations et actions de liaison 

Les sollicitations dans la travée hyperstatique / sont déterminées par super- 
position des sollicitations dues au chargement extérieur et celles dues aux 
moments sur appuis. 

Soit, pour le moment fléchissant, on peut écrire : 

Mi(x) = M 0i (x) + M f _i (\ - + Mi (j^ (5.23) 
De même pour l'effort tranchant : 

M- 1 M- 

Vi(x) = V 0i (x) + -p - -± (5.24) 

On déduit les actions de liaisons de l'appui A,- des valeurs de l'effort 
tranchant à droite et à gauche de celui-ci : 

Y Al = Vidù- Vi + i(0) (5.25) 
ce qui donne pour des travées quelconques : 

Y Ai = Y O i+ Mi^i_Mi^l (5 .26) 



5.3 Poutre continue 
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et pour des travées de longueurs identiques : 

Mi-i- 2Mt + Aff+i 



/ 



(5.27) 



5.3.5 Formulaire des rotations usuelles 

Cas d'une poutre sur deux appuis simples de rigidité El constante. 



Géométrie et chargement 



a 0 



i,9 



Y4- 



V A i X 



Pab(l + b) 
ËËfl 



Pab(l + a) 
ËËfl 



YA 



Ai.] 



ft 



Ai X 



pi 2 

16F/ 



PI 2 
16F/ 



Ai.i 



qP_ 
24F/ 



24F/ 





< 

q 




A|-i 


lllllllllli 


Ai X> 




^ a 





qa 2 (2l-a) 2 
24£// 



qa 2 (2/ 2 - a 2 ) 
24£// 
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a) Exemple d'application : 

Soit la poutre continue à quatre travées de la figure 5.23. Supposons les 
moments d'inertie suivants : I\ = \, h = h = 2 et Z4 = 1,5. 



2q 



q 

1111 lu- 



! l l l l l l ! 



4a 



A 2 



4a 



A 4 




-l,7835qa' 



2,4577qa" 



Figure 5.23 Exemple d'application de la formule des trois moments 



La formule des trois moments nous donne : 

(5.28) 



5.3 Poutre continue 
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avec : 



M 0 
lL l,d — 



lL 2,d 



1L 3,d 
Ma 



0 

ff(3a) 3 
24^ 
2q(4a) 3 
24EI 2 
2q(4af 
24EI 2 

2AEh 

2\Eh 
0 

0 



(5.29) 



Ce qui donne après résolution du système de 3 équations à 3 inconnues : 



Mi = -1,7835^ 2 
M 2 = -2,4577 'qa 2 
M 3 = -0, 3856qa 2 



(5.30) 



5.3.6 Formulaire de la poutre continue à 2 travées égales 



Géométrie et chargement 


Moment sur appui et réaction d'appui 


p 

A 0 ^ 




Mi = — —PI 
1 32 

13 11 3 

Y An = —P YV = —P Y Aj = P 

0 32 1 16 2 32 


1/2 




P 

A 0 i 


P 


f A 2 


M-, = — — PI 
1 16 

5 11 5 

Y An = —P VV = —P Y A , = —P 
0 16 1 8 2 16 


1/2 


1/2 


4k 
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Géométrie et chargement 


Moment sur appui et réaction d'appui 


p 


p 

V V 


Ai A 2 


M 1= -lp, 

5 8 1 

Y An = -P Y A , = -P Y A , = --p 
0 6 1 6 2 6 


4^1/3 
l«-> 


1/3 


^1/3^^ 4rs 


P 

A 0 


P 

V V 


P 

Ai 


P 

V V 


A 2 


Ya 0 = \p Ya, = \p Ya 2 = \p 


4^1/3 


1/3 


^1/3^^ 1/3^ 


1/3 




P 

A 0 V 


p 


P 

f 


y 


Ai A 2 


My = -—PI 

1 64 

81 126 15 

Y An = —P Y a, = P Y a-, = P 

0 64 1 64 2 64 


1/4 


1/4 


1/4 


1/4 


P 

A 0 V 


P 


P 

f 


V 


P 

Al V 


P 


P 

f 


\ A 2 


15 

M^=--P, 
Y Ao = ^P Y A , = ™P Ya 2 = ^P 

A 0 32 A, 32 A 2 32 


1/4 


1/4 


1/4 


4u 

1/4 1/4 


1/4 


1/4 


1/4 




c 


1 




M-, = - — g/ 2 
1 48 W 

7 15 3 

Y An = —ql Y a, = —ql Y Ai = ql 

0 16 1 24 2 48 


N^N N^N 

A 0 Ai A 2 


q 












ql 2 
1 8 

Y Ao = \ql Y A ,=\ql Y Al = \ql 


N^»N N^N 

A 0 Aj A 2 
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5.3.7 Formulaire de la poutre continue à 3 travées égales 



Géométrie et chargement 



Moment sur appui et réaction d'appui 



A 0 



1/2 



Ai A2 A 3 



Mi = - — PI M 2 = —PI 

1 10 1 40 

Y An = -P VV = — P 

A ° 5 1 40 



Y A2 = 


-Ap 
20 


Y A2 = 


Ap 
40 


M,=- 


-A P , 

40 


M 2 = - 


40 


Y A 0 = 


40 


Y A, = 


40 


y a 2 = 


40 


Y A 3 = - 


Ap 

40 



A 0 



Ai 



A 2 



a 3 

lu 



A« 



1/2 



Ai 

"F 



1/2 



M, = — -PI M 2 = — —PI 
10 20 



A 2 

HT 



a 3 

lu 



Ya 0 


= lp 
5 


K», = 


Ap 

20 


Ya 2 - 


= Ap 

10 


^3 = 


-Ap 

20 


M, = 


-A P/ 

20 


M 2 = 


-Ap, 

20 


Ya 0 


= Ap 

20 


Ya, = 


Hp 

20 


Ya 2 


= 2 -l P 
20 


Ya 3 = 


Ap 

20 



A 0 



Ai 



A 2 

~3Û~ 



a 3 

"3u 



1/2 



1/2 



1/2 
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Géométrie et chargement 


Moment sur appui et réaction d'appui 


q 

N^N \^.N 

A 0 Aj A 2 A 3 


Ya 0 = f Q ql Y A , = f o C 

Y A =- — ql Y A =—ql 
2 10 W 3 60 W 


4 

\^\ \^\ 

A 0 Aj A 2 A 3 


11 1 

Y «> = To ql Y ^ = -To ql 


a a 

\^.\ N^N \^.N \^.\ 

A 0 Aj A 2 A 3 


Yao = |q/ = 

9 1 

Va, = — ql Y A, = ql 

A 2 2Q W A 3 3Q W 


q q q 

N^s\ N^N N^v\ N^N 

A 0 A x A 2 A 3 


1 10 w 1 10 w 

2 1 1 

Y Ao = - 5 q' Y A , = -ql 

^ = ^/ = §1/ 
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5.3.8 Formulaire de la poutre continue à 4 travées égales 



Géométrie et chargement 


Moment sur appui 


\^.\ 

A 0 Ai A 2 A 3 A 4 


Mt = -0, 080P/ 
M 2 = -0, 054P/ 


ri p n 

\^s\ 

A 0 Ai A 2 A 3 A 4 


Mt = -0, 090P/ 
M 2 = 0, 095P/ 
M 3 = -0, 090P/ 


llïll llîll 

\WS \^v\ 

A 0 Ai A 2 A 3 A 4 


M-, = -0,071g/ 2 
M 2 = 0, 036g/ 2 
M 3 = -0,071g/ 2 


q q 

\W\ \W\ 

A 0 Ai A 2 A 3 A 4 


Mt = -0, 054g/ 2 
M 2 = —0, 036g/ 


q q q 

N^N X^>X \NAX \NAX N^AN 

A 0 A x A 2 A 3 A 4 


Mt = -0,121g/ 2 
M 2 = -0,018g/ 2 
M 3 = -0, 058g/ 2 
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Formulaire de la poutre continue à 5 travées égales 



Géométrie et chargement 


Moment sur appui 


p I p I 

i i 

N^N N^vN \^sN \^s\ \^sN N^N 

A 0 Ai A 2 A 3 A 4 A 5 


M-, = -0,1 06P/ 
M 2 = 0, 048P/ 
M 3 = -0, 088P/ 
M 2 = -0, 072P/ 


PI PI PI PI PI PI 

Y Y Y Y . .Y Y . 


M-, = -0, 140P/ 
M 2 = -0,1 05P/ 
M 3 = -0, 105P/ 
M 4 = -0,1 40P/ 


\^v\ \WS 

A 0 Ai A 2 A 3 A 4 A 5 


q 

| | | | | | | | | | || | | | | || 

\W\ \^v\ \WS \WS 

A 0 Ai A 2 A 3 A 4 A 5 


Mt = -0,105g/ 2 
M 2 = -0, 079g/ 2 
M 3 = -0, 079g/ 2 
M 4 = -0,1 05g/ 2 


q q q 

\W\ 

A 0 Ai A 2 A 3 A 4 A 5 


M-, = -0, 053g/ 2 
M 2 = -0, 039g/ 2 
M 3 = -0, 039g/ 2 
M 4 = -0, 053g/ 2 



5.3 Poutre continue 
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5.3.10 Poutre continue sur appuis élastiques ponctuels 
a) Cas général : formule des cinq moments 

Par hypothèse, les dénivellations d'appui sont dans ce cas proportionnelles 
aux réactions d'appui : 



Vi = -kiY Ai 
On obtient donc à l'aide de la formule (5.26) : 



(5.31) 



M t Mi - M M 



U 



+i 



(5.32) 



En remplaçant les expressions de v t dans la formule des trois moments 
avec dénivellation d'appui on obtient la formule des cinq moments : 



U-iU 



6 EL 



kj-if 1 | l\ _ VI + _L 

h \h-l h ) h \h k+1 



6EIi ■ 6EI M J " If {li + l M ) + lf +l 



k+l 



h (\ 1 
+ 



k i+ i ( \ + 1 



6EI i+ iJ \li J li+2 

k i+ i 



U+iU 



-Mi_ 2 

Mi 

M M 
M i+2 



(5.33) 



+2 



On utilise cette formule notamment dans l'étude des systèmes de poutres 
croisées. 
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b) Cas d'une poutre infiniment rigide 

L'application de la formule des cinq moments à une poutre dite infiniment 
rigide mettrait en évidence des déplacements v t et en conséquence des réac- 
tions d'appui Y Ai fonctions linéaires des abscisses x. Soit : 



Pki 



1 ^ e ^ ki • 
1 + ^ „ — xi 



U = ËU 1+ Ë?rv (534) 

en supposant l'origine O choisie de sorte que ^^ Xjkj = 0 et en désignant 
par P la résultante de l'ensemble des charges appliquées et par e sa position 
(voir figure 5.24). 



i 



X 



F £ o ? î î~ 

////// ////// 17/777 777777 777777 



Figure 5.24 

Les sollicitations en une section 2 de la poutre sont obtenues en faisant 
une coupure dans cette section et en écrivant l'équilibre des forces de gauche 
(ou de droite). 



5.4 SYSTÈMES DE POUTRES CROISÉES 
5.4.1 Principe 

Les systèmes de poutres croisées sont calculés par l'une des méthodes sui- 
vantes : 

- la méthode directe, c'est-à-dire l'écriture de l'égalité des déplace- 
ments verticaux aux point de croisement des poutres, sous l'action des 
charges appliquées, 



5.4 Systèmes de poutres croisées 
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- la résolution de l'équation des cinq moments (équation 5.33), 

- l'utilisation de codes de calcul appliqués sur un modèle faisant inter- 
venir les liaisons à chaque croisement des poutres. 



5.4.2 Cas particulier des poutres de même inertie 

Si on suppose que toutes les poutres sont simplement appuyées, de même 
inertie, de raideur de torsion négligeable et soumises à un chargement P à 
chaque intersection, l'application des méthodes précédentes conduit à des 
résultats simples dans les cas suivants : 



Géométrie 



Poutres A 



Poutres B 



M - 



b 3 



4(a 3 + b 3 ) 

.. b3 / 

2(a 3 + b 3 ) 



Pis 



M - 



4(a 3 + b 3 ) 

-. * 3 / 

2(a 3 + b 3 ) 



PI h 



M - 



b 3 



3(5a 3 + b 3 ) 
5a 3 + b 3 



Ma 



M-- 



5a 3 



4(5a 3 + fa 3 ) 
2(5a 3 +b 3 ) 



PIh 



M - 



b 3 



3(a 3 +b 3 ) 

= -^-1 
a 3 + b 3 



Pis 



M - 



3(a 3 + b 3 ) 
a 3 +fa 3 



Mh 
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5.4.3 Cas particulier des poutres infiniment rigides dans une 
direction 



g g g g g 









































3 £ 



Figure 5.25 

Soit le cas particulier du système de poutres croisées représenté à la figure 
5.25. Les poutres B sont supposées infiniment rigides (admissible lorsque 
lb < 21 a et si la hauteur des poutres B ~ hauteur des poutres A). Dans ce cas 
particulier, on peut considérer qu'une poutre B peut être considérée comme 
rigide et appuyée sur des appuis élastiques ponctuels. Nous pouvons donc 
utiliser l'équation (5.34) pour déterminer les réactions 1%. Dans la mesure 
où les poutres A suivent la même loi de variation d'inertie, on peut substituer 
aux coefficients k t les inerties des poutres . Nous obtenons alors : 



Pli 

EitV ' Mit 



1 + ^ o xi 



(5.35) 



5.5 POUTRE SUR APPUI ÉLASTIQUE CONTINU 

5.5.1 Définition et paramètres 

On considère la poutre de largeur b, d'inertie constante, reposant sur un 
appui continu élastique (de type sol) de la figure 5.26. 

Y t q(x) 

■""m * 

r 

Figure 5.26 



5.5 Poutre sur appui élastique continu 
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Si y désigne la flèche de la poutre au point d'abscisse x, la réaction d'ap- 
pui en ce point est équivalente à une densité d'effort r : 

r = -kby (5.36) 

où b est la largeur de la poutre supposée constante, k est la raideur de l'appui 
élastique. Si q(x) est la charge répartie appliquée à la poutre, la flèche doit 
satisfaire l'équation différentielle d'équilibre suivante : 

EI^+kby + q(x) = 0 (5.37) 
Pour résoudre cette équation, on pose : 



1 A kb 



h V AEI 



(5.38) 



La longueur l e est appelée longueur élastique de la poutre. Les solutions 
y(x), y'(x), M(x) et V(x) font intervenir des fonctions W(x), A'(x), y(x) et 
Z(x) suivantes : 



pour x > 0 : W(x) = -e A *(cos(Àx) + sin(Ax)) 
X(x) = e~ Xx sin(Ax) 
y(x) = e~ ÀX (cos(Xx) - sin(Ax)) 
Z(x) = -e~ Kx cos(Ax) 



pourx < 0 : W(x) = -^ Ax (cos(Ax) - sin(Ax)) 



(5.39) 



X(x) = e ÀX sin(Ax) 

y(x) = e ÀX (cos(Ax) + sin(Ax)) 

Z(x) = e ÀX cos(Ax) 



(5.40) 
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5.5.2 Formulaire de la poutre infinie 



Géométrie et chargement 


Solutions 


l , 

0 


/ PA^ , » PA , » 
M=^y(x) V=^Z(x) 


Y> 


< 

C 

X 


Ta3 TA 2 

M=^Z{x) V = -yW(x) 


0 






pour 0 < x < a : 

y' = ^[W(a+x)-W(a-x)] 

y = -^77- [2 + Z(a +x) + Z(a - x)] 
2/CD 

M= ^_ [AT(a+x) + Af(a-x)] 

^=-^ [-Ma+x) + Ma-x)] 
pour x > a : 

/ = ^ [W(x + a)-W(x-a)] 

y=-jj^ [Z(x + a)-Z(x-a)] 

M=£^ [X{x + a)-X{x-aj\ 

V=~ [-y(x + a) + y(x-a)] 
4A L J 


0 

a 


a 

^ H 


H > 
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Géométrie et chargement 



Solutions 



YA 



/(O) : 

/ (a) : 

y(0): 
y(a) 



2kba 

g 

'ikba 

q 



4\kba 

g 

' 4\kba 



[1 + AaW(a) + Z(a)] 
[1 - Aa + Z(a)] 
[1 - y(a) + 2\aZ(a)] 



[2Aa - 1 + y (a)] 



pour 0 < x < a : 

M= [WM- W(a-x) 
bA 3 +2Aa*(a-x)] 

V = t^- [*M + #(a-x) 

4A 3 -AaMa-x)] 

pour x > a : 

M = [WM - W(x - a) 
bA 3 -2Aa*(x-a)] 

V = [Af (x) - *(x - a) 

4A 3 -AaMx-a)] 



5.5.3 Formulaire de la poutre semi-infinie 

Pour obtenir les déformations et les sollicitations dans une poutre semi- 
infinie, on utilise les résultats de la poutre infinie en appliquant une force 
Po et un couple Co juste à gauche de O (voir figure 5.27) de telle façon que 
sous l'ensemble du chargement M et V soient nuls en O. On peut écrire : 

P 0 = 4(ÀM 0 + V 0 ) C 0 = y(2AM 0 + V 0 ) (5.41) 

A 

avec M 0 et Vb le moment et l'effort tranchant dus à q(s) pour la poutre 
infinie juste à gauche de O. 
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YA 



q(x) 
I | l 1 1_ 



x 



<=> 



P Y t i (x) 
1 1 1 



Figure 5.27 



On obtient alors pour la poutre semi-infinie : 

PoA 2 



y'(x) = y[(x)+-^-X(x) + 



CqA 3 
kb 
CnA 2 



y(x) 



(5.42) 



y(x) = yi (x) + ^W(x) + ^-X(x) 
2kb kb 

M(x) = M,(x) + ^y(x) + ^Z(x) 

V(x) = V^x) + ^Z(x) - ^-W(x) 

avec y[,yu M\ et V\ les solutions obtenues pour la poutre infinie soumise 
au chargement q(x). 



Géométrie et chargement 



Solutions 



Y-f 



X 



Po = P[y(c) - 2Z(c)] 
C 0 = -\y{c)-Z{c)] 

A 

Cas particulier c = 0 : 
, 2PA 2 , . 2Pa 2 

y 1 = -77- Wx) y <°) = ~nr 
kb kb 

y = -T7-Z(x) y(0) 



kb 



kb 



M = X{x) M(0) = 0 

À 

v = -py(x) v(o) = -p 
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Géométrie et chargement 


Solutions 


Vi 

O 


c 

< ► 


C 

X 


P 0 = -2CA [Z(c) + W(c)] 
C 0 = -C[2Z(c) + W(c)] 

Cas particulier c = 0 : 
, 4CA 3 . . „ n . 4CA 3 

y = — n-^w y (o) = -T7- 

kb kb 
2CA 2 , , , n , 2CA 2 

y = - fcb Mx) y(0) = - kb 

M = CWM M(0) = -C 
V = 2CA^(x) \/(0) = 0 


( — 


— — > 


O 


< c ^ ^ ^ 


P 0 = y [Af(2a + c) - Af(c) 

-M2a + c) + Mc)] 

C 0 = [2*(2a + c) - 2* (c) 

-M2a + c) + Mc)] 
Cas particulier 2a = oo : 

= t [MO - X(0] 

A 

Cas particulier c = 0 : 
y'(x) = 0 y(x) = --5- 
M= 0 V = 0 




r 

2a 

< H 



a) Poutre semi-infinie articulée à l'origine 

Pour obtenir les déformations et les sollicitations dans une poutre semi- 
infinie articulée en son origine (voir figure 5.28), on utilise les résultats de 
la poutre semi-infinie. La réaction verticale sous l'articulation s'écrit : 



(5.43) 
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avec y\(0) le déplacement vertical du point O de la poutre semi-infinie non 
articulée sous l'effet du chargement extérieur q(x). 



Y A 



q(x) 



0 I l l 



X 



Figure 5.28 

On obtient alors pour la poutre semi-infinie articulée en son origine : 
y'(x) = y[(x)+—W(x) 

y{x) = yi{x) -— Z{x) (544) 

M(x) = Mi(x)+—X(x) 
À 

V{x) = V l {x) + Ry{x) 

avec y[,yi, M\ et Vi les solutions obtenues pour la poutre semi-infinie sou- 
mise au même chargement q(x). 



5.5.4 Formulaire de la poutre de longueur finie 

Pour étudier une poutre A B de longueur finie / soumise à un chargement 
extérieur T exU il suffit de superposer les efforts obtenus dans une poutre 
infinie de même inertie, résultant des deux cas suivants : 

- poutre infinie soumise à T ext . Détermination de Va, M a et V B , M B 
respectivement dans les sections A et B. 

- poutre infinie soumise en A à la force et au moment ponctuels Pa et 
Ca ainsi qu'à Pb et Cb en B. 

Pour obtenir un système équivalent à la poutre finie, il faut écrire que 
les efforts tranchants et les moments fléchissants doivent être nuls en A et 



5.5 Poutre sur appui élastique continu 
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B. On obtient un système de 4 équations à 4 inconnues P A , Ca, Pb, C b à 
résoudre : 



(5.45) 



P A + P B ZQ) - C A A + C B XW{1) - 2V A = 0 
P A Z(l) + P B - C A ÀW(/) + C fi A + 2V fi = 0 
Pa + P B yO) - 2C A A - 2C fi AZ(/) + 4AM A = 0 
PaW) + + 2C A AZ(/) + 2C fi A + 4M B = 0 

On obtient alors pour la poutre finie de longueur / : 

A:# A:# 

y(x) = + ^W) + ^X(x) + W(Z - x) - ^*(Z - x) 

M(x) = M\(x) + ^y(x) + ^Z(x) + ^y(Z - x) - ^Z(Z - x) 

V(x) = V x (x) + ^Z(x) - ^W(x) - ^Z(Z - x) - ^W(Z - x) 

(5.46) 

avec y{, y\, M\ et Vi les solutions obtenues pour la poutre infinie soumise 
au même chargement extérieur T ext avec l'abscisse x mesurée à partir de A. 



a) Cas particulier d'un chargement symétrique par rapport 
au milieu de la poutre AB 

La symétrie entraîne V B = —Va, M B = M A , P B = Pa et C B = — C A . Le 
système précédent de 4 équations à 4 inconnues devient un système de deux 
équations à deux inconnues Pa et C A : 

P A [l + Z(l)] - C A A[1 + W(Z)] - 2V A = 0 
PaU + y(Z)] - 2C A A[1 - Z(Z)] + 4AM A = 0 (5 ' ' 



118 



5 • Systèmes hyperstatiques 



y\x) = y[(x) + [X(x) - X(f - x)] + - y il - x)] 

ko ko 

y(x) = yi (x) + ^ [W(x) + W(/ - x)] + ^ [*(*) + *(/ - x)] 
2kb ko 

M(x) = M x (x) + ^ + y(l -x)] + ï± [Z(x) - ZQ - x)] 

V(x) = V,(x) + ^ [Z(x) - ZQ - x)] - ^ - W(Z - x)] 

(5.48) 

Remarque : 

- Lorsque la longueur de la poutre est telle que / < l e , il n'y a pas 
lieu de tenir compte de la déformation des appuis. La poutre peut 
alors être calculée comme si elle était infiniment rigide. 

- Lorsque la longueur de la poutre est telle que À/ > 5, on peut 
considérer les couples (Pa-, Ca) et (P Bj C B ) indépendants l'un de 
l'autre. Chacun peut alors être déterminé en supposant la poutre 
semi-infinie. 



5.6 PORTIQUE 

Pour la suite, on définit le paramètre suivant : 



5.6 Portique 
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5.6.1 Portique à un seul montant et à deux extrémités 
articulées 



Géométrie et Chargement 



Diagramme des moments 



Moments 



A 



M B 



Pab/IQ +b/l) 
2(1 +k) 



M P = (Pa + M B ) 



B llUlJUUUlUUUUl c 
3k 




M B - 



ql 2 



8(1 + k) 



5.6.2 Portique à un seul montant et à deux extrémités 
encastrées 



ii 

A 
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Géométrie et Chargement 



Diagramme des moments 



Moments 



A 




M A 

M b : 

M c : 
(2 



Mb 
2 

Pab b/l 
I 1 + k 
Pab 

Z r 

b/l)k + 2(1 - b/l) 



Mp 



2(1 +k) 

+ ' 1 ~ 7 



M c 



M(s) 



b I u U U 1 1 U I U I U 1 1 1 Ut 



A 



Ma : 
M B : 

Mç : 



Mb_ 
2 



ql 2 



-\2C\+k) 
ql 2 (3k + 2) 
' 24(1 + k) 



5.6.3 Portique à un seul montant et à une extrémité 
encastrée et l'autre articulée 



A 
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Géométrie et Chargement 


Diagramme des moments 


Moments 


B 


1 C 


u 


(s) 


M A = — 

h 

3Arb/A7 + 2(1 +b/h) 
4 + 3k 

M B = — 

h 

3/c(1 - 

y 7 

4 + 3/c 

M P = — + — M/\ 

+ (i-£)m. 


A 

\v 


< a > | < b ^ 


w 


] 


3 C 




(s) 


/^M — — 

Q/ 2 

Mfi = -12(1 + « 
C 24(1 + /c) 


-c 






^ 


B 

A 


uniunuiuiuiunc 


]V 


(s) 


M A = — — 
ql 2 

MB = -2(A + 3k) 






B 


c 

A 

v\ 


M 

.( 

2-^ 


h (s) 






4(4 + 3fc) 
Ms= 4(4 + 3k) 
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5.6.4 Portique à deux montants articulés 



A 

4k 



D 

4k 



Géométrie et Chargement 



Diagramme des moments 



Moments 



A 



D 




M B = M c 

Pab 



I 2(3 + 2k) 



M P = + M B 



B IUIUUIUUUUUIU C 





M B = M c 



Mrr 



ql 2 



4(3 + 2k) 



ql 2 



+ M B 



M B - 



Mr 



Ph 

Y 

Ph 

Y 



5.6 Portique 



123 



Géométrie et Chargement 



Diagramme des moments 



Moments 



A 
4u 



D 

4v 



+ 1 



Pa 
2 

k(2 - b/h)b/h 
3 + 2k 

Pa 

y 

k(2-b/h)b/h 
3 + 2k 

M P =U- ^\ (Pb + M B ) 



M B = — 



Mr 



A 
4k 



M B - 



Mr 



qh 2 



4 V 6 + 4A: 



+ 1 



qh 2 



4 V 6 + 4A: 



5.6.5 Portique à deux montants encastrés 



A 



D 
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Géométrie et Chargement 



Diagramme des moments 



Moments 



4L 



D 



m! 



M A 



M B - 



Pab 



1 2b /I - 1 



2(2 + k) 2(1 + 6k) 
Pab 

: r 

1 2b/ 1 - 1 

x | + 1 



ï 



2 + k 2(1 + 6A:) 
Pab 
T 

1 2b /I - 1 



2 + k 2(1 + 6A:) 
Pab 

~T 

1 2b /I - 1 
x I — + 



Mp : 



2(2 + k) 2(1 + 6A:) 



P B 



A 



D 



M A 



M b 



a 



Ph 1 + 3k 

~~2~1 +6A: 
Ph 3k 

YTVëk 
M c = -M B 

Mq = -M A 



A 



D 



«1 




a 2 = 
a 3 -- 
M A 



Pab 1 +b/h + kb/h 
~h 2(2 + k) 
Pab (1 -b/h)k 
~h~ 2(2 + k) 
3Pa(1 - b/h)k 



2(1 + 6/:) 
Pa 



+ «3 



Mg = —Qij + «3 
Mç = —dj — «3 

M D = + y - 



«3 
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Géométrie et Chargement 



Diagramme des moments 



Moments 



BUUUUUUUUUUUC 



A 



D 




M A =M D : 

M B = Mç - 



ql 2 



/ \2(2 + k) 

ql 2 
6(2 + k) 



D 



1 



M A 



M b 



M C : 



qh 2 



3 + k -\+4k 



6(2 + k) ï+6k 

qh 2 



k 2k 
+ 



6(2 + k) *[+6k 
qh 2 

k 2k 



6(2 + k) 1 + 6A: 

Mq = -M A 



5.7 ARCS HYPERSTATIQUES 

5.7.1 Arc circulaire à deux articulations sans tirant 

Soit l'arc circulaire hyperstatique de la figure 5.29. La géométrie de l'arc 
est un cercle de centre O et de rayon R. Nous avons également les données 
suivantes : 

- E : module d'élasticité du matériau de l'arc, 

- / : moment d'inertie d'une section d'arc (considéré comme constant), 

- S : aire d'une section d'arc (considérée comme constante), 

- G : module de cisaillement du matériau de l'arc. 
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Y 




O 

Figure 5.29 Arc circulaire bi-articulé 



Géométrie et Chargement 


Efforts et moment de flexion 


Y 


.llllllllllllllllllllllll 


Y= 1 

1 + \ x 4a 

SR l 2a — 3 sin 2a+Aa COS 2 a 

Xa = X B = %^ 

(9 - 4a 2 ) sin 2a - 1 0a COS 2a - 8a 


/ \ X 


x ± 
2a — 3 sin 2a + 4a COS 2 a 

N(cp) = qRcp sin cp + X A cos cp 
M((p) = qR 2 (a sin a - cp sin cp + cos a - cos cp) 

- X A R(C0S (p - COS a) 

V((p) = -qR(p cos ç> + X A sin ç 


aJ 




Y 


fi Aj 


X A = -X B = -qRa 

N(cp) = q/?(sin cp + cp cos cp - a cot a sin <p) 

M((p) = qR 2 (a cot « sin cp - cp cos cp) 
V(ç) = qR(a cot a cos (p + cp sin <p — cos <p) 


A4 
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5.7.2 Arc parabolique à deux articulations sans tirant 

Soit l'arc parabolique hyperstatique de la figure 5.30. Nous avons les don- 
nées suivantes : 

- E : module d'élasticité du matériau de l'arc 

- / : moment d'inertie d'une section d'arc (considéré comme constant) 

- S : aire d'une section d'arc (considérée comme constante) 

- G : module de cisaillement du matériau de l'arc 

Les équations de la parabole sont les suivantes : 



y = —z x(L -x) 




(5.49) 



L 



Figure 5.30 Arc parabolique bi-articulé 
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Géométrie et Chargement 



Efforts et moment de flexion 



X A = X B 



qL 2 



UUUIUUUUIUUUIU 



8/7(1 + 15//8S/7 2 ) 




Y A = Y B = y cos 0 - X A sin 6 
N(cp) = — x^j sin <p + X A cos<p 



1111111111111 



qL 2 



16/7(1 + 15//8S/7 2 ) 




Y A = q(~^- -x^J cos # -X A s\nd 
QL 

Y B = -— cosfl -X^sinfl 
8 

Nac(<p) = q(^- -x^J s\n(p+X A œs(p 
qL 

Nçb(<p) = - — - sin <p + X A cos <p 
8 

M cb (<p)= ^(L-x)-X A y 

O 



Chapitre 6 

Plaques et coques 



6.1 PLAQUES 

Une plaque est un élément prismatique d'épaisseur h petite devant les deux 
autres directions de l'espace. Le plan moyen sera le plan (0,x,y), le déplace- 
ment trans verse étant la direction E On supposera dans la suite l'hypothèse 
des petits déplacements vérifiée. 




Figure 6.1 
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6.1 .1 Formules générales 

On appelle v(x, y) la flèche de la plaque (suivant l'axe z). Les solutions sont 
données pour des plaques élastiques de Kirchhoff-Love, valables pour les 
plaques minces. La flèche vérifie l'équation : 

_ / d 4 v „ d 4 v d 4 v\ , ^ ^ 

avec q(x, y) la densité de charge et D la raideur en flexion de la plaque : 

Fh 3 

D 



où E et v sont les paramètres élastiques du matériau, respectivement le 
module d' Young et le coefficient de Poisson. 

> Contraintes généralisées 

Moments de flexion : 

r /d 2 v d 2 v\ _ f d 2 v d 2 v 



dx 2 dy 2 ) y \dy 2 dx 2 

d 2 v 



dxdy 



Efforts tranchants : 



M xy = M yx = D(l - v)- 

dM x | dM xy _ dM y | ÔM xy 

dx dy y dy dx 



> Contraintes 

Contraintes normales : 

12M X , , 12M 



o- x (z) = ^^-z o- y (z) = ^p^z 
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Contraintes de cisaillement : 



l2M xy 

T^Z 




6.1.2 Méthode de résolution pour les plaques rectangulaires 

On considère une plaque rectangulaire de dimension a x b, avec a < b. La 
méthode de résolution se base sur la solution analytique associée au charge- 
ment réparti de la forme : 



y) = qo sin yn — J sm J 



La solution pour tout type de charge peut alors être obtenue en décomposant 
celle-ci sous la somme : 



q(x,y) = ^2 ^2 a mn sin sin ('?- 1 



m = 1 n- 



7ry_\ 

b J 



La flèche v(x, y) solution de l'équation (6.1) s'écrit alors : 



v(x, y) = ^2^2 A mn sin ( m ^-) sin (n— 1 ] 



b J 



Les coefficients A mn sont déterminés en fonction des conditions aux limites. 
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6.1.3 Plaques rectangulaires 

a) Simplement appuyées sur les quatre côtés 
> Charge ponctuelle au centre de la plaque 



y* - 
U 



Figure 6.2 



Flèche maxi (point 0) 


Moment maxi (point O) 


Fb 2 


4 



a/b 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


oo 


a 


0,1265 


0,1478 


0,1621 


0,1714 


0,1769 


0,1803 


0,1849 




0,149 2 


0,2101 


0,261 6 


0,3062 


0,345 6 


0,380 8 


0,420 5 



> Charge uniformément répartie 



y* - 



Figure 6.3 



Flèche maxi (point O) 


Moment maxi (point O) 


qb 4 


f3qb 2 



a/b 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


oo 


a 


0,044 3 


0,061 6 


0,077 0 


0,0906 


0,101 7 


0,1106 


0,1422 




0,047 9 


0,062 6 


0,075 3 


0,0862 


0,094 8 


0,101 7 


0,1250 
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b) Encastrées sur les quatre côtés 

> Charge ponctuelle au centre de la plaque 



////s///////////// t 
B 




Figure 6.4 



Flèche maxi (point O) 


Moment maxi (point B) 


Fb 2 


4 



a/b 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


oo 


a 


0,061 1 


0,070 6 


0,075 5 


0,077 7 


0,078 6 


0,078 8 


0,079 2 




-0,1257 


-0,1490 


-0,1604 


-0,165 1 


-0,1667 


-0,1674 


-0,1680 



> Charge uniformément répartie 




Figure 6.5 



Flèche maxi (point O) 


Moment maxi (point B) 


qb 4 


flqb 2 



a/b 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


oo 


a 


0,0138 


0,0188 


0,022 6 


0,025 1 


0,026 7 


0,027 7 


0,0284 




-0,051 3 


-0,063 9 


-0,072 6 


-0,078 0 


-0,081 2 


-0,082 9 


-0,083 3 
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6.1.4 Plaques circulaires 

a) Bord simplement appuyé 

> Charge ponctuelle P au centre de la plaque 



2a 

< ► 



2a 



Figure 6.6 



Contraintes : 



Flèches : 



3P a 
o- rr (r) = - — y(l + v) In - 



lire 2 



3P / a 
<roo(r) = - — r (1 - v) + (l + v)\n- 
Lire L V r 

v(r) = -J— (—(a 2 - r 2 ) + 2r 2 ln - 
16e D \1 + v v y a 

_ 3Pa 2 (3 + ^)(l - i/) 

Umax " ÏËTôTT) 

> Charge uniformément répartie g sur la plaque 



2a 

< > 



2a 



Figure 6.7 
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Contraintes 



Flèches : 



(Too(r)=^((3 + v)a 2 -(l + 3v)r 2 ) 

3q 2 

O-max = O-rr(O) = <r 00 (O) = ^(3 + v)d 



-fl 0 



64 D \\ + v 

> Charge sur la circonférence d'un cercle concentrique de rayon c 



2</ 



Figure 6.8 



_£ — * 



< 2a > 



Contraintes : 

sir < c, (T rr (r) = a 00 (r) 

si c ^ r ^ a, cr rr (r) 
3P 



3P 
Aire 1 

3P 



(1 - *0 1 



2(1 + ln- 



477^ 2 *\r 2 a 2 / 2(l + ï>)ln^ 



477é> 2 



d - *0 2 



2 2 



2(1 + ^) ln 
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Flèches : 



si r < c, u(r) — - (a 1 - c 1 ) — — + (r 1 + c 1 ) ln - 

877Z) V 2(1 + ^)a 2 v 7 a 



v(r) 



8ttD 



(a 2 - r 2 ) 1 + 



si c ^ r ^ a , 

1 — j> a 2 — c 2 



8ttD 



2(1 + i/) a 2 

( ( 2 2\ (3 + *0 2l C A 

(a z - c z ) — + c 2 ln - 

V v J 2{\ + v) a) 



+ (r 2 + c 2 ) ln 



> Charge uniformément répartie à l'intérieur d'un cercle 
de rayon c 



'U 



Contraintes : 



2a 

< > 



Figure 6.9 



q 



2a 



si r < c, cr rr (r) = ^ ( 4 - ( 1 - v)- - (3 + f )- - 4(1 + f) In ■ 



<ree(r)=^- ( 4 - (1 - v)^ - (1 + 3r£ - 4(1 + p)\n ° 



8e 2 



si c < : r < : a, a n -(r) = f (1 - v) ( ^ - 1 j ^ - 4(1 + v)ln ^ 



(r "" lr) 's^ I " J,) ( 4 



4 + l ) 4) -4(l + f)ln- 
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Flèches : 



si r < c, 



, N 4c 2 / r 4 4a 2 r 2 - (1 - ^)c 2 r 2 4(3 + v)a 2 - (7 + 3^)c 2 

TTTT I — — — + 



16D V4c 2 



2(1 + j/)a 2 



4(1 + v) 
+ (c 2 + 2r 2 ) ln 



v(r) 



qc 
Ï6D 



si c ^ r $C a, 
3 + ^o 1 — z> 



a 2 y V 1 + ^ 2(1 + 1/) 



c z ) + (c 2 + 2r 2 )ln 



b) Bord encastré 

> Charge ponctuelle P au centre de la plaque 




Figure 6.10 



4 



2 a 



o 



Contraintes 



Flèches : 



arr(r) = ^ ((l + f)In^-l) 
3P 



o-eeir) 
v(r) = 



lire 2 
P 



!6eD 



[a 



2 r 2 + 2r 2 ln 



-) 

a/ 



3Pa 2 (l - y) 
4Ee 4 
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> Charge uniformément répartie q sur la plaque 




mm 



Figure 6.1 1 



Contraintes 



Flèches : 



°rr{r)=^({\ + v)a 2 -0 + v)r 2 ) 

cTsoir) = |J (d + *V - (1 + ^)r 2 ) 

3g 2 

CTmax = CT rr (0) = O- e0 (0) = — (1 + Z^)c7 



v(r) 



64D 



(a 2 



qa 

Vmax = 64D 

> Charge sur la circonférence d'un cercle concentrique de rayon c 




> 


< > 


r 


<- 




-> 


2a 



Figure 6.12 
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Contraintes 



3 P f c 1 c 

si r < c, (T rr (r) (r ee {r) = + ^ ( ^2 ~ 21n ^ ~ 1 



3P / / r 

sic < r < a, o> r (r) = [ (1 - *0— - 2 + (1 + I — - 21n 



<roo(r) 



3P 
Aire 1 



-(1 - i/W - 2i/ + (1 + - 2 ln - 



Flèches : 

si r < c, i>(r) 



8ttD 



(£(.»-e»)(l ♦£)♦<«» ♦,>)„£ 



si c ^ r $C a, 



> Charge uniformément répartie à l'intérieur d'un cercle 
de rayon c 



o 

Oh 




Figure 6.13 



2c/ 



Contraintes : 



si r < c, cr rr (r) 



3<?c 2 



r 2 / c 2 

-(3 + v )-+(l + v ) 41n- 
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cTeeir) = f -(1 + 3v)^ + (1 + v) [ ^ - 4 ln ° 



si c ^ r ^ a, cr rr (r) 



<roo(r) 



3qc 2 



3 ^/ 2 £l_ 4 -4(l + ,)ln'" 



2^— - 4^-4(1 + ^)ln- 



Flèches : 

si r < c, »(r) = — - — + a 2 - -b 2 + (c 2 + 2r 2 ) ln - 

si c ^ a, v(r) = (f - r 2 - + \ b 2 + (c * + 2r 2) ln £ 



6.1.5 Plaques annulaires 

a) Bord extérieur simplement appuyé 

> Charge uniformément répartie q sur la plaque 



Contraintes : 




Figure 6.14 



a rr (r) 



3qa 2 
~8e r 



(l-v) 



2c 2 



+ B- 



(3 + v)- 



-4(1 + v)(a- ln -) 4 
V a) a z 
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3qa 2 ( ( 2c 2 a 1 \ 

<roo(r) = --^r (1 - v) — + + (1 + 3v)- 



8e 2 



+4(1 + v)(a — ln -) - 
\ a) a 



Flèche : 



v(r) 



q (r 4 - c 4 - 8(A + 1) (a 2 - r 2 ) c 2 + 4 (#a 4 - 2b 2 r 2 ) ln ^) 



64£> 
Constantes : 



A = ~2 2 111 

7 z — 



a 2 - c 2 c 4(1 + ^) 



1 1 (l + 3v) + (3 + v)- 1 



5 



1 - v 

> Charge P sur le bord intérieur 



1 / 6 ,z " c \ 
A{\ + v)- 0 -in-+(3 + j0 



V o 



2« 

« > 



Figure 6.1 5 



o Contraintes 



o 



cr rr (r) = 



3P 
4ttD 

3P 



2(1 + 1/) 



A-ln- )-{l-v) 1-5-r 



3P / / r\ (a 

-««=^(2(l + ,)(A-In-) + (l-,)(l-fl. 
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Flèche : 



v(r) 



Constantes 



8ttD 



((1 + A) (a 2 - r 2 ) + (Ba 2 + r 2 ) ln ^ 



B 



2 \l + v 
1 fl + v 



a 2 + c 2 



2 \\ — v 1 \a 2 + c 2 



ln- 



ln - 



b) Bord extérieur encastré 

> Charge uniformément répartie q sur la plaque 




Figure 6.1 6 



*2 



Contraintes : 



cr rr (r) 



3qa 
~^ë 2 ~ 



2c l a 1 
+ B- 



-4(1 + 1/) (A-ln-) C — 

\ aJ a 1 



3qa 2 ( (2c 2 a 2 \ r 2 



+4(1 + y) f A - ln -) 4 
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Flèche : 

V(r) = 64D 
Constantes : 

1 



[r 4 - c 4 - 8(A + 1) (a 2 - r 2 ) c 2 + 4 {Ba 4 - 2b 2 r 2 ) ln ^ 



4(1- v)a 2 + (1 + v)c 2 



(\-v) 2 + 



B 



(1 - v)a 2 + (1 + v)c 2 
> Charge P sur le bord intérieur 



/ a\ a 2 

+ ( 1 + 3^-4(1 + ^)ln-) — 

V a) c l 



(\-v)\ l + 44ln-V(l-^)4 




Figure 6.17 



Contraintes 



cr rr (r) 



3P 
AttD 



2(1 + v) ( 



A — ln 



(1 - v) 1 - B 



Flèche : 



^W=^(2(l + .)(A-Inî) + (l-v)(l-B^ 
((1 + A) (a 2 - r 2 ) + {Ba 2 + r 2 ) ln ^ 



u(r) 



8ttD 
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Constantes : 

A = -l -J (l- V )(l-^)-2(l +V ) C - 2 l n C -) 
2 (1 - v)a 2 + (1 + v)c 2 y V a ) a 2 a) 

2b 2 ( c\ 

B = - - ( l + (l + i/)ln - ) 

(1 - v)a 2 + {\ + v)c 2 V v } a) 

c) Bord intérieur simplement appuyé 

> Charge uniformément répartie q sur la plaque 




Figure 6.18 



Contraintes : 

= ^ ( (1 - v) ( 2 - 4) - (3 + - 4(1 + ^ ( A - in D) 

<r„(r) = ((1 - v) (l - B Ç) + (1 + 3^ + 4(1 + v) (a - In ^ 

Flèche : 

+4 (Ba 4 + 2b 2 a 2 ) ln ^) 

Constantes : 

c 2 a 1 / c 2 \ 

A = "2 2 ln --I7T^ (l + 3r) + (3 + r)- 

a 2 — c z c 4(1 + ^) \ a z y 
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B = — ^- f 4(1 + y) ^ 9 ln- + (3 + ^) C 



d) Bord intérieur encastré 

> Charge uniformément répartie q sur la plaque 



2a 

< ► 



Figure 6.19 



: cl 



Contraintes : 



(Trrir) -- 

a 

S Flèche : 

0) 

or 

I "(0 : 



8e 2 



3ga 2 
"8^2" 



(1-^) 2-5- -(3 



2 

, + -4(1 + */) (A -In-) 



64Z) 



(1 - v) [2 - B ^ J + (1 + 3^ + 4(1 + *0 (A - ln ^) 

r 2 -c 2 ) ( r 2 + c 2 + 8a 2 (l + A-ln^)) 

+4 (5a 4 + 2b 2 a 2 ) ln - 
v J a 



p Constantes : 
1 



o 

Oh 

^ A = ~A (\ + v)a 2 + (\ - v)c 2 



l + 3^ + (l - ^) ( 2 + -41n C ' ! ( 



Q 



B = yr- : ~ n + + ^-41n- 

(1 + + (1 - v)c L V V <2 <2 , 
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6.2 COQUES 

6.2.1 Cylindriques verticaux 

a z : contrainte suivant l'axe du cylindre 

do : contrainte perpendiculaire à l'axe du cylindre 

8 : déformation radiale 

a) Pression interne uniforme 
> Extrémités fermées 



Figure 6.20 



p r 
27 



ce 



p r 

e 



p r 
2Ë~e 



(2-v) 



> Extrémités ouvertes 



Figure 6.21 



(i 7 = 0 



pr 



8 = 
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b) Pression interne linéaire 



p(z) = Po : 



V 1 



Figure 6.22 




c) Pression externe due au vent 




Figure 6.23 



d 2 

Œ z = p—COS(p 



as = —pr sin cp 
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6.2.2 Cylindres horizontaux remplis par un liquide 

p : densité du liquide 

g : accélération de la pesanteur 

Po — Pg2r : pression maximale 

a) Extrémités simplement appuyées 




Figure 6.24 



cr z = -^z(L-z) cos0 
2e 



o-o = — (pç>r - pgr cos 6) 

e v 7 



b) Extrémités encastrées 











L 




Figure 6.25 



2e e 
+ ( — - vr 2 ] — cos 6 



&0 = - (por - pgr 2 cos/ 
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6.2.3 Coupole sphérique fermée 

e : épaisseur de la coupole 
r : rayon 



a) Charge p uniforme (poids propre) 




Figure 6.26 




b) Charge p uniforme par unité de surface projetée (neige) 



III 



III 

È 





Figure 6.27 



7?r 



pr 

as = — cos(2<p) 
2e 
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c) Pression p extérieure 









r 




Figure 6.28 



pr 


pr 




2e 



d) Charge P au sommet 



e z A 




Figure 6.29 
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6.2.4 Coupole sphérique ouverte 

e : épaisseur de la coupole 
r : rayon 

cpo : angle d'ouverture de la coupole 

a) Charge p uniforme (poids propre) 



-J' (\ 


a i 


fa 












te* 


r 







Figure 6.30 



pr(cos ç>o — cos (p) 

<p = r~2 " 

esm cp) 



pr (cos cpo — cos cp) 

(70 = = COS (p 

esm cp) 



b) Charge p uniforme par unité de surface projetée (neige) 



II II 




Figure 6.31 
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c) Pression p extérieure 




Figure 6.32 



a = EL ( i _ sm ^° 

^ 2e l sin 2 <p 



pr sin <p 0 

a 0 = — I 1 + , 2 



2e 



sm <p 



d) Charge P en bordure du sommet 





Figure 6.33 



P sin cpo 


P sin <£>o 


°> = . 2 




e sm <p 


e sm <p 
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6.2.5 Coque sphérique 

e : épaisseur de la coque 
r : rayon 



a) Pression interne p 



e r 





Figure 6.34 



(T r = 



p r 
27 



p r 
2^ 



p r 

ô = — (1 - v) 
2Ee y } 



b) Pression interne p<\ et pression externe p 2 





Figure 6.35 



_ rlp 2 + r\pi p 2 


- Pi 


^ _ rlp 2 + r\p l ( p 2 


- Pi 


' 2 (r 3 — r|) r 3 


- r 3 


2 (r 3 — r|) r 3 


- r 3 



Chapitre 7 



Formulation des éléments finis 



7.1 INTRODUCTION 

La méthode des éléments finis permet de résoudre un problème dont la solu- 
tion analytique ne peut éventuellement pas être déterminée. Elle fournit une 
solution approchée de la solution exacte. Le milieu étudié est discrétisé en 
plusieurs éléments reliés entre eux par des nœuds. La géométrie d'un élé- 
ment est caractérisée par un nombre fini de nœuds sur son périmètre. La 
résolution d'un problème par la méthode des éléments finis consiste à trou- 
ver les déplacements (par exemple les translations et les rotations) de ces 
nœuds. Le champ de déplacement en tout point est déterminé par interpola- 
tion entre les valeurs déterminées aux nœuds. L'interpolation est basée sur 
l'utilisation de fonctions déformes. 

7.2 PRINCIPE DES ÉLÉMENTS FINIS 

Pour un problème de mécanique en élasticité, le principe des travaux virtuels 
s'écrit en un point M d'un domaine H chargé par des efforts volumiques / 



7.2 Principe des éléments finis 
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et des efforts F sur une partie du bord d£l. 

s 



f(M).u(M)dV+ [ F(M).u*(M)dS 

JdVL 



dV = 0 



pour le champ de déplacement u*. La structure étudiée est discrétisée en 
un ensemble d'éléments : le maillage. Chaque élément possède n nœuds. La 
résolution du problème se base sur l'interpolation du champ de déplacement 
de chaque élément suivant la forme : 





~u(M) 




u{M) = 


v(M) 






w(M) 





<Pll <£>12 <Pl3 
<P2l <P22 <P23 
<P3l <P32 <P33 



<Pl(3n) 
Ç2{3n) 
<P3{3n) 



U\ 
V\ 
W\ 

u n 

Vn 
W n 



[<p] [u] 



Les fonctions cpij sont les fonctions de forme. Le tenseur de déformation 
s'obtient par les équations de la cinématique pour s'écrire sous la forme : 



[s(ÏÏ(M))] = [8] [u] 



(7.1) 



où [S] désigne la matrice des dérivées des fonctions de forme. 

Enfin, le tenseur des contraintes s'obtient grâce à la loi de comportement 
sous la forme : 

[a] = [E] [ô] [u] (7.2) 

Pour un matériau élastique isotrope de module d' Young E et de coefficient 
de Poisson v, 

1 



[E] = 



(1 + v){\ - 2v) 



— V 


V 


V 


0 


0 


0 


V 


1 - V 


V 


0 


0 


0 


V 


V 


1 - V 


0 


0 


0 


0 


0 


0 (1 


- 2*0/2 


0 


0 


0 


0 


0 


0 (1 


-2*0/2 


0 


0 


0 


0 


0 


0 (1 


-2*0/2 
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Finalement, le principe des travaux virtuels s'écrit : 

( I [/int] M dV+ ( [F ext ] [<p] ds) [u] 

-[uf Qf [S]' [E] [8]dV^j [u]=0 

En considérant que ce principe doit être vérifié pour tout [u] et en introdui- 
sant la matrice de raideur : 

K = / [8]' [E][Ô]dV (7.3) 

et le vecteur des efforts : 

F = / [/int] M dv + / [^ext] M dS (7.4) 
le principe des travaux virtuels se réduit au système : 

KU = F 

Le problème est résolu par inversion du système. A partir du champ de 
déplacement, le champ de déformation est calculé avec la relation (7.1), puis 
le champ de contrainte est calculé avec la relation (7.2). 

7.3 ÉTAPES DE LA RÉSOLUTION D'UN PROBLÈME 

1 . Définition de la géométrie de la structure étudiée 

2. Maillage de la structure (discrétisation en éléments) - Génération de 
la table de connectivité entre les éléments 

3. Définition du modèle associé à la structure 

4. Choix des paramètres du modèle (deux pour un matériau élastique 
isotrope) 



7.3 Étapes de la résolution d'un problème 
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5. Construction de la matrice de raideur globale par assemblage des 
matrices de raideur locales des éléments en fonction de la table de 
connectivité 

6. Construction du vecteur effort global par assemblage des vecteurs 
d'effort locaux suivant la table de connectivité 

7. Prise en compte des conditions aux limites 

8. Prise en compte du chargement 

9. Résolution du système pour obtenir les déplacements nodaux des élé- 
ments 

10. Traitement du résultat (calcul des déformations, calcul des contraintes) 

> Assemblage de la matrice de raideur et du vecteur efforts 
globaux 

Pour l'élément 1-2, les grandeurs locales s'écrivent : 



K 



K\2 K 2 2 



et F 



Fi 

F 2 



Dans le système global, et pour les nœuds / et j dans la table de connec- 
tivité : 





••• K u ■ 


■ Ky ■■■ 




Fi 


K = K + 






et F = 






••• K U ■ 


■ K n ••• 




F} 
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7.4 APPLICATION A L'ETUDE D'UNE POUTRE 
SOLLICITÉE EN FLEXION 

7.4.1 Description du problème 

Nous illustrons la méthode des éléments finis à la résolution d'un problème 
de résistance des matériaux classique, une poutre encastrée sollicitée en 
flexion. 




Le matériau est un acier de module d'élasticité E = 200 GPa. La section 
de la poutre est rectangulaire, de largeur b = 2 cm et de hauteur h = 4 cm. 
Sa longueur est i = 1 m. Le chargement est un effort ponctuel F = 1 kN, 
appliqué à l'extrémité droite de la poutre. 

Pour ce problème de flexion simple, nous utilisons une discrétisation élé- 
ments finis de type poutre d'Euler-Bernoulli. Nous comparerons la déformée 
de la solution calculée avec celle de la solution exacte qui s'écrit : 




7.4.2 Construction de la matrice de raideur locale 

L'élément utilisé pour la résolution du problème est un élément poutre à 
deux nœuds, comportant deux degrés de liberté par nœud, un déplacement 
vertical et une rotation. 



7.4 Application à l'étude d'une poutre sollicitée en flexion 
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a) Fonctions de forme 

Pour un problème de poutre en flexion, dont la déformée est un polynôme 
de degré trois, on propose la fonction d'interpolation suivante : 

v(x) = do + d\X + <22X 2 + CI3X 3 

Les constantes ^0,^1,^2 et sont obtenues en fonction de la valeur des 
degrés de liberté aux nœuds 1 et 2, soient v\,0\ et v 2 , 0 2 : 

v(0) = v x i/(0) = 0i 

v(£) = v 2 v'(£) = 0 2 

Ces quatre conditions permettent d'écrire le système de quatre équations à 
quatre inconnues, permettant de calculer les constantes aQ,a\, a 2 et : 

a 0 = v x 
ci\ = 0\ 

ao + a\£ + ci2t 2 + a^i 3 = V2 
ai + 2a 2 £ + 3a 3 f = 0 2 

soit : 

cio = V\ a\ = 0i 

v 2 -vi 20i + 0 2 vi - v 2 0i + 0 2 

G2 = 3 — 1 ^ = 2^^ + ^^ 

La fonction d'interpolation s'écrit finalement : 

U 2 - V\ 20i + 0 2 \ 2 ^ A^l ~ U2 ^ 0l + 6,2 A 

l>(x) = v\ + 0ix + 3 — — ) x + 2 — -r — + 



i 1 i \ P i 1 
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Sous forme matricielle, nous voyons apparaître les fonctions de formes <p 



v(x) 
0(x) 



3x 2 2x 3 2x 2 x 3 3x 2 2x 3 x 2 x 3 

l r + ~f x ~~f + ¥l 2 ~~l 3 ~ + ¥ 

6x 6x 4x 3x 2 6x 6x 2 2x 3x 2 

■ ~l 2+ ! 3 ~ ~T + 1 2 ~1 2 ~1 3 ~~T + 1 2 





Vl 








A 



Le calcul de la déformation repose sur l'hypothèse d'Euler-Bernoulli : 

d 2 v(x) 



[e] = en 



dx 2 



-y 



et l'écriture sous forme matricielle donne : 



M = -y 



6 12x 4 6x 6 12x 2 6x 



^2 



La relation de comportement est réduite à : 

[a] = [E] [s] = E [s] 

b) Construction de la matrice de raideur locale 

Les termes de la matrice de raideur locale se calculent à partir de l'expres- 
sion (7.3) : 



K = [ [S]< 



[E][ô]dV 



K 2 2 K 2 3 K 2 4 

(Sym) K 33 K 34 

K44 



7.4 Application à l'étude d'une poutre sollicitée en flexion 
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Les différents termes étant : 
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-b/2 J-h/2 . 

r f" /2 f h/2 2, /"V16 48x 

E LL/ dydz L 



El 



\6x 24x 2 12x 31 

— ^ 37- + 



36x 2 

£ 4 



4EI 



dx 



-K 



12 



6EI 



rb/2 rh/2 ,i 2 / 4 gx / 

J-t/2 J-H/2 Jo y { £ + £'){ 



-b/2 J-h/2 

E UL/ dydz l Kë~- 

£7 



2 6x 
36x 2 



8jc 18jc 2 12jc 3 



2£7 



dxdydz 
dx 



o 



12£7 



6£7 



£ 2 

/b/2 çh/2 ni 
\ I y 
-b/2 J-h/2 JO 

.b/2 .h/2 ,i / 

' / / y 2 dydz 

J -b/2 J-h/2 JO \ 



~~£*% ' dxdydz 



El 



-b/2 J-h/2 

4x 12* 2 12x 3 



24x 36x 2 \ , 
— ^- + — 77- dx 



4EI 



7.4 Application à l'étude d'une poutre sollicitée en flexion 



La matrice de raideur pour l'élément poutre est : 



K„ 



4 

e 2 - 
4 

6 EIu 

4 



6^1 

4 



'7 





6 ? 1 




4 
4 
4 


k 

k J 



7.4.3 Implantation et résolution dans Matlab 

Le programme permettant la résolution du problème est le suivant 



°/o Paramètre matériau 
E = 200. e9 ; % [Pa] 

°/ 0 Longueur de la barre 
lb = 1. ; °/ 0 [m] 

°/ 0 Caractéristique de la section 
b = 0.02 ; °/ 0 Largeur [m] 
h = 0.04 ; °/ 0 Hauteur [m] 
A = b*h ; °/ 0 Section [nT2] 

I = b*h~3/12. ; % Moment quadratique [m" 4] 
nel = 5 ; °/ 0 nombre d' éléments 

ndof = 2 ; °/ 0 nombre de degrés de liberté par noeud 
ntotdof = ndof*(nel+l) ; °/ 0 nombre total 
de degrés de liberté du système 

1 = lb/nel ; °/ 0 Longueur d'un élément 
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°/ 0 Conditions aux limites 
ddlfixe(l) = 1 ; valfixe(l) = 0. ; 
ddlfixe(2) = 2 ; valfixe(2) =0. ; 

°/ 0 Déclaration de K et F 
vecforce = zéros (ntotdof , 1) ; 
kglobal = zéros (ntotdof , ntotdof ) ; 

°/ 0 Chargement 

vecforce (ntotdof -1) = l.e3 ; °/ 0 [N] 

°/ 0 Construction de la matrice de raideur globale 
for iel=l:nel 

°/ 0 Définition de la matrice de raideur locale 
klocal = [12.*E*I/1~3 6.*E*I/1~2 
-12.*E*I/1~3 6.*E*I/1~2; 

6.*E*I/1~2 4.*E*I/1 -6.*E*I/1~2 2.*E*I/1; 
-12.*E*I/1~3 -6.*E*I/1~2 12.*E*I/1~3 
-6.*E*I/1~2; 
6.*E*I/1~2 2.*E*I/1 -6.*E*I/1~2 4.*E*I/1] ; 
°/ 0 Assemblage de la matrice globale 
for idof=l:2*ndof 
for jdof=l:2*ndof 
igl = ndof *(iel-l)+idof ; 
jgl = ndof *(iel-l)+jdof ; 
kglobal (igl, jgl) = kglobal (igl, jgl) 
+ klocal(idof , jdof ) ; 

end 

end 

end 

7o Prise en compte des conditions aux limites 
for idof =1 : length(ddlf ixe) 



7.4 Application à l'étude d'une poutre sollicitée en flexion 



c = ddlf ixe(idof ) ; 
for jdof =1 : size(kglobal) 
kglobal(c, jdof ) =0. ; 

end 

kglobal (c,c) =1. ; 
vecforce(c) = valf ixe (idof ) ; 

end 

°/ 0 Inversion du système 
vecdisp = kglobal\vecf orce ; 

°/ 0 Déplacement calcule 

disp( [' Déformée max calculée : 

' ,num2str (vecdisp (ntotdof-1) )] ) 

°/ 0 Déplacement théorique 

disp( [' Déformée max théorique : 

' ,num2str(vecforce(ntotdof-l)*lb~3/(3.*E*I))] ) 

°/ 0 Construction de la solution approchée 

abscisse = zéros (nel+1 , 1) ; déformée = zéros (nel+1 , 1) 

i = 0 ; 

for idof =1 :ndof :ntotdof 
i = i+1 ; 

abscisse(i) = (i-l)*l ; 
deformee(i) = vecdisp(idof ) ; 

end 

°/ 0 Construction de la solution exacte 
npref = 100 ; 

abscisse_ref = zéros (npref , 1) ; 
deformee_ref = zéros (npref , 1) ; 
for idof=l: npref 

abscisse_ref (idof ) = lb/(npref-l) *(idof-l) ; 
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def ormee_ref (idof ) = vecf orce (ntotdof -1) * . . . 

(lb*abscisse_ref (idof ) ~2/2-abscisse_ref (idof) "3/6)/ 
(E*I); 
end 

plot (abscisse_ref ,def ormee_ref , abscisse , def ormee) 
legend ( ; exacte ; , ' approchée ; , 2) 
xlabel ( ; abscisse (m) ' ) 
ylabel( ; def ormee (m) ' ) ; 

Comme attendu, la solution éléments finis aux nœuds de chaque élément 
est exacte : la fonction d'interpolation est en effet du degré polynomial de la 
déformée. La comparaison entre la solution exacte et la solution numérique 
obtenue avec cinq éléments est représentée ci-dessous. 



0.016 




abscisse (m) 

Figure 7.1 



7.5 Éléments isoparamétriques 
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7.5 ÉLÉMENTS ISOPARAMÉTRIQUES 

Les éléments isoparamétriques reposent sur l'utilisation d'un change- 
ment de repère entre deux systèmes de coordonnées. Ils permettent 
l'écriture unique des fonctions de forme pour un type d'élément. Les 
fonctions de formes s'écrivent dans le repère naturel 77, £), le pro- 
blème de référence s 'écrivant dans le repère physique (x,y,z). La 
transformation géométrique pour passer d'un repère à un autre uti- 
lise les fonctions de forme de l'élément. Ainsi, pour un élément à n 
nœuds, 

n 

i=l 

n 

i = \ 

n 

i = \ 
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7.6 FONCTIONS DE FORME DES ÉLÉMENTS 
ISOPARAMÉTRIQUES COURANTS 

7.6.1 Élément barre à deux nœuds 



V 

> 






1 




2 


• 




— — • > 


(-1,0) 




(1,0) f 


Figure 7.3 


*>!(£, 17) =^d-f) 





7.6.2 Élément barre à trois nœuds 



(-1,0) 



(1,0) 



Figure 7.4 



6 Fonctions de forme des éléments isoparamétriques courants 169 
6.3 Élément triangulaire à trois nœuds 




6.4 Élément triangulaire à six nœuds 




(1,0) 



ç>i(f,i7) = (l-f-7j)(l-2f-2 î? ) 


? 2 (£,T7) = £(2f-l) 


<P3(£, t?) = t?(2t? - 1) 


= 4^(1-^-7?) 


<P5(Ç, v) = Hv 


<P6(^, T?) = 4T?(1 - ^ - 7?) 
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7.6.5 Élément quadrangulaire à quatre nœuds 





71 A 




(1,1 




4 


3 


-> 




1 


2 




(-1,- 


-1) 







1) 





<P2(£,v) = \a+ex\-v) 







7.6.6 Élément quadrangulaire à huit nœuds 





71 > 




(1,1) 




4 


7 3 
6 






1 


5 2 




(-1,- 


-1) 







7.6 Fonctions de forme des éléments isoparamétriques courants 
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<Pi(£,v) = \a-SM-vX.-i-S-v) 

4 


<P2(Ï,V)= J(l + f)(l-r/)(-l + f-r?) 


<P3&V)= i(l + £)(l + »7X-l + £ + 


Mt,v) = i(W)(i + *?)(- W + rç) 




P6^,'))=^ + ai-'/ 2 ) 




^,'?)=^(l-f)(l-'? 2 ) 



7.6.7 Élément quadrangulaire à neuf nœuds 













4 


7 


3 






( 8 




6 ( 
















1 


5 


2 




(-1, 


-1) 









«P2a,î?) = ^(f 2 + I r )(î? 2 -T7) 


<P3(€,v) = \(€ 2 + t)(v 2 + v) 




V5(è,v)=\(\-è 2 )(v 2 -v) 


«P6(l : ,r ? )=i(f 2 + ^)(l-r 7 2 ) 




?>8(f,T?)=^(f 2 -^(l-i7 2 ) 


« P9 (^7 ? ) = (l-f 2 )(l- 7 7 2 ) 





Chapitre 8 



Instabilité des structures 



Dans ce chapitre sont exposés les principaux résultats concernant le calcul 
mécanique de structures élancées soumises à des risques d'instabilités. Pour 
de tels éléments, les classiques équations de la RDM relatives à la méca- 
nique des milieux continus en petites transformations ne s'appliquent plus. 
Il faut dorénavant tenir compte de la déformation de la structure pour évaluer 
les efforts internes à cette dernière. 



8.1 INSTABILITÉ DE POUTRES 
8.1.1 Poutre d'Euler 

Soit une poutre bi-articulée de longueur /, de section S et de moment qua- 
dratique /. La colonne est chargée par une force P, considérée positive 
en compression. Nous supposerons que la colonne est parfaitement alignée 
avec l'axe de la charge avant que cette dernière ne soit appliquée. La plus 
petite charge pour laquelle l'instabilité élastique est détectée est dénommée 



8.1 Instabilité de poutres 
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la charge critique d'Euler P E et a pour expression : 

tt 2 EI 

La déformée correspondante forme une demi sinusoïde, donnée à un facteur 
amplificateur près q : 

v(x) = q sin(77x)// 

en notant i = y/TJS, le rayon de giration de la section, on peut définir 
l'élancement À de la poutre par la relation À = l/I. On peut ainsi obtenir la 
valeur de la contrainte critique de flambement : 

tt 2 E 
(Te = — — 



\v(x) 



Figure 8.1 



La charge critique d'Euler correspond au plus petit effort possible per- 
mettant de maintenir l'équilibre sous la déformée v(x). Le flambement 
d'une telle poutre peut être obtenu pour des charges critiques P CYn plus éle- 
vées correspondant à des harmoniques v n (x) différentes. La charge critique 
de flambement du mode n ainsi que son allure sont données par : 

n 2 7r 2 EI 

Pcr n = p et v n (x) = q n sin(/î77x)// 
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8.1.2 Solutions générales des poutres comprimées 

Du fait des conditions limites particulières liées à la résolution du problème 
de la poutre d'Euler au paragraphe précédent, les efforts internes à la struc- 
ture pouvaient être déterminés sans prise en comptes des efforts de liaison. 
Dans le cadre le plus générale, ces derniers doivent être pris en compte, et 
l'équation différentielle suivante doit être résolue, avec les effets du second 
ordre dans l'équilibre d'un tronçon de poutre : 

(EIv"{x))" + (Pv'(v))' = 0 

La solution générale pour E, I et P constants est donnée par : 

v(x) = A sin kx + B cos kx + Cx + D 

A, B , C et D sont des constantes à déterminer en fonction des conditions 
aux limites de la pièce étudiée. Les types de liaison les plus fréquemment 
rencontrés sont, pour un appui O : 

- L'encastrement :L>o=0eti/ o =0 

- La rotule : vo = 0, M 0 = 0 et uo =0 

- L'appui simple : vo = 0 et M 0 = 0 

- L'appui continu de raideur k 

- Libre : V () = 0 et M () = 0 

- L'encastrement partiel, induisant une raideur en rotation : &o 

8.1.3 Solutions particulières pour des poutres de section 
constante 

Nous appelerons /o la longueur de flambement, correspondant à la longueur 
qu'il faudrait appliquer à la solution de la poutre d'Euler pour obtenir les 
mêmes charges critiques de flambement pour tout autre type de conditions 
limites. 
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Cas de colonnes supportant une charge axiale en tête : 



Type de liaisons 


Schémas 


Charges critiques 


Longueur de 
flambement 


Appui simple - 
Rotule 




\ 

\ 

I 
1 
1 
1 
]l 


7T 2 EI 

Pcr — — — 

£ 2 


to=£ 


Appui simple - 
Encastrement 




\ 
\ 
1 
/ 

/ 


2 ri 

cr (0.699^ 2 


£ Q = 0.699^ 


Encastrement - 
Encastrement 




) 

/ 


_^EI 

Pcr ~ fi 


£ 0 = 1/2 


Encastrement - 
Libre 




/ 


_^EI 


£ 0 = 2£ 


Encastrement - 
Encastrement 

et v 0 / 0 




TT 2 EI 


£o = £ 


Encastrement 
partiel - libre 


/ 

A 


/ 

1 


i. 2 ci 

cr ~ 4k£ 2 + tt 2 EI£ 


1 ici 


Rotules - 
Appui continu 


I 77 V El 


Pcr = 

7T 2 EI ( 2 ki A \ 

l 2 [ m + m 2 ^El) 


£/\m 2 + k f A 
1 V m 2 7r 4 EI 
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Cas de colonnes tenues par des encastrements partiels, les raideurs en 
rotation sont caractérisées par des paramètres &a et : 

- Pour des noeuds d'extrémités non déplaçables : 



A W 



Nous posons 



1 1 

k A Kg 



1 , El 

—— et £ = — 
k A k B l 



etK 



I 
I 
I 



2>k 2 p + ks + Qi32 
3k 2 p + 0.1ks +0.16 



/ 



,ir 2 EI 



Figure 8.2 

- Pour des noeuds d'extrémités déplaçables : 



1 1 7 El 
s = — + — et £ = — 



I 
I 
I 
I 
I 
I 
I 



k k 
Pa = 7— et p fl = — 

0.4 + p fi 0.4 + p A 



0.8 + ^ 



0.8 + ^ 



/ 
/ 

rB A 

Figure 8.3 



2 77 p A a A 2 77 

= ol a + , K B = a B + 



et ^ = Max(X A et K B ) 
tt 2 EI 



4K1 2 



et /o = 2/vT 



8.1 Instabilité de poutres M7 



Cas de colonnes supportant des charges sur toute leur hauteur : 



Type de liaisons 


Schémas 


Charges critiques 


Longueur de 
flambement 


Encastrement - 
Libre 


t 


J 


\\\ 


/ .o37T tl 


4 = 1.12^ 


Encastrement - 
Libre 


\ fi\ 

j Xi 
WWW 


Pcr = 
n 2 EIY;Pi 


4 = 

/EMl-cos^) 


4**£P ; (1 -cosîS) 


V E^/ 


Encastrement - 
Libre 


t 




P 

w\ 


^ ^ + 1-20 


/ 4P 

£ 0 = £\ +1.26 


Rotule - Rotule 


t 


î * 


p 

'p 


_ 2^0 





8.1.4 Prise en compte d'un défaut initial 

Le cas des colonnes parfaites étudiées jusqu'à présent est un modèle 
idéalisé. En réalité, différentes sortes d'imperfections inévitables doivent 
être considérées. Par exemple, des poteaux peuvent être soumis à des 
charges latérales parasites ou présenter une géométrie courbe induisant une 
excentricité initiale. Dans ces deux cas il est important de pouvoir estimer 
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l'amplification de ce défaut initial à l'approche de la valeur critique de la 
charge de flambement. Si on assimile un défaut initial à une courbure de la 
colonne ayant pour déplacement maximal vq (ou moment parasite maximal 
M 0 ), le déplacement (ou moment) atteint en ce même point, fonction de la 
charge appliquée P est amplifié du paramètre jul : 

1 

M " l-P/P C r 



a: 




déplacement transversal v 
Figure 8.4 



8.2 CALCUL DES MOMENTS DANS UNE POUTRE 
COMPRIMÉE FLÉCHIE 

En considérant une poutre comprimée, sujette au risque de flambement, et 
supportant des charges latérales induisant une flexion dans le plan de flam- 
bement, l'analyse des contraintes requiert de prendre en compte l'amplifi- 
cation des efforts due à l'instabilité. 

Le moment fléchissant maximal et déplacement transversal maximal 
dus uniquement aux charges latérales sont dénommés respectivement M 0 

et Vq. 



8.3 Déversement latéral de poutres 
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Poutre sur deux appuis simples d'extrémité supportant une charge 
transversale symétrique : le moment est maximal en son centre et une 
bonne approximation de ce dernier est obtenu en considérant : 



M = M 0 1 + 



7T 2 U 



V(l - U) 



P 7T 2 EI M 0 l 2 

Avec u = — , P cr = — — et ^ 



P c / l 2 EIvo 

Poutre en console 



M = M 0 I 1 + 



? 

TT U 



4^(1 - U) 



P 7T 2 EI Mol 2 

Avec u = —, P cr = -rjï- et * = -f- 
P cr M 1 EIvo 



8.3 DÉVERSEMENT LATÉRAL DE POUTRES 

Le déversale latéral d'une poutre est un phénomène d'instabilité élastique 
consistant en un flambement accompagné de torsion. Ce phénomène se pro- 
duit pour des poutres étroites, non maintenues latéralement, pour des valeurs 
de contraintes de flexion bien inférieures à la limite d'élasticité du matériau 
considéré. 



8.3.1 Déversement latéral de poutres à section rectangulaire 

En considérant une poutre en flexion, il est possible d'obtenir la valeur limite 
des moments maximaux critiques menant au déversement. Le tableau ci- 
après récapitule les diverses conditions classiquement recontrées. On note 
G, le module de cisaillement et /, le moment d'inertie en rotation autour 
de ~z . 
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Description 


Schémas 


Moments critiques 


Poutre sur appuis simples 
soumise à un moment 
constant 


M x m 

v 

T J—l 


a/j ^ /ci r~ i 
Mer = —ytlyKjJ 


Poutre bi-rotulée soumise 
à un moment constant 


fP- Ù M 

r* a 1 

1 t \ 


M cr =^-^EI y GJ 


Poutre sur appuis simples 
soumise à une charge 
constante 


x p 

H H H [ 

r — ï — i 


M Cr =^^^EIyGJ 


Poutre sur appuis simples 
soumise à une charge 
concentrée 


à 


-3 \* 

* — t — 1 




Poutre en console soumise 
à une charge constante 


www 


-3 p 

' t 1 


M cr = —^y/ElyGJ 


Poutre encastrée soumise 
à une charge concentrée 


/ 


~x 

^ i p 

' e 1 


1 28tt 
M cr = —^y/ElyGJ 



8.3.2 Déversement latéral de poutres à section en I 

Dans le cas particulier des profilés à section ouverte, ou en I, les conditions 
d'attache aux extrémités peuvent limiter le gauchissement de la section et 
donc engendrer des augmentations de raideur en cisaillement. Les simples 
résultats précédents doivent être modifiés pour prendre en compte ces varia- 
tions. 



8.4 Instabilité et voilement de plaques 
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Description 


Schémas 


Moments 
critiques 


<*1 


a 2 


a 3 


Poutre sur appuis 
simples soumise 
à un moment 
constant 


r — i — i 


Mer = jy/ElyGJ 


1 


1.83 


0 


Poutre sur appuis 
simple soumise 
à une charge 
concentrée 


13 f 

r e 1 


M cr = jy/ElyGJ 


1 


1.35 


0.55 


Poutre bi-rotulée 
soumise à une 
charge concentrée 


13 i p 

1' e 1 


M cr = jy/ElyGJ 


0.5 


1.07 


0.42 


Poutre bi-encastrée 
à déplacement axial 
libre 


±3 f c 

* — 1 — r 


M cr = jy/ElyGJ 


1 


1.7 


1.42 



Avec a = — 77 



I 7T 2 EI M , 2 v tf 3 77 /E/û, 



et en 



notant /z la hauteur du profilé, I œ = -+-h . J est le moment d'inertie polaire. 



8.4 INSTABILITÉ ET VOILEMENT DE PLAQUES 

L'instabilité par voilement se rencontre pour des éléments de structures plan 
élancés de type plaque. De façon analogue au flambement des poutres, il 
convient d'étudier l'équilibre d'une plaque en position déformée. Afin de 
déterminer les équations d'équilibre d'un élément de plaque soumise à une 
compression dans son plan, il convient de faire l'hypothèse d'une surface 
« presque plane ». De telles hypothèses nous permettent de prendre les axes 
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du repère général comme ligne de coordonnées portant les abscisses curvi- 
lignes ~5t et ~y de la plaque. En projetant les équations vectorielles d'équi- 
libre d'un élément de plaque dans le repère lié à la plaque, nous obtenons 
cinq équations d'équilibre dont une seule sera utile pour notre cas. En intro- 
duisant les expressions des efforts tranchants en fonction de la déformée 
transversale de la plaqué notée co, nous obtenons l'équation de Lagrange 
généralisée. 



a , T d 2 co , T d 2 co _ d 2 co 
D^&co + N x — + N y — + 2N xy —— 
ox z oy l oxoy 



0 



Pour une plaque d'épaisseur h, la rigidité de flexion de l'élément vaut : 
D = Eh 3 / 12 (l — v 1 ). Nous traiterons d'une plaque simplement appuyée 
sur ces quatre côtés, soumises à un effort normal le long du plus grand 
coté a. 




Figure 8.5 

En dénommant les différents modes d'instabilité possible par le nombre 
de demi-ondes metn transversales de cloquage développées respectivement 
dans les directions x et ~y, l'allure de la déformée transversale (à un facteur 
d'amplification q près) est de la forme : 

. mirx . niry 

co xy = q mn sm sm — — 

a b 

La charge fondamentale de flambement (la plus petite) est obtenue avec n = 
1 et m = a/b et vaut F cr = Ai^Djb 1 . Si les dimensions de la plaque sont 
fixées, m n'est pas un entier. Nous pouvons tracer sur un schéma l'évolution 
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de la charge critique en fonction du rapport a/b pour différentes valeurs 
de m : 




Figure 8.6 



En fonction des dispositions d'appuis le long des faces latérales, les 
efforts critiques de voilement prennent différentes valeurs. 



Description 


Schémas 


Efforts critiques 


Plaque sur appuis simples 


a 


Fer = 4ir 2 D/b 2 


Plaque encastrée-appuis simples 


a 


F cr = 5A27T 2 D/b 2 


Plaque bi-encastrée 


' a 


F cr = 6.977r 2 D/b 2 


Plaque encastrée-libre 


a 


F cr = 1.2877- 2 D/b 2 


Plaque appuis simples-libre 


a 


F cr = 0A257T 2 D/b 2 
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8.5 FLAMBEMENT DE STRUCTURES NON PLANES 
INITIALEMENT 

8.5.1 Flambement d'arc et d'anneaux 

Considérant des pièces non planes supportant une pression extérieure p, 
les résultats ci-après récapitulent quelques résultats importants de pression 
critique avant flambement. 

- Arc circulaire simplement appuyé : 

a représente l'ouverture angulaire de l'arc de rayon R considéré. 



477 2 

a 2 



1 R3 



Arc circulaire encastré à ses extrémités : 



(k 2 ~ 1) 



El 



avec 



a 


60° 


120° 


180° 


240° 


300° 


360° 


k 


8.62 


4.37 


3 


2.36 


2.07 


2 



8.5.2 Flambement de tubes minces 

Le tube est d'épaisseur e, de longueur L et de rayon R. On note : 



R 



La contrainte critique de flambement a„ est donnée ci-après. 



8.5 Flambement de structures non planes initialement 



185 



Tube mince en compression simplement appuyé à ses extrémités 

P 



-V 




7T^ 



(j C y — kç 



7T Z E 



12 (1 - v 2 ) (l/e) 2 



Figure 8.7 

Tube mince en compression encastré à ses extrémités 

4 + 3/3 2 




7T 2 E 



12(1-p 2 ) (l/e) 2 



Figure 8.8 

- Tube mince sous pression externe 

La pression critique de flambement du tube soumis à une pression uniforme 
vaut : 



Chapitre 9 



Calcul non-linéaire, analyse 
limite, plasticité 



9.1 INTRODUCTION 

Les problèmes précédemment traités étudiaient le comportement élastique 
des matériaux. Cependant les lois de comportement élastiques linéaires ne 
sont pas toujours valables. Lors de sollicitations exceptionnelles ou acciden- 
telles, des contraintes supérieures à la limite élastique du matériau peuvent 
se produire et entraîner des déformations permanentes tandis que la struc- 
ture peut continuer à résister sans s'effondrer. Nous étudions alors la stabilité 
globale de la structure suite à l'apparition locale de la plasticité. 
Cette étude de stabilité s'appelle analyse limite en RDM. 



9.2 Modèles de comportement des matériaux 
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9.2 MODÈLES DE COMPORTEMENT 
DES MATÉRIAUX 

Si l'on trace les courbes de comportement cr-s pour un matériau ductile on 
obtient une courbe de même type que celle de la figure 9.1 (a). A partir de 
cette courbe de comportement, on peut adopter plusieurs modélisation : 

- figure 9.1 (b) : modèle élasto-plastique avec écrouissage, 

- figure 9.1 (c) : modèle élasto-plastique parfait, 

- figure 9.1 (d) : modèle rigide plastique ou plastique parfait. 




I 



(c) 



(d) 



Figure 9.1 Lois de comportement : matériau ductile (a), modèle élasto-plastique 
écrouissable (b), modèle élasto-plastique parfait (c), modèle rigide plastique (d) 



9.3 PLASTIFICATION EN FLEXION : NOTION DE 
MOMENT PLASTIQUE ET ROTULE PLASTIQUE 

9.3.1 Hypothèses 

Nous considérerons par la suite que : 
- le matériau est isotrope, 
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- le comportement est élasto-plastique parfait, 

- l'hypothèse de Bernoulli est vérifiée. 

9.3.2 Section symétrique 





Figure 9.2 Poutre en flexion de section symétrique 



(9.1) 



Considérons la poutre de la figure 9.2 de section symétrique majoritaire- 
ment soumise à de la flexion. En élasticité, nous pouvons écrire la relation 
suivante entre la contrainte et le moment fléchissant : 

M M 
<j = — v = — 
/ W 

avec W = I/v le module de flexion de la section considérée. 

Si nous continuons à augmenter le moment de flexion, la limite élastique 
va alors être atteinte en premier pour y = v. Nous obtenons alors les phases 
de comportement de la figure 9.3 : 

- phase élastique : comportement « classique » d'une section en flexion 

- État 1 : limite élastique. Le moment sollicitant est alors égal au 
moment maximum élastique M e . On a alors cr max = a e = M e /W 
d'où M e = Wa e . La courbure élastique maximum vaut alors 
Xe = /(El) et la déformation élastique s e = \e-v 

- État 2 : plastification progressive de la section (M > M e ). Dû à la 
symétrie de la section et de la courbe de comportement du matériau, 
l'axe neutre reste à mi-hauteur. On peut montrer que Xejx — 2y e /h 
ce qui implique que la courbure augmente tandis que y e diminue. 
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L 

Phase élastique 



7 




Etatl 




8 -de /° 

> — s^- — ► 



y e =o 



État 2 État 3 

Figure 9.3 Différentes phases de comportement d'une section en flexion 



- État 3 : plastification complète de la section. Formation d'une rotule 
plastique, lim y = oo. A cet état correspond le moment plastique 

y e ->0 

noté M p i . On définit le module plastique W p i = M p \ /cr e = 2M X avec 
M x le moment statique de la moitié de la section / x. 

Pour une section rectangulaire b x h on obtient : 

W = bh 2 /6 et W p i = bh 2 /4. 

Soit M e — Wcr e le moment élastique maximum et M p[ = W p icr e le moment 
plastique. On définit le gain de résistance par plastification : 



G =100(^-1 



ce qui donne pour une section rectangulaire, G = 50 %. 



(9.2) 
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9.4 ANALYSE LIMITE D'UN SYSTÈME DE POUTRES 

9.4.1 Enjeux 

Soit un système de poutres soumis au chargement extérieur T ext . Sous l'hy- 
pothèse d'un comportement plastique parfait du matériau constitutif de la 
structure, celle-ci sera mise en ruine pour un chargement T cr u = Kit^ext 
par l'apparition d'un nombre suffisant de rotules plastiques transformant 
la structure en mécanisme. On aura donc un accroissement proportionnel 
de toutes les charges par le facteur de charge À jusqu'au facteur de charge 
ultime Kit . Pour déterminer celui-ci on a deux approches : la méthode STA- 
TIQUE et la méthode CINÉMATIQUE. 

9.4.2 Théorème statique 

À tout état de contrainte statiquement admissible (\M\ ^ M p i partout) véri- 
fiant les conditions cinématiques (CL. et moment M p \ de même sens que les 
rotations plastiques au droit des rotules) correspond un facteur de charge : 

Ktat ^ Kit (9.3) 

> Application 

Soit la structure de la figure 9.4. 

• On calcule le moment M dans la structure et on trace le diagramme 
associé. 

• On suppose que ±M p i est atteint en un point de la structure (point de 
moment maximum) pour A = Ai. On a donc Aig/ 2 /12 = M pt soit 

Ai = \2M pl /ql 2 . 

• On place deux rotules en A et C (voir figure 9.5) qui reprennent alors un 
moment égal à — M p \ et on vérifie que \M\ ^ \M p \ | pour toute la structure. 
Ici, M B = \ql 2 /8 - M p i avec : 

A > Ai => A<?/ 2 /8 > Ai<?/ 2 /8 = 3M pl /2 > M pl . 

on a donc bien M B > 0 et \M B \ ^ M pi . 



9.4 Analyse limite d'un système de poutres 
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Figure 9.4 Poutre bi-encastrée en flexion 




Figure 9.5 Poutre bi-encastrée en flexion avec deux rotules plastiques en A et C 
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• On continue jusqu'à la ruine de la structure qui correspond à la forma- 
tion de n rotules où n = H + 1. Ici, À croît jusqu'à avoir M B = M p \. Ce qui 
donne À stat ql 2 /8 — M p \ = M p \ soit À stat = \6M p i/ql 2 . On obtient alors le 
mécanisme de la figure 9.6 



M, 



Xq 



Sa 

I < 



M pl 



m pi\b/ m pi 



i 



Figure 9.6 Poutre bi-encastrée en flexion avec trois rotules plastiques en A, B et C 



9.4.3 Théorème cinématique 

A un mécanisme de ruine arbitraire, où tous les éléments non plastifiés sont 
supposés rigides (entre deux rotules successives), correspond un facteur de 
charge 



^cin ^ ^ult 



(9.4) 



> Application : 

Soit la structure de la figure 9.7. 



API 



X 



Figure 9.7 Poutre simplement appuyée en flexion 



9.4 Analyse limite d'un système de poutres 
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• On imagine un mécanisme de ruine, par exemple celui de la figure 9.8. 




Figure 9.8 Poutre simplement appuyée en flexion 



• Pour chacun des mécanismes on applique le PTV à des champs rigidi- 
fiants par morceaux plastiquement admissibles (\M p i\ ^ M p i). Dans notre 
exemple, cela donne : 

-2M p iO + ÀPO 1 - = 0 soit A = (9.5) 

• On aurait pu placer la rotule en une section quelconque — > À/ . Le théo- 
rème cinématique nous dit alors que X cin ^ min {Xi}. 

a) Points particuliers de la méthode cinématique 

Sections potentiellement critiques : 

- point d'application des charges ponctuelles 

- point de moment maximum pour les charges réparties 

- assemblage 

- encastrement 

Mécanismes de ruine possibles : 

- mécanisme de poutre 

- mécanisme de panneau 

- mécanisme des nœuds 
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Soit la structure de la figure 9.9. 



Figure 9.9 Structure de poutres 

Les figures suivante illustrent les différents modes de ruine que l'on peut 
rencontrer. 





Figure 9.10 Mécanisme de poutre : ruine partielle {n < H + 1) 



Figure 9.11 Mécanisme de panneau 




Figure 9.12 Mécanisme de nœud : ruine partielle {n < H + 1) 



Chapitre 10 



Dynamique et vibrations 



Pour un problème dynamique, le chargement et les paramètres de réponse 
sont fonction du temps. La principale caractéristique d'un problème dyna- 
mique est la présence d'efforts d'inertie non-négligeables qui, d'après le 
principe de d'Alembert, s'opposent au mouvement imposé par le charge- 
ment appliqué. 

Une charge est dynamique si sa norme, sa direction ou son point d'ap- 
plication varie avec le temps. Si l'évolution de ces grandeurs est prescrite 
de manière unique, l'analyse est dite déterministe. Au contraire, si la varia- 
tion du chargement en fonction du temps n'est pas connue et ne peut être 
définie qu'en termes statistiques, le chargement est dit aléatoire. Parmi les 
différentes charges dynamiques, on peut distinguer : 

- les charges périodiques de période T (p(t + T) = p(t)), parmi les- 
quelles nous pouvons distinguer les charges harmoniques simples et 
les charges périodiques quelconques, 

- les charges non périodiques : charges impulsionnelles, charges arbi- 
traires de longues durées (séisme, vent, houle, ...). 
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10.1 SYSTÈME À 1 DEGRÉ DE LIBERTÉ 

En dynamique, on appelle le degré de liberté d'un système la somme du 
nombre de mouvement possible de chaque masse du-dit système. Toute 
la difficulté du modélisateur réside donc dans l'étape de discétisation per- 
mettant de choisir le nombre minimum de masses ponctuelles décrivant au 
mieux la réponse dynamique de l'ensemble. 



10.1.1 Équation du mouvement 

Soit un système masse-ressort-amortisseur soumis à la gravité et à un char- 
gement extérieur fluctuant dans le temps p(t). Les différentes forces agissant 
sur un tel système peuvent être associées à : 

- un ressort supposé linéaire dont le déplacement u(t) développe une 
force de rappel élastique f e (t), 

- un amortisseur, supposé linéaire, produit une force dissipative vis- 
queuse proportionnelle à la vitesse fd(t), 

- une accélération subie par la masse m, dans un référentiel fixe, se tra- 
duit par une force d'inertie (selon le principe de d'Alembert) fi(t), 

- des efforts extérieurs f ext (t) appliqués à la masse m proviennent de 
deux origines : un effort statique dû au poids propre mg et une sollici- 
tation dynamique p(t). 



www 

I fe(t) = -ku(t) 

k\ [ylc fé(t) = -cù(t) 

\ fi(t) = -mu(t) 



\m~g B I \ p{f) 
Figure 10.1 



10.1 Système à 1 degré de liberté 
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L'équation d'équilibre dynamique de la masse m s'écrit en invoquant le 
principe de d'Alembert par la prise en compte des efforts d'inertie sous 
forme d'efforts extérieurs : ^ F(t) = 0, soit : 

ku(t) + cù(t) + mu(t) = mg + p(t) 

En décomposant le déplacement u(t) en un terme statique u s = mg/k 
et un terme dynamique u d (t), il s'avère que les composantes statiques s'éli- 
minent entre elles dans l'équation du mouvement précédente pour ne conser- 
ver que la partie dynamique du déplacement : 

ku d (t) + cù d (t) + mû d {t) = p(t) 

Pour des systèmes linéaires, il y a découplage entre les effets statiques et 
dynamiques. Un problème de dynamique linéaire ne devra pas comporter 
de composantes statiques aussi bien en effort qu'en déplacement. Seule sera 
prise en compte la partie dynamique du déplacement assimilée à une pertur- 
bation dynamique autour de la position d'équilibre statique. Dans la suite, 
dans un souci de simplicité nous noterons u d (t) = u(t). 

On peut différentier l'étude d'un système dynamique en fonction de dif- 
férents régimes vibratoires : libre, forcé ou permanent. 

10.1.2 Le régime libre 

Un régime libre correspond à la solution générale de l'équation différen- 
tielle sans second membre (pour p(t) = 0). Les conditions intiales à t = 0 
sont un déplacement et une vitesse initiale notée uç> et ùo. On note la pul- 
sation co 2 = k/m, l'amortissement critique c c = 2m co = iVkm = 2k/(o 
et l'amortissement relatif è, — c/c c . En recherchant des solutions sous la 
forme u = C z st , l'équation caractéristique du système est : 



s 2 + 2Çcos + co 2 = 0 
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a) Régime libre conservatif 

Un régime conservatif, caractérisé par un amortissement c = 0 ne dissipe 
pas d'énergie par amortissement. La solution en déplacement de l'équation 
différentielle est : 

Ùq 

u(t) = — sin o)(t) + uo cos (o(t) 

Ci) 

Nous pouvons composer ces deux fonctions en les projetant sur l'axe des 
réels afin d'obtenir un vecteur résultat de longueur p — \ju^ + ^— ^ et 
d'angle de phase 6 = tan -1 (zi 0 /^o^X ainsi : 

u(t) = pcos(cot — 6) 

b) Régime libre dissipatif 

Pour c 7^ 0, la solution de l'équation caractéritisque est : 



S = i;ù) dz <j}\J g 2 — 1 

En fonction de la valeur donnée à l'amortissement =, > ou < à 1) , 
trois types de comportement peuvent être distingués : 

• Amortissement critique, i; = 1 
La solution, apériodique, est : 

u(t) = (mo(1 + o)t) + ùot) exp~ Mt 

• Amortissement surcritique, è, > 1 



en notant co = coy (g 2 — l), la solution apériodique, est : 

k(0 = cosh(&0 + ^° S inh(û>0^ e - ^ 



10.1 Système à 1 degré de liberté 
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• Amortissement souscritique, £ < 1 

C'est le cas le plus courant rencontré en dynamique des structures. La 
pseudo-pulsation de l'oscilateur amorti co D est égale à co D = cû^J (l — £ 2 ). 
L'amortissement a donc tendance à diminuer les pulsations propres d'un 
système. Le déplacement s'exprime comme suit : 

u(t) = \uq cos(co D t) + ^ 0)110 + ^° sm(û) D t)] e - ^ 

v y 

Ces deux fonctions trigonométriques périodiques peuvent se combiner pour 
exprimer le déplacement sous la forme simple suivante : 

u(t) = pcos(co D t - d)z~ èo)t 



, 7 /^Wo + WoV ^ + _ 1 {Ç(DU 0 + Ù( ) 

Avec p = \ ui + et ^ = tan 

( ' ù) D J V (Wd^O 



10.1.3 Le régime forcé sinusoïdal 

Dans ce cas de figure, nous allons étudier la réponse du système à un degré 
de liberté soumis à une excitation harmonique du type sinusoïdal. En effet, 
nous verrons par la suite que les systèmes les plus complexes, soumis à 
des charges les plus variées peuvent se ramener à ce cas précis, le second 
membre de l' équations différentielle du mouvement s'exprime maintenant 
par p(t) = po sin(ù)t). Nous noterons fi, le rapport de la pulsation excitatrice 
à la pulsation propre du système : /3 = cb/co. 

a) Régime forcé conservatif 

La solution complète s'écrit : 



u(t) = — 



1 

1 -f3 2 



( sin(ù)t) — fi sin(ù)t)) 
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On définit le facteur de réponse dynamique comme le rapport entre la 
réponse dynamique et le déplacement statique (u stat = po/k). Ici : 



À la résonance, (3 = 1 et R(t) — > oo. 
b) Le régime forcé dissipatif 

La résolution de l'équation du mouvement nous fournit un déplacement 
ayant pour forme la somme d'un terme en régime transitoire et d'un terme 
en régime permanent forcé. Si on ne s'intéresse qu'à la réponse en régime 
permanent après disparition des termes transitoires par amortissement, le 
déplacement du système s'écrit : 



R(t) 



u(t) 1 



( sin(o)t) — fi sin((ot)) 



u stat 1 — fi 2 




1 



2 ((1 - 0 1 ) sin(ûtt) - cos(û>0) 



u(t) = 



£ 2 ) 2 +(2^) 2 




o 



o 



0,5 



1,5 
(3 = cb/ù) 



2 



2,5 



3 



Figure 10.2 



10.1 Système à 1 degré de liberté 
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Le facteur de réponse dynamique maximal, noté R mâx vaut dans ce cas : 

1 



^max 



10.1.4 Régime permanent sous une charge périodique 
quelconque 

a) Réponse en fréquence 

Une méthode de calcul plus efficace consiste à remplacer les efforts et dépla- 
cements par des nombres complexes. En régime forcé harmonique, on peut 
exprimer l'effort extérieur p(t) en fonction d'une force complexe p(t) : 

PO) = Poe lû)t = poicoswt + / sincôt) 

Dans ce cas, la partie réelle ou la partie imaginaire du chargement pourra 
être utilisée. Le déplacement u(t) est alors la partie réelle ou imaginaire du 
déplacement complexe û(t), ainsi, 

û(t) = Ge im 

En appelant Z l'impédance complexe et H l'admittance complexe, nous 
définissons la réponse complexe en fréquence comme : 

1 1 1 

H(ô)) = -, = - 7 

Z(û>) k(l-/3 2 )+2i£/3 

Le déplacement (complexe) peut donc s'écrire : 

û(t) = H(â))poe iibt = \H(ô))\e m - e) 

avec : 



et 

tan 6 = — hJ --, 
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b) Représentation d'une charge périodique en série complexe 
de Fourier 

On peut représenter une charge périodique, de période 7> par décomposition 
en série de Fourier de la forme : 

oo oo 

p(t) = ci{) + ^ a n cos ô) n t + ^ b n sin co n t 

n—l n—\ 

La charge est ainsi décomposée en une somme d'une charge constante 
représentant la moyenne et d'une série de charges harmoniques de fré- 
quence â) n ^ d'amplitude a n et b n , représentant la variation par rapport à la 
valeur moyenne. Par les équations d'Euler, nous pouvons relier ces fonc- 
tions trigonométriques aux fonctions exponentielles complexes. La série de 
Fourier pour la charge p(t) devient : 

oo 

p(t) = Y J Pne i ^ 

— OO 

Avec : 

Pn = ^r I ' piOe-^dt 
l p Jo 

c) Réponse à un chargement décomposé en série complexe 
de Fourier 

Du principe de superposition, la réponse totale permanente d'un oscillateur 
élémentaire à une excitation périodique peut s'écrire : 

oo 

u{t) = Y, H{îb n )P n {(b n )é ii '" t 

— OO 

10.1.5 Réponse à une charge arbitraire 

Dans ce parapraphe, nous décrivons deux manières différentes (temporelle 
et fréquentielle) permettant d'obtenir la réponse d'un système à un ddl sou- 
mis à une action dynamique quelconque. 



10.1 Système à 1 degré de liberté 
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a) Réponse simplifiée à une impulsion 

Considérons une charge impulsionnelle p(t). Si la durée du chargement 
(t\) est très courte par rapport à la période de vibration du système 
T (t\ <C T), on peut faire l'hypothèse qu'il n'y a pas de changement notable 
du déplacement durant le temps d'application du chargement mais par 
contre qu'un changement de vitesse peut être considéré. Soit la relation 
entre impulsion et quantité de mouvement : 

/ = mM ~ / p(t)dt 
Jo 

Si l'impulsion (unitaire de Dirac) a lieu à t = r, le déplacement du sys- 
tème au temps t = t — r se met sous la forme : 

u(t) — — !— sin col = h(t) 
inco 

oùh(t) est la réponse à une impulsion unitaire. 
Pour un système amorti : 

h(t) = —^— e ~ €œt sin co D t 
mcoo 



b) Réponse temporelle à un chargement arbitraire 

Le chargement est ici considéré comme une suite d'impulsions unitaires. La 
réponse totale est obtenue par la somme des réponses différentielles : 



/' 

Jo 



u(t) = / p(r)h(t - r)dr 



soit : 

1 f f 

u(t) = / p(r) sin co(t — r)dr 

mco J 0 

Cette dernière équation est appelée intégrale de Duhamel. 
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c) Réponse en fréquence à un chargement arbitraire 

Il s'agit d'une extension aux fonctions non périodiques de la décomposition 
en séries de Fourier. Une charge quelconque (ou fonction quelconque, p(t)) 
peut être décrite par une intégrale représentant les contributions des com- 
posantes harmoniques ayant un spectre de fréquence continu de — oo à +00. 
Les intégrales : 

/+00 t p+00 

p(t)e~ iG)t dt et p(t) = — / P{(b)e iM d(b 
-00 ^ j — 00 

constituent une paire de transformée de Fourier. P(cô) est appelée transfor- 
mée de Fourier de p(t) et p(t) est appelée transformée inverse de Fourier 
deP(û>). 

Par analogie à la réponse d'un système soumis à une charge périodique 
décomposée en séries complexes de Fourier, la réponse d'un système à un 
ddl est obtenue en calculant la somme des composantes de la réponse sur 
l'étendue des fréquences, on obtient ainsi : 

+00 1 p+00 



u(t) = —- I U{(b)e mt dâ) = -3- / H(cb)P(ù))e i5)t dù) 



-[ 

/+OO 
U 
-00 



avec : 

Û(ô)) = j u{t)e~ l(0t dt = H(ù))P(cb) 



Remarque : La réponse impulsionnelle h(t) est la transformée 
inverse de Fourier de la réponse en fréquence H(â)). 



10.1.6 Réponse à des chargements impulsionnels simples 

La tableau ci-après fournit le coefficient d'amplification dynamique maxi- 
mal obtenu pour un système à un degré de liberté soumis à des charges 
impulsionnels sinusoïdale et rectangulaire. Pour des impulsions de courtes 
durées, nous considérons ci-après que le temps d'application de l'effort 
t\ est petit devant la période propre de l'oscillateur T. Ainsi, le déplace- 
ment maximum est atteint pendant la vibration libre, à la fin du chargement 
(/3 = ùj/ùj = T/T > 1). 



10.2 Système à N degrés de liberté 
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Description 


schémas 


/^max 


impulsion sinusoïdale 




< 

l \-+ t 


2/3 TT 

— ^ cos 

1 - /3 2 2/3 


impulsion rectangulaire 


P(t) À 


< 

1 P0 

' >t 


2sin TT — 
T 



10.2 SYSTÈME À N DEGRÉS DE LIBERTÉ 

10.2.1 Équations du mouvement 

Après discrétisation du problème mécanique, de manière analogue par 
exemple à la méthode des éléments finis en statique, nous pouvons 
construire les matrices de masse M, d'amortissement visqueux C, de 
raideur K et le vecteur des efforts extérieurs F. L'équation du mouvement 
sous forme matricielle peut se mettre sous la forme : 

MU + CÙ + KU F 

U est le vecteur des déplacements des degrés de liberté considérés. 

10.2.2 Signification des modes propres et fréquences 
propres 

Nous supposerons dans cette partie un système discret, conservatif (C = 0) 
et sans chargement extérieur. L'équation discrète du mouvement prendra la 
forme suivante : 

MÛ +KU =0 
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Prenons l'exemple d'un bâtiment de trois étages. Concentrons la masse 
au niveau des planchers, nous obtenons ainsi un système à trois ddl dyna- 
miques. Soumis à un chargement quelconque, le mouvement horizontal de 
chaque plancher n'est plus harmonique. Par contre, si la vibration libre de la 
structure est initiée par des vitesses et des déplacements appropriés imposés 
aux masses, la structure oscillera suivant un mouvement harmonique simple, 
trois allures caractéristiques de la déformée existent pour un système à trois 
ddl. Ce sont les modes propres de vibration. 




mode 1 mode 2 

Figure 10.3 



WXXW WXXW XXtXW XXXXW 



mode 3 



À chaque mode nous pouvons associer une fréquence propre de vibra- 
tion : 

f _ 1 co n 

Jn — rj 



27T 



Avec : coi < û)2 < C03 < ù) N 



Le vecteur <j> n définit l'allure de la déformée de la structure répondant à 
la pulsation propre co n . <p n définit seulement le rapport des déplacements des 
planchers entre eux (u\/u2,u\/u?). 



10.2.3 Détermination des fréquences propres de vibration 

La réponse d'un mode propre est harmonique, ainsi le vecteur des déplace- 
ments peut s'écrire : 

U(t) = Û sin(û)(0 + 0) 



10.2 Système à N degrés de liberté 
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Û représente les modes de déformations possibles qui ne changent pas avec 
le temps (seule l'amplitude varie). L'analyse du système nous montre qu'une 
solution n'est possible que si : 

det(M- l K-co 2 l) = 0 

C'est l'équation aux fréquences du système. En développant le déterminant, 
nous obtenons une équation polynomiale de degré N en o? . Les N racines 
sont les fréquences propres des N modes de vibration du système (o>\, o? 2 , 

2 2 \ 

Ù) 3 ,...Ù) N ). 

10.2.4 Détermination des modes propres de vibration 

Pour compléter le problème aux valeurs propres, on peut calculer les N 
formes de vibrations correspondant aux N modes de vibration. Si nous rem- 
plaçons les valeurs de co n : 

[K-û^M](J7 n ) = 0 

L'ensemble des vecteurs modaux de vibration propre peuvent ête mise sous 
forme matricielle : 

[(/)] = [ c f) l d> 2 </>3 c/>n] 

10.2.5 Propriété d'orthogonalité des modes 

Si (o n ^ (o m alors, 

et 

0mM0„ = 0 

Les vecteurs modes de vibration sont dits orthogonaux par rapport aux 
matrices M et K. Ainsi la matrice [</> r ]M[</>] = [M] est diagonale et 

cf> T n Mcf> n = M n 
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est appelée masse généralisée du mode n. De façon analogue pour la matrice 
de raideur, nous pouvons définir la matrice de raideur généralisée [K] par : 

cf> T n Kcf> n = K n 

10.2.6 Normalisation des vecteurs modes de vibration 

Il est avantageux de normaliser U n de façon à obtenir : 

Vn^ n = 1 
1 T 

Ainsi, si <p n = — 777 C/ n alors [<fi ]M[<fi] = I. 

10.2.7 Équations modales du mouvement - Superposition 
des modes 

a) Coordonnées normales 

Les modes constituent un ensemble de coordonnées qui ne sont pas couplées 
dynamiquement. Pour un système à N ddl, ils constituent N allures indépen- 
dantes. Leur amplitude sert de coordoonnées généralisées. Tout vecteur U (t) 
peut être trouvé par superposition des amplitudes des modes de vibration. 

U(t) = faYtf) + 02*2(0 + 03*3(0 + .... + cf> N Y N (t) 

N 

U(t) = Y J <t>nY n {t) = WY(t) 

n=\ 

Ainsi les coordonnées normales appréciant l'amplitude de chacun des 
modes peuvent être déterminées : 

_ <f> T n MU(t) 
Yn[t) - cf>T n M<f> n 



10.2 Système à N degrés de liberté 
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b) Équations découplées du mouvement non amorti 

Soit un système à N ddl : 

MÙ + KU = F(t) 

Dans le domaine linéaire, [<fi] reste constant au cours du temps. Par pro- 
jection dans l'espace des modes propres, l'équation du mouvement peut se 
mettre sous la forme : 

4>lM4> n Y n + <ftK4> n Y n = cf> T n F(t) 

Appelons M n — <filM<fi n la masse généralisée, K n = <filK<fi n la raideur 
généralisée et F n = <f>%F(t) la force généralisée. L'équation du mouvement 
se retrouve sous sa forme scalaire en Y n (t) pour le mode n : 

M n Y n + K n Y n = F n 

Soit encore : 

On résoud les Y n pour n = 1,2,....,^. Connaissant Y(t), on calcule la 
réponse dans le système de coordonnées géométriques : 



U(t) = Y t 4> n Y n (t) = [4>]Y(t) 



n=\ 



c) Équations découplées du mouvement amorti 

Si les 4> n ne sont pas orthogonaux par rapport à la matrice d'amortissement 
C alors on ne peut pas découpler les équations du mouvement. Dans le cas 
contraire, l'équation discrète du mouvement peut se mettre sous la forme 
suivante pour le mode n : 



M n Y n + C n Y n + K n Y n — F n 
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avec : 

tâCfin = C n 

et : 

tâCfa = 0 

Si de plus nous introduisons l'amortissement relatif C n — 2g n co n M n alors 
nous obtenons n équations découplées en Y n : 

Y n +2Ç n a>J n + (olY n = ^- 

Il est toutefois plus facile et plus physique de définir l'amortissement relatif 
% n par mode que le coefficient C n . 

d) Condition d'orthogonalité par rapport à C 

Lord Rayleigh a montré qu'une matrice C proportionnelle à la matrice de 
masse et/ou de raideur satisfait les conditions d'orthogonalité : 

C = aM + ySK 

Connaissant la valeur de l'amortissement à imposer à deux modes m etn, 
nous pouvons calculer les deux paramètres a et fi par l'équation suivante : 

1 1 

ù) n ù) m _ 

Il existe un nombre infini de ce type de matrice, définie par la série de Cau- 
ghey : 

oo 

C = M a b (M- l K) b 

b= — oo 



a 



10.3 Vibration des systèmes continus 
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10.3 VIBRATION DES SYSTEMES CONTINUS 
10.3.1 Vibration axiale des barres 

Considérons une barre d'axe x, de longueur L, de section A(x), de densité 
linéique p, composée d'un matériau élastique linéaire de module young E 
soumise à une densité linéique d'effort q(x, t). 





0 








h * 



Figure 10.4 

Quelle que soit l'abscisse x, l'équation générale de vibration axiale de la 
barre est 



d_ 
dx 



E.A(x) 



du 
dx 



d 2 u 

pA(x)—+q(x,t) = 0 



Les conditions limites, en effort ou en déplacement sont données ci-après : 
du 

en x = 0, E.A(x)— (0, t) + q(0,t) = 0 ou u(0,t) = 0 et en x = L, 
dx 

du 

E.A(x)—(L, t) + q(L, t) = 0 ou m(L, t) = 0. 
dx 

Pour une barre de section constante A, en vibration libre (sans effort exté- 
rieure q(x,t) = 0) l'équation dynamique se ramène à l'équation des ondes : 



d 2 u 1 d 2 u 

dx 2 = 



avec c = y E/p, la célérité des ondes uniaxiales de traction-compression 
dans le milieu élastique. 
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10.3.2 Vibration transversale des poutres 

Soit une poutre isostatique d'axe x, de longueur L, de moment quadratique 
constant /, de densité linéique p, composée d'un matériau élastique linéaire 
de module young E. Les pulsations propres de flexion lorsque la poutre 
repose sur deux appuis simples sont : 



10.3.3 Détermination du mode fondamental de vibration : 
méthode de Rayleigh 

Le but de la méthode de Rayleigh, par application du principe de conser- 
vation de l'énergie, consiste à déterminer le mode fondamental de vibration 
d'une structure. Il faut être capable d'estimer la forme du mode propre 
recherché. Pour cela, en considérant qu'un mode de vibration correspond 
à un état de déformation de la structure pour lequel les efforts de rappel 
élastiques équilibrent les efforts d'inertie, la déformée du mode propre est 
calculée par une analyse statique, en considérant un chargement inertiel 
appliqué à l'ensemble des masses de la structure. Considérons la structure 
simple ci-dessous, de masse linéique variable m(x), dont on cherche à 
estimer la première pulsation propre de vibration horitonzale. Le char- 
gement équivalent permettant de calculer la déformée modale est obtenu 
en appliquant une accélération uniforme à l'ensemble de la structure. La 
déformée obtenue est notée if/(x). 




Les n modes propres associés sont : 




10.3 Vibration des systèmes continus 
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Îx 



F (x) = m(x) i 



r\ \ \ \ \ 

Chargement équivalent Déformée 
Figure 10.5 



En égalant les énergies cinétiques et potentielles maximum calculées sur 
la base d'un mouvement harmonique, la pulsation fondamentale de la struc- 
ture se ramène à : 



/ EI(x)[ip // (x)] 2 dx 
Jo 



m(x)i// (x)dx 



Pour un système discret, le quotient de Rayleigh va s'exprimer comme suit 
permettant de déterminer le mode n de vibration : 



10.3.4 Modes propres de vibration des poutres 

Le premier mode propre de vibration est donné dans le tableau ci-après pour 
des poutres de longueur /, de module d' Young E, de masse linéique m et de 
moment quadratique /, supportant différents types de conditions limites. La 
solution générale est de la forme : 
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Type de liaisons 


schémas 


a 


Appui simple - Rotule 




9.87 


Appui simple - Encastrement 


WWW 




15.4 


Encastrement - Encastrement 


www 




22.4 


Encastrement - Libre 


/ 
/ 

/ 
y 
y 




3.51 



10.3.5 Modes propres de vibration des plaques 

Le premier mode propre de vibration est donné dans le tableau ci-après pour 
des plaques d'épaisseur constante h, de masse volumique p, de grand côté 
a et de petit côté b, supportant différents types de conditions limites. La 
solution générale est de la forme : 

V ma 4 



Type 
de liaisons 


schémas 


b/a 
= 1 


b/a 
= 1.5 


b/a 
= 2 


b/a 
= 2.5 


b/a 
= 3 


oo 


Appui simple 


«r t «r <r 
X X X X 


5.7 


4.17 


3.56 


3.3 


3.17 


2.85 


Appui simple - 
Encastrement 


T T T T. 

1 4 J 

X X À. X ' 


6.83 


5.46 


5.0 


4.8 


4.69 


4.45 


Encastrement - 
Encastrement 


\ f 


8.36 


7.23 


6.88 


6.71 


6.64 


6.46 


Encastrement - 
Encastrement 




10.38 


7.79 


7.09 


6.86 


6.69 


6.46 



Remarque : Pour une plaque circulaire, de diamètre a, simplement 
appuyé sur son pourtour, j3 = 1 1 .84 et encastrée sur son bord (5 = 5.9. 
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